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HE SOLUT IONS of general canonical systems show a 
very complex behavior; various aspects of the subject 
have been considered in the literature (1-4).2. The nature of 
the first integral sometimes serves to indicate some proper- 
ties of the solution; for example, if every line U(a,y) = 
constant is bounded, then every solution of the system 


dx _ OU(x,y) dy _ _ OU (x,y) 
dt by dl Ox 


is also bounded. Atkinson (5) has examined the equation 
+ =0 [*] 


which is readily given the canonical form for f(t) = 0; then 
the first integral is <2 + 2?"/n. Atkinson has formulated the 
conditions on f(t) sufficient to give the function #°(t) + 
x*"(t)/n a finite limit as t ~ o, if z(t) is any solution of 
Eq. [*]. 

The method adopted here is applicable to any system of 2n 
equations very similar to a canonical one. In particular, 
Theorem 2a states restrictions on the right sides of the equa- 
tions 


— 
Ox; pA fe; 


sufficient to make all solutions of this system bounded as t > 
o. Theorem 1 states the condition sufficient to make 


lim H(x(t),y(), = 
to 


for each bounded solution? of 


dary _ OH (x,y,t) t) 

dt oy; + gi(z,y,l) hu(z,y,t) 

d _ oH (z,y,t 
= ) + gio(x,y,t) + his(2,y,t) 


Translated from Vestnik agro Universiteta (Bulletin of 
Moscow University, Series 1, Mathematics-Mechanics), 1960, no. 
3 Pp. 26-36. Translated by Research Information Service, New 

or 

‘Some of these theorems have been presented previously 
leg. ref. (10)], but without proof. 

? Numbers in parentheses indicate References at end of paper. 

8 Here and elsewhere the ig a er of the functions are given 
in condensed form; H(z, y, t) denotes H(m, ..., 2n, Ya, 

Yn, t), and 2(t), y(t) denote .. &n(t) and y,(t), 
and so on. 
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This is a generalization of Ascoli’s (6) theorem. 

In Section 2, some conditions are given sufficient to insure 
that ¢ + f(z,#,t) = 0 has solutions that resemble the solutions 
of the corresponding canonical equation, namely ones that 
tend to constant limits ast ~ ©. In Section 2, the equatior 
% + x = f(z,2,t) is used to demonstrate the results given by 
the theorems in this simple case. 


=> N 


Consider the real range of the system of Eq. 1; all fune- 
tions appearing on the right are defined in the space XY YT, 
which is the product of the space YY of the variables 2, 
, and y1, ..., Yn by the half-axist > T. It is assumed 
that H(z,y,t) and the solutions and 
yn(t) are continuous. It is possible that system 1 has a 
definite and continuous solution provided that T < t <7 < 
©, but the solution cannot be extended continuously beyond 
the point 7. This is an inextensible solution; the equation 
& + 2 + t-y?2 = 0 has an infinite number of such solutions, 
for example. 
The method presented here is applicable to any solution 
_ without change in the form given to the theorems, so in the 
~ following solutions to Eq. 1 will be dealt with that are defined 
and continuous for T <7 < ~, and also with the asymp- 


totic behavior fort—~>7. 


Let R be some region in XY, and in RT let the functions 
H, 0H/d2;, OH /dy; = 1, ..., n) be bounded (as regards 
modulus) by a constant C; further, let 


oH 
(24 ha + 5 is) 


i=1 


Theorem 1 


Finally, let there be a function ¢(t) 


<7 we have 


be of constant sign. 
such that 


and that in the range T < 


| H'(2,y,t) | (| gu(x,y,t)| | gio(a,y,t) | ) [2] 
i=1 
forallzand yin R. Then there exists a finite limit { 


lim H(x(é),y(0),t) 


for each solution x(t), y(t) of Eq. 1 for R in the range T : 
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Proof 


Let us suppose that 


"(oH 


H 
hi ha) <0 


The conditions imposed imply that 


i=1 
H(x(t),y(O,0)| 
T < t <7, so the set of H(2( 
point for t > 7. 
tvo such points; then the continuity in H, x(t) and y(t) im- 
plies that there are numbers @ and 6B (8 > a) and a mono- 
tonic sequence t — 7 such that for all k 


(xbox 41), = [3] 


The total derivative dH/dt along x(t),y(t) in the intervals 
(/y,t¢41) is found by taking 


(t),y(t),t) 


E(x (tog) (lor) =a 


d/1(a t),t oH oH 


i=1 


< for 
has at least one limit 
Let us suppose that there are not less than je 5 i (= 
U(x,y) 


oH oH OH (OH of 
(2 ha + oy; + ot =. ( + gi 


Lemma | 


We are given a set R C XY, 


a solution z(t), y(t) and an 
infinite monotonically increasing sequence i, ~ 7 (T< &) 


such that the points (x(é),y(t)) lie in R for ty < 
Conditions that are obeyed are 


| U(x,y) e< [5] 
ha + = hs) 2 0 [6] 
Ovi 


in which ¢ takes the value +1 or —1; let there be a function 
g(t) € L;,,such that for (x,y) € Rand T < t <7 we have 


t < 
Ox; 


oH 
| We(a,y,t) | >> (2 Jia +o 


+H < 
| 


Ox; OY; 


i= 


) + |H’,) < 2nCe(t) + = (2nC + 1)¢(t) [4] 


which is to be integrated the limits and in 
conjunction with Iq. 3 to give 


B-a< (2nC + 1) ¢(t)at 


in which 8 — a@ is a positive constant. The integral is a 
function of k and decreases without bound as & increases; 
in fact, g(t) © Lm., when t, > 7, so Eq. 4 is contradictory if 
k is large. This means that we were wrong to suppose that 


H has several limiting points. Q.F.D. 
Further, Eq. 2 implies that lim H(z,y,t) exists for each 


t—r 


point (z,y) € R; let that limit be U(z,y). 
true, then all points Q of the limiting set of the loci x(é) 
lie on a surface U(zx,y) = constant. 

In what follows it shall be supposed that U(x,y) is continu- 


If Theorem 1 is 
» 


cannot have simultaneously for all & 


[8] 


assert that we 


U (lox) = (tex +1) 


Then we 


ifq > (1 + c)/(1 — isa positive number. 
Proof 
We produce along x(t), y(t) the derivative 
1 dH 1 oH oH 
( Ji + 


and integrate the two sides of Eq. 9 from ty to tx41; let € 
equal +1, and suppose that Eq. 8 applies to that interval. 


s; the difference H(x,y,t) — U(a, ill be denoted by : : ; 1 
= ee Sen el Then we have for the left-hand side that 
L = In = In < = 
+y (tox) (tox) / U (tex) 
q 


and for the right that 


L = 
tons jl (oH oH (oH oH 
1 v U VU (24 Gil + + (3 + + U ( dt 


Now we estimate a lower limit to this expression. 
sum in the integral is not negative (« = +1), so 


2 dt > 
Si 1+ ¥/U + dy 92) + uf dt > 
1 


i=1 


vix,y,t). We shall now prove several theorems analogous to 
Theorem 1, except in that no bounds are imposed on H, 
0/1 /dx;and 0H/dy;. These theorems can be used to examine 
the unbounded solutions of Eq. 1 that tend to zero. The 
mode of proof is the same in all cases, which are formulated 
here as a lemma. The notation which shall be used is 
U‘a(t),y(t)) = Ut), = and so on. 
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The condition implies that 1/(1 + y’//U) > 1/1 — ce) > 0; 


the second 


(The last inequality follows from Eqs. 5 and 7.) Now ¢(t) 
is integrable, so the right part becomes negative for t; — 7, 
but the modulus can be made as small as may be desired 
when k is large. This estimate of L conflicts with that given 
by Eq. 10, so Eq. 8 is contradictory for k large. 

In what follows we shall examine the solution 2(¢),y(é) to 
Eq. 1, which we suppose to be defined for T < t < 7 and to 
lie in some manifold K C XY. The case K = XY corre- 
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sponds to an arbitrary solution to Eq. 1; the case AK c XY 
enables us to examine at least some of the solutions to Eq. 
1, if the restrictions on the complement are not complied 
with throughout the whole of XY. 

The following theorems are presented in condensed form, 
and similar proofs are omitted, with a view to economy in 
space. 


Theorem 2 


We define a set Rei{Rc:] of points (x,y) by means of the 
relations (ry) € K, U(z,y) > C [<C]. Let there be a C 
such that Eqs. 5, 6 and 7 apply to this set, and let ¢ 
equal —1. Then 


(1 + c) 
2) 


(in which ¢ is defined by Eq. 5) and use the fact that U(¢) is 
continuous to select a new monotonically increasing sequence 
> 7, such that 


U (tex 1) = qU (tex) [11] 
(tex) < U(t) < U (tex 41) for toy. < t < [12] 


12 shows that (2(8),y(t)) Rei for ty < t < so 
we have the conditions of the lemma. But « = —1, so Eq. 
11 contradicts the assertion in the lemma. The theorem is 
thus proved. A simple consequence is as follows. 


Theorem 2a 
We have the system 


_ 


Now Eq. 


n 
Let V(x) be continuous; let 7 a 


Vay 


i=1 


. [14] 


Vis) = 
sy) = 


= 


and let us suppose that outside a certain sphere S within the 
space of x and y we have for t > T that : 


s=1 


+) Yo < ¢(t- 


U(xy) 


[15] 


Then all extensible solutions of Eq. 13, together with their 
derivatives, are bounded for t 

Note: If the right sides of Eq. 13 are continuous and obey 
the conditions, the system has no unbounded solutions; the 
absence of unbounded inextensible solutions is demonstrated 
as for unbounded extensible ones (Eq. 15 is applied to a finite 
range in 


4 We prove the first assertion in the theorem. 
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< const [> const] for T < t<7 
Proof 
We assume the converse, i.e., that U(te) ~ + © along 
some sequence fg ~ 7.4. We take some number 2 


Proof of Theorem 2a 

The conditions are those of Theorem 2. We introduce 
variables y; = #;; then Eq. 13 takes the form of Eq. 1, with 
oy = ha = gin = 0, and 


= H(xyt) = Via) + + 


+ 


Then the Re, of Theorem 2 is found by putting K = VY 
and using the fact that the set of points (x,y), as defined by 
V(r) + yi? > const 
t= 


for values of that constant sufficiently large, lies outside th: at 


sphere. The assertion of Theorem 2 amounts to the assertion 
that along any solution 

l n 

V(r) + 5 < const [16 

so V(2(t)) < constant, which in conjunction with condition 
14 gives® for t ~ +0 a 
n 


i=1 


The continuity of V(x), in conjunction with Eq. 14, implies 
that V(x) > constant for alla. This, with Eq. 16, gives us 
n 


< const for 1 + 


i=1 


em 3 


Let conditions 2, 6 and 7 be complied with for some C in 
the set Ro and let be +1. Then if 


= + flim... = 


lim U(a(0,y(0) 


lim U( 


the more exact relation 


is also obey ed. 

Now several theorems will be given that can be used to 
examine the solutions to Eq. 1 near the equilibrium point of 
the abbreviated canonical equation. The proofs of Theorems 
3 to 5 are omitted because they are merely simple deductions 
from the lemma. 


We define a set Res[Rea] by the relations 
(vy) ERK (0 > U( 


Let conditions 5, 6 and 7 be complied with for some C in the 
set, and let e be +1. Then if 


Theorem 4 


U(z,yy) < C 


(lim 


lim U(x(),y(O) < 0 (r(),y() > 0] 


5 Condition 14 is used only at this point in the proof. The 
solutions to Eq. 13 can be shown to be bounded when less rigorous 
conditions are imposed on V(x), for example 


lim min Vie) = +o 

xc? = 
1 
which takes the form in Vir) = + ~ forn = 1; it is the 

necessary and suffice baat condition for all solutions of the abbre- 
viated canonical equation # = —V'(x) to be bounded. 


ARS JouRNAL SUPPLEMENT 


> 
‘ 
| 
( 
« 
J 
ul 
the 
|| 
a. 


we have for all ¢ sufficiently close to 7 that 


> 0) 


U(a(t),y(t)) 


Theorem 5 


Let conditions 5, 6 and 7 be satisfied for some C in the set, 
and let ebe —1; the fact that x(t), y(t) obeys 


lim U(r(t),y()) < 0 flim U(r(t).y(Q) > 0] 


sr 


implies that a more exact relation 


lim U(x(t),y(d) lim > 0] 

is also obeyed. 

Theorem 4a 
We have system 13 


ind we suppose that in a region near the point 0 in the xy 
plane the system obeys, for ¢ > 7, the conditions that V(x) 
is continuous, that (0) = 0, that > 0 for ~ 0), and 


n 
= 
j We + < U(ay) 
i=1 
Vir) < Ury)! 


1 


gl(di) < « 


Then Eq. 13 has no solution that tends to zero along with its 
derivatives as ©, 


Proof 


Any region about the point (0,0) € NY contains at least 
one of the regions defined by the relations 0 Uiey) < C. 
This means that set Re. of Theorem 4 exists, so Theorem 4 
ipplies to Eq. 13. Then if a solution to Eq. 13 exists that 
tends to zero as t + ©, for é > T that solution satisfies 


Vat <0 


Now V(x) > 0, so it is necessary to have V(x) = 0 and #;? = 
0% But V (x) is positive definite, so.v = = 0. 


Given here are theorems on the existence of solutions that 
‘end to constant values for a nonlinear second-order equation ; 
this is an extension of Villari’s (7) theorem. The proof is 
mitted because it follows the lines of Villari’s proof exactly. 

We call the function f(x,y,f) g-continuous in the region S 
f the variables x and y if for any ¢ > 0 we can find a 6 > 0 
uch that if 

(rin) ES 


| Lal 


(2,2) 
<4 


then 
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for all ¢ within the region in which f exists. We define a 


Lipschitz g-condition in a similar way 


Let f(c.0,.) 4 0, and let 
> Ta function g(t) such that 


l. > 0, 


2. {f(c,0,0) 


r, We are given for 1 > T that 
t+ = 


together with a real number c. 
there be for ¢ 


ty(tdt << 


€ ¢(b) 


3. Ina certain region S around the point (c, 0) in the VY) 
plane the function f(.c,y,0) is g-continuous. 
Then Eq. 17 has a solution whose properties are 


wer 


For any O there is a such that 


| | < (1 + €) s — t)y(s)ds 


That solution is unique if f satisfies a Lipschitz g-condi- 
tion in S. Villari had 


= A(t) - f(x) 


< ise 


pr +2 


Winter and Hartman (8) have shown that if for every 
region S around (¢,0) on the YY plane we have 


max < —¢s(t) < 0 tes(t)dt = + 
(x,y) ES 


then Hq. has no solution whose asymptote is Hq. 19. 


Theorem 6 is, in this sense, correct. 


= g(xyy,b) [20] 


Theorem 7 
We are given for > 7 that 
= f(a,y,t) 


and numbers c; and such that + 0. 
Let us suppose that we can select a function g(t) such that 


2. | f(er,¢2,t)| + < 


3. f(ayt) and are g-continuous in some region S 
around the point (¢1,¢2). 

Then Eq. 20 has a solution whose asymptotic behavior is 
This solution is unique if f 
and g satisfy a Lipschitz g-condition. This may be proved 
by means of Villari’s method of successive approximations. 

Villari (9) has given similar theorems for nth order equa- 
tions; these may be reformulated in terms of integral esti- 


mates for the right sides. See 


> y(t) co fort > 


§3 


In this section we shall discuss application to the equation 


E+ = [21] 
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The foregoing theorems, applied to Eq. 21, give rise to the 
following conclusions. Suppose that g(t) € LT 


if 
\f(xy.t)| < + [22] 


for some region S in the XY plane fort > 77, then all extensible 
solutions of Eq. 21 are bounded as t > ~. 

2. If f(z,y,t) is also continuous, then all solutions of Eq. 
21 are extensible. 

3. If for every R > 0 we can select a gp(t) € Ly» such 
that 2? + y? < Rif(z,y,t)| < ¢r(t), then any bounded solu- 


tion of Eq. 21 has the property that 
constant. 
4. If a and yw are such that 


lim {2x2(t) + y%(t)} = 
ta 


f(x cos t + yo sin t, — xo sin t + yo cos t,t) =O 
then Eq. 21 has as solution 


x(t) = x cost + ysint 
= —2 sin t + yo cost 
5. If a and y are such that ao? + yo? ¥ 0 
f(%o cost + yo sin t, — ao sin t + yo cos t,t) 4 0 


and that in some region around the point (20,yo) this f is ¢- 
continuous and g-bounded, then Eq. 21 has at least one 


solution whose asymptotic behavior is 7 
x(t) cost + yo sini + 


ry i(t) = —2 sin t + yo cos t + O(1) 


> & 


“a. 

6. If the condition of g-continuity in 5 is replaced by 
Lipschitz ¢-condition, we may say that one solution having 
the form of Eqs. 23 corresponds to the point (x,y). 

7. If Eq. 22 applies to some region around the point 
(0,0) in the Y ¥ plane, then Eq. 21 has no solution that tends 
to zero. 

Statements 1, 2 and 3 are direct applications of Theorems 
2a, 1 and 4a; statement 4 is obvious. We may prove state- 
ments 5 and 6 by performing a substitution in Eq. 21, namely 

x= £&cost+ nsint 
—ésint+ neost 


This substitution is mutually continuous throughout the 
space YY 7; then Eq. 21 becomes 


— sin t-f(E cost + nsint, — Esint + 7 cos t,t) 


Theorem 7 is applicable to Eqs. 24 if (20,0) satisfies the 
requirements listed in statement 5. This gives us the solu- 
tion 


x(t) = [x + O(1)] cost + [yo + O(1)] sin t = x cost + 
yo sin t + O(1) 


Now we assume that there is a solution to Eq. 21 of the 


form 


x(t) = cost + yo sint + = &() cost + n(t) sin t 
z(t) = —x sin t + yo cost + = —é(t) sint + cost 


We solve these identities for &(t) and n(t) to get 


E(t) = a + 0:(1) cos t — 02(1) sin t 


nt) = yo + 0:(1) sin + 0.(1) cos 


23 cor- 


Thus two different solutions having the form of Eqs. 
respond to two different solutions of Eq. 2 that tend to a 

single point (.r0,yo) [if 0:(1) cos t — 0.(1) sin t = 0,*(1) cos t — 

—0.*(1) sin ¢t and 0,*(1) sin ¢ + 0.*(1) cos t = 0,(1) sin t + 
0.(1) cos t, then 0,(1) = 0,*(1) and 0.*(1) = 0.(1)]. This is 
impossible if statement 6 applies. 


—Received February 25, 1959 
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HE REGRES SSION analysis of mathe matic al statistics 
makes it possible to estimate the unknown parameters in 
i functional dependence on the basis of a set of experimental 
results concerning the quantities in this dependence. 

In the past, the amount of computing involved in making 
a regression analysis meant that it was suitable only for cases 
of stationary functional dependence. 

Owing to the development of high speed automatic com- 
puters, regression analysis is now being applied more widely. 
In this article we describe the method of sequential, i.e., con- 
tinuously refined, regression analysis and consider some of its 
possible applications to problems of automatic control and 
regulation. 

Consider the dependence of the random variable 
arguments 21,2, ..., 7, of the form 


Y with 


= Wa, OY [1] 


where the function f(a, v2, ... , vz) is defined and éy is a ran- 


dom quantity with expectation zero. In other words, 
X2,..., te) 1s the conditional mathematical expectation 
of Y 


..., te) = f(t, te, [2] 
We can use Eq. 2, which is usually called the regression 
equation of Y on a, 2%, ..., 2%, to approximate to the ex- 
pected value of Y corresponding to given values of x21, x2, .-., 
- It is obvious that the relative accuracy of this estimate of 
‘ will ae on he relation between the standard devia- 
Further, taking to be independent of 21, 2, ..., 
shall consider the regression equation in the form 


Ue, WE 


(Y /21,2¢, 


= f(x » V2, 
= 


a 
l 


i=1 


[3] 


where the W(x, 22, , t,) are given single valued functions, 
and the a; are constant coefficients, the re gression parameters. 

The regression analysis of the random variable Y satisfying 
conditions 1 and 3 consists in finding estimates of the un 
known coefficients a;. 

If the regression parameters a; are unknown, then the ir 
‘stimates a; can be found on the basis of some numbe or 
n(n 2 1) of actual values of Y given in the form of n hg in 
k + 1)-dimensional space: 22;,..., = 1,..., 
1). The estimate of the a; is usually found by the me of 


Translated from Jzvestiia Akademii Nauk SSSR, Otdelenie 
ekhnicheskikh Nauk, Energetikai Avtomatika (Bulletin of the Acad- 
my of Sciences USSR, Div. Tech. Sci., Power and Automation), 1960, 
10. 2. Translated by Research Information Service, New York. 
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Sequential Regression Analysis 
and Its Application to Some 
Problems of Automatic Control 


I. I. 


least squares, which consists in finding the a; which make the 
function Q 


n l 2 


j=1 i=1 


a minimum. 

It is easy to show (1)! that when éy has a normal distribu- 
tion the method of least squares gives the maximum likelihood 
estimates a;. 

Equating the partial derivatives of Q to zero, we have 


where | 


From Eq. 5 we obtain a system of / equations to determine 


» Uj) 


4 
os + a2 + 
j= 


a>, = Vii 
j=l j=1 


j=1 
j=1 j= 
j=1 j=l 


j=1 


Wuyi 


j=l 


[6] 


— The determinant of this system is Gram’s determinant, 
which is equal to the sum of the squares of the /th order minors 
of the matrix A (2) 


when < n. 

It follows from this that system 6 gives a unique solution for 
a(t = 1,...,/) if A contains at least / linearly independent 
rows. This is equivalent to the condition that there is at 


1 Numbers in parentheses indicate References at end of paper. 
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BS: 


2(m—1)k+ i(t) 


least one point among the n points of the /-dimensional space 
generated by the axes Yi, Yo, ..., Ww which does not lie on the 
same hyperplane of the (J — 1)th order as all the others. 


Consider the process }(t) which changes over time, and 
depends on the arguments 2;(¢), x(t), . , a(t) and the 
parameters po, ..., p, and is of the form 


Pr] + dy(t) 


Sequential Regression Analysis 


Y(t) = flai(t), 2o(t), , Sme(t), Po, [8] 


where 

z(t) = x(t), z(t) = Mell), z(t) = x,(t) 
Zr+1(t) = a(t — At), = — Ad), 
24(t) = At) 


x(t 1) At], Zm-1)k+2(t) = 
(m — 1)At],..., = — (m — 1) At] 


At is a given time interval; 6y(t) is a random function of time 
with zero expectation 


flat), 2.(t), Di, pr] = 
ai(pr, po. DrlWiler(t), [9] 
i=1 
where y; are given single valued functions, and a;(p;, ps, .. .? 
pr) are the regression coefficients of V(t) on 2:(t), g(t), ... > 


Zni(t), to be determined. 
We note that any function 


which can be represented in the neighborhood of the point 
Zo(210, 220, -- +; Zmio) by a finite number of terms of its Taylor 
series expansion, satisfies condition 9 in the neighborhood of 
this point. We must also note that, to the necessary degree 
of accuracy, relation 9 expresses the dependence between the; 
output variable f(t) and the input variables (disturbances) 
x(t) (¢ = 1,..., k) ina dynamic nonlinear system without 
feedback consisting of a set of linear dynamic stages W to- 
gether with static (memoryless) nonlinear stages F’, where 
the relation between the output of a given (nth) dynamic 
stage, u,(t), and the input signal of this stage, v,(é), is defined 
by the equation 

= halt — [10] 

where /,,(f — 7) is a continuous function satisfying the condi- 


tion 


— 7) | < 


a>0o A>0O O<€—7)< 


and the relation between the output of a given (mth) non- 
linear stage, u,,(t), and its input signals v,,;(4) (7 = 1, ... 
is defined by the equation 

= Vaslt), , [12] 

Here F,, is a single valued function which has a convergent 
Taylor series expansion in the neighborhood of any point 
lying in the region of possible values of the arguments. 
These statements will be proved in the Appendix. 

Let the process }(é) satisfy conditions 8 and 9 at integral 
moments of relative time 7(7 = 1,2,...,m,...,2 — 1, n) 
where 7 = t/At, and At is a given, sufficiently small time in- 
terval, and suppose that observations are made during the 
process, giving data about the n — m + 1 sets of values of 
y(t), 2(7), 22(7), 2me(j). This is with the convention 
that knowing the functions y; [2;(2), 22(2), .. ., Zmu(2) ] is equiva- 
lent to knowing the n — m + 1 sets of quantities 


y(9), (j 
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duces to the form 


ay > + ae +... + 
j=n-(N-1) j=n-(N-1) 
n n 
j=n—-(N—-1) j=n—-(N-1) 
a1 > + 
j=n-(N-1) j=n-(N-1) 
n n 
j=n—-—(N-1) j=n—-(N-1) 
j=n—-(N-1) j=n-(N-1) 
n n 


where 


Using these latter quantities, the estimates of ai(~,pe, . . ., 
Pr) = a@; can be found by solving a system of equations with / 
unknowns analogous to system 6 (with n — m + 1 > mk). 
The a; thus obtained can be used to give an estimate of the 
expected value of y(n + 1) with the given z2:(n + 1) (@¢ = 1, 

It is obvious that for fixed At the accuracy of these esti- 
mates increases as n increases. 

Before, we began from the assumption that the regression 
coefficients were constant. In real processes occurring in 
technical operations, the process parameters (pi, po, . Dr) 
are not strictly constant, but depend on changes in the sur- 
rounding conditions, on wear of the apparatus, and so on. 
Hence the coefficients a;(p;, po, ..., pr) Should be thought of 
as functions of time. However, it is assumed in our division 
of the quantities affecting )(¢) into arguments and parameters 
that the quantities a;(p, po, . .., pr) have considerably less 
relative variation with time than have y;[a(t), z(t), . ., 
Zmx(t) J. 

Next, assuming that a;(i = 1,..., 7) vary with time, we 
consider the problem of obtaining estimates of the sequential 
value of a; for? = n. The values of y(j) (7 = m,..., ) ob- 
tained at time 7? = n are in this case not of equal value for 
calculating the estimate of a;(n), and this must be taken into 
account in Eq. 4 by introducing an additional factor g = 
y(n — j) into each of the terms of 


2 


i=1 


Q= gn- jl yj - 


j=m 


The function g(n — 7), which we shall call the reliability 
function of the jth reading with respect to the moment n, 
must be defined beforehand on the basis of the known sta- 
tistical properties of the process Y(t) in such a way that the 
estimates obtained give the best fit. It is obvious that the 
reliability function can always be normalized so that g(0) 
= 1. Theng(n — Jj) < 1 for (n — 7) > 0. 

Both theoretical and practical difficulties arise in the de- 
termination of the reliability function. However, by 
taking into account the slow variation of a; over time, we can 
select a time interval 7 = N At which is large enough from the 
point of view of the variation in y; during which a; can be 
taken as constant to an accuracy which is sufficient in practice. 


With these conditions, g(n — j) canbe ofthe form 
h(n —j) = 1lforj>n—WN 
h(n —j) = Oforj <n 


The system of equations giving the estimates of a; then re- 


j=n—-—(N—-1) 


j=n-(N-1]) 


ta Yi [15] 
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When continuous tests of the process Y(é) are made, i.e., as Al 0, system 15 will be transformed in accordance with the 
relations 


lim = T (i,m =1...0) 
4 


—1) 


n 


—0 j=n ) 


The uniqueness condition for the solution of system 15 then means that the segment of the trajectory which the process describes 


n /-dimensional space W over the time interval f — T < 6 < t does not lie completely on a hyperplane of the (J — 1)th order. 
We note two particular cases of Y(é). 


1 For a linear dynamic system with one input variable x(t), Eq. 9 reduces to the form 


; 

f(p = h(j — U7] 


i=j-l 


sa result of formula 6’ of the Appendix and using dimensionless time. 
The unknown values of the weight function 4(0), A(1), ..., A(Z) can be estimated from the results of a sufficiently large 
umber V(N 2 1) of observations of the process, using the solution of system 15, here of the form 


n n n 
j=n- 


(N-1) g=un—(N—1) j=n—-(N-1) j=n—(N—-1) 


n 


(QO) x(j — + >. 
n n 
j=n—-(N-1) j=n—-—(N-1) 
n n 
j=n-(N-1) j=zn-(N-1) 


2. = —Dy(y) [18] 


jJ=n-(N-1) j=n-(N-1) 


where a(j — i) = Owhen(j 
When X(t) and ¥(¢) are stationary random functions with zero expectation, as the number of observations V increases we 


have 
l n 
lim: — — t)a(j — m) = jt m=1...) 
n 
j=n—(N-1) 


for n < N where ¢,, is the autocorrelation function of x(t) and ¢g,, is the correlation function of the random variables NX (¢) 
and Y(é). 

Using 19, system 18 can then be transformed into a well-known form which is used in the experimental determination of the 
transfer function of a system during normal operation ( 3). 

2 Fora nonlinear static system | with one input variable Eq. 9 reduces to the form 


fi =X ale 


i=0 


System 15, determining the estimates of the unknown a;, becomes 


=. aN + a >. a(j) + > a(j) = 


j=n-(N-1) jJ=n-(N-1) jJ=n-(N-1) j=n—-(N—-1) 
n n n n n 
j=n-(N-1) j=n—-(N-1) j=n-(N—-1) j=n—-(N-1) j=n—-(N-1) 
n n n n n 
Q > xj) + + ae +... +a = 


j=n—-(N-1) j=n-(N-1) j=n—-(N-1) j=n-(N-1) j=n—-(N-1) 
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Sequential regression analysis can be carried out with spe- 
cial, sufficiently fast, computing units, which receive at their 
input the data of successive experiments on the process being 
studied. These units must include memory blocks for stor- 
ing the results of the last NV experiments as well as the func- 
tional transforms for calculating the given functions y;, and 
a unit for solving the system of Eqs. 15. An analog com- 
puter is sufficient for problems involving one or two unknown 
coefficients, but in more complicated cases it is necessary to 
use digital electronic machines. 


Some Applications of Sequential Regression 
Analysis to Problems of Automatic Inspection 
and Control 


Automatic sequential regression analysis can be used to 
solve a number of problems in the field of automatic control. 
We will not spend time on the obvious applications of this 
method to the study of object and process parameters, but 
shall consider some other problems. 


(a) Application to Control Systems With a Disturbance 


Consider the control object OP (Fig. 1), at whose input act 
the disturbances ro(t), .. . , 2-(t), causing the reaction 
f(t). This reaction, together with the action dy(t), means that 
there will be a disturbance of the output coordinate of the 
object y(t). Let 2x,(t) be a corrective action, and let 2x,(t) 
(¢ = 2,...,k) be the disturbances to be measured. We are 
required to find a law for the corrective action 2x;(t) which 
stabilizes the output coordinate of the object y(t) from the ef- 
fect of the disturbances 2;(t). 

If the equation of the object is known and its parameters do 
not vary with time, this problem can only be solved when the 
corrective action .7;(f) is a function of the disturbances 2;(t) 
ee k). As we know (4), control as a function of 
disturbance can completely stabilize y(t) from the effect 
of the disturbances .x,(t), ie., can make the reaction f(t) 
identically equal to zero. Regulation in terms of disturbance 
has a special role in problems where the output coordinate of 
an object with a pure lag equal to 7 has to be freed from vari- 
able disturbances with components whose period of change 
is less than 7. The usual methods of control proportional to 
deviation have proved ineffective in solving problems of this 
type. However, an essential defect in systems of control by 
a disturbance is the absence of feedback control, which causes 
output errors when the parameters of the object or of the 
regulator change. 

Because of this, it has been suggested (5,6) that for objects 
with large position lags, a combined regulation should be 
used, consisting of a control system as a function of disturb- 
ance, whose control parameters are corrected with the help 
of special computing units which receive information from— 
both the input and the output ends of the process. In (7), 
a system operating according to this principle was developed 
for the control of a memoryless linear object with onc meas- 
urable disturbance at the input. . 
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Fig. 2 


Using the relations established previously, combined con- 
trol can be applied, in the more general case, to the stabiliza- 
tion of y(t) from the effect of disturbances 2x;(t)(¢ = 2,..., k) 
in the control object, described by Eqs. 8 and 9. The control 
system in this case consists of two computing units, BY-1 
and BY-2 (Fig. 1). BY-1, which receives the data of 
discrete or continuous readings of the quantities y(0), 2(t), 
xo(t), . . . , w(t) calculates the sequential estimates of the 
coefficients a; in Eq. 9 by solving the system of Eqs. 15. 
The results which are obtained by BY-1 are sent to BY-2, 
in which the equation f(t) = 0 is solved with respect to 2x,(¢), 
and the necessary sequential value of x,(t) is determined and 
sent to the input of the object. 


(b) Application to Extremal Control 
ry 

A similar method can be used to solve problems in ex- 
tremal control. For in many such problems the dependence 
of the optimized function y on the controlled arguments 2, 
X2,..., a can be given theoretically in the form of an equa- 
tion with unknown or slowly varying coefficients. If this 
dependence satisfies Eqs. 8 and 9, then computers BY-1 
and BY-2 can be used to find the extremal value of y. In 
BY-1, just as in the previous case, is calculated the esti- 


mate of the sequential value of the coefficients a; in the equa- 


tion of the process, and these quantities are sent to BY-2, 
in which are computed the values of 2:(i = 1,..., &) which 
will give y an extremal value, using the process equations. 
When noise é6y is present, the estimates of the coefficients 
a; which have been calculated will vary within definite limits, 
thus making the value of y deviate from itsextremum. These 
deviations are used by the system as test steps. 


In a number of problems in automatic inspection the exact 


value of the controlled quantity x can be obtained only by 


special discrete control measurements, which may take a 
long time to obtain (for example, chemical analysis for each 
test). However, it is possible to obtain a continuous estimate 
of the sequential value of x by means of approximate indirect 
methods which measure some other quantity y which is 
connected with x by the regression equation 9. 

If the coefficients a; in this equation are known, then for 
any given measurement of y the best estimate of the value 
of x which corresponds to it can be found. If the regression 
coefficients change sufficiently slowly in comparison with the 
possible intervals between exact control measurements of 2, 
then the estimates of a; can be found by using sequential 
regression analysis on the difference between the exact meas- 
urements of x2 and the measurements of y which correspond 
to them. 


Appendix 
We wish to show that in a dynamic nonlinear system with- 
out feedback, consisting of stages which satisfy conditions 10 
and 11, or condition 12, the dependence of the output quan- 
tity f(t) on the disturbances 2,(t), re(t), . . . , can be 
put in the form of Eq. 9 to the necessary accuracy. 
Consider the system of Fig. 2, where the static (memory- 


X, 


x(t) 
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less) nonlinear stages F; and /’, are described by the equa- 


tions 
u(t) F, [ai(t), r2(t), x(t) [1’] 
f(t) = F.[v(t)] (2’] 


and where F, and F2 are single valued continuous functions 
which can be expanded in a convergent power series in the 
neighborhood of any point lying within the region of possible 
changes in their arguments. 

The linear dynamic stage W is described by the equation 


= fi u(t)h(t — r)dr {3’] 


where the pulse transfer function A(t — 7) is continuous and 
satisfies condition 11. 

It is obvious that if the stage W is stable, then condition 11 
can always be satisfied by a suitable choice of the correspond- 
ing values of the quantities A and a, 

Expanding F; in the neighborhood of the point 22, ..., 
x. = 0 ina finite power series with enough terms for the re- 
quired accuracy, we have 
J 
u(t) = ayo + lxi(t), 
i=1 
where dp and a; are values of the function F; and its first and 
higher partial derivatives at the point x, 2,..., 2% = 0, and 
gilai(t), xo(t', ..., a,(t)] is the product of corresponding in- 
tegral and zero powers of 2,(t), we(t),..., z(t). 

It is easy to see that when condition 11 is satisfied and the 
functions A(t — 7) and u(r) are continuous, the integral in 
Kq. 3’ can be replaced to the necessary degree of accuracy 
by the finite sum 


[4’] 


M-1 
h(mAtyu(t — mAt) At [5’] 
m=0 


where the quantities At and ./ are chosen to give the neces- 
sary accuracy. 

Transforming to relative time f = t/ At, when? = 7 we have 
M-1 


h(m)u(j — [67] 
m=0 


v(7) = 
Putting Iq. 4’ in Eq. 6’, we obtain 


M-1 


m=0 1=1 


v(7) = 


lai(j — m), — m),..., mit [7’] 


Introducing the 7k new variables z;(j) (¢ = 1, ..., Wk) 
7) = J), 9) = X2(J), 


= — 1), 


= — (VM — 1), J) = 
te [j — (M — 1)], = —(M—1)]  [8’] 


and putting Eqs. 8’ in 7’, we have 


22(3), 


i=1 i=l 
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v(j) = bo + 


+ boi gi 2e+2(J), - 


where 
=I 


h(m)ao 


Relation 9’ can be put in the form 


v(j) = bo + Av(j) 


where Av(j) = 0 when 2;(j), 22(j),..., zus(j) = = 0. 


Further, expanding F, in a power series in the neighborhood ik. 
of the point v = bp, and taking sufficient terms of the ex- | 


pansion to give the required accuracy, we obtain 


Q 
= co+ cglAv(j)]* 
q=1 


where co, ¢g are corresponding coefficients in the expansion. — 
Putting the value of Av from relation 9’ into 11’ and simpli- — 


fying, we find 


J) ] 


L 
=e + anila(s), 
l=] 


where y; is the product of corresponding integral and zero 
powers of 2;(7), 22(]), 
coefficients and b,,:. 

It is easy to see that Eqs. 4’, 9’ and 12’ correspond to condi- 
tion 9. Moreover, the functions ¢; and in Eqs. 4’, 9’ 
and 12’ are of the same nature; i.e., they are products of in- 
tegral and zero powers of the input variables 2()). 

Thus we have seen that when a disturbance, described by Eq. 
9, is supplied to the output of the stage F or W, the reaction 
at the output of stage F or W can be expressed with an arbi- 
trary accuracy in the form of Eq. 9. Quantities y; in Eq. 9 
are finite products of zero and integral powers of the variables 
Zi, ..., a& corresponding to the given as well as to some 
previous discrete moments of time. 

It is clear therefore that the reaction at the output of a 
nonlinear dynamic system without feedback, and consisting 
of a finite number of stages F and W, which satisfy the condi- 
tions given, can be expressed arbitrarily accurately in the 
form of Eq. 9. 

—Submitted om 2, 1959 
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‘| HIS article examines several problems related to the con- 
trol of an interplanetary missile. By an interplanetary 
missile we mean an apparatus moving outside the limits of 
the atmosphere. This may be an artificial Earth satellite, 
a stage of a space rocket, ete. The control system for such 
a missile must permit it to change the motion of the center 
of mass of the missile and execute arbitrary rotations about 
the center of mass; it must also be possible to keep the axis 
of the missile in a given position. 

To move the missile about the center of mass, Tsiolkovsky 
proposed two methods: Creation of a reactive force for which 
the line of action does not pass through the missile’s center 
of mass, and the use of reactive flywheels. Control of the 
missile by means of reactive engines can be analyzed relatively 
simply. However, analysis of the control with the aid of 
reactive flywheels cannot. 

For purposes of generalization, let us suppose that the num- 
ber of stabilizing flywheels in the body of the missile is equal 
ton. Use of an arbitrary number of flywheels does not intro- 
duce any complications in writing the equations of motion. 
We will consider the positions of the flywheels and the orien- 
tations of their axes within the body of the missile as arbi- 
trary, but fixed with respect to time; let us also suppose 
that the flywheels are statically and dynamically balanced 
upon their axes; in other words, their axes of rotation are the 
main central axes of inertia. This assumption is obviously 
correct for flywheels with axial symmetry. 


| Equations of Motion About the Center of Mass 


In order to set up the equations of motion, we introduce 
two coordinate systems, taking as the general starting point 
for the calculations the point O, the center of mass of the 
material system consisting of the body of the missile and the 
reactive flywheels. In Fig. 1, ONXo¥oZo is the coordinate 
system whose axes always remain parallel to the axes of a 
certain inertial system; OX YZ is the coordinate system which 
is rigidly fixed with respect to the body of the missile. The 
problem of orientation of the OX, OY and OZ axes within 
the body of the missile will be left open for now. The in- 
stantaneous angular velocity of system OX YZ relative to 
system ON o¥oZ> will be designated by w; the projections of 
this angular velocity on the OX, OY and OZ axes will be 
designated respectively by p, q and r. 

Let positive directions along the axes of suspension of 
each of the flywheels be given; we will designate the orienta- 
tion of the flywheel axes by the cosines of the angles formed 
by the positive directions of the flywheel axes with the posi- 
tive directions of axes OX, OY and OZ. Then for the ith 
flywheel, this direction cosine will be a;z, ai, and a;.. The 
angular velocity of the ith flywheel relative to the ON YZ 
system will be designated as @;, and the projection of this 
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vector on the axis of rotation of the flywheel by w;. The 
projection of vector w; on the OX, OY and OZ axes in the 
foregoing notation will be respectively wi@iy and 

The angular momentum of the material system under con- 
sideration in its relative motion about the center of mass re- 
duced to the point O will be 


K = [r X vjdm 


where 
v = velocity of an element of mass dm with respect to the 
system 
r = radius vector of the observed element of mass dm 


Integration is carried out over the entire system 8. 

Let us represent the vector K as the sum of the angular 
momentum of the body of the missile and the angular mo- 
menta of the reactive flywheels. The angular momenta of the 
flywheels can be divided into the angular momenta arising 
from translational velocities of points on the flywheels in the 
OXo¥oZ system and the angular momenta arising from the 
intrinsic rotational motions of the flywheels relative to 
the OX YZ system 


K = lr X vjdm + fi (r X Verans]dm + 
[Ir X Vreildm 


where in the first integral, integration is carried out over the 
body of the missile C, and in the other integrals it is carried 
out over the ith flywheel M;. 

The velocity of any point of the body v (including the trans- 
lational velocity of every point of the flywheel Virans) is ex- 
pressed by the angular velocity of the ON YZ system in ac- 
cordance with the Euler equation v = [ X rj. This makes 
it possible to group the first n + 1 addends 


X vjdm + [r X Virans]dm = 


Ir X (w X r)|dm 


Consequently 


x 
K = Ir X X r)|dm + [Ir X Vrei]dm 


Now let us examine the second term of the expression ob- 
tained. 

As a result of the fact that the axes of suspension of the 
flywheels in the body of the missile coincide with their prin- 
cipal axes, the angular momenta of relative motion of the 
flywheels in the OX YZ system—if we assume that they are 
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reduced to centers which are the centers of mass of the fly- 
wheels—will be @J 1, Here @2,... , de- 
note the angular velocities of the flywheels relative to the 
ONYZ system, and J;, Jo, ..., J, are their axial moments 
of inertia. Relative to the OX YZ system, the velocities of 
the centers of mass of the flywheels are equal to zero (the 
flywheels are balanced on their axes of suspension within the 
body of the missile), which results in independence of the 
corresponding angular momenta from the centers to which 
they are reduced. Taking point O as the general center of 
reduction for all the flywheels, we obtain 


Su: Ir X Vreildm = wid; 


‘rom which the expression for K will be 


x-f, ir xX 


By representing the angular momentum of the system in 
this form, we can select the orientation of the OX, OY and 
OZ axes in the body of the missile in such a way that the ex- 
ression for the projection of the angular momentum on the 
ixis of the ON YZ system will appear in its simplest form. 
Che determination of such a system of axes is important not 
nly from the point of view of simplifying the writing of the 
‘quations of motion, but also from the point of view of im- 
proving the dynamics of the stabilized missile. Locating 
planes in which the signals from the sensitive elements of the 
stabilization system are produced parallel to the coordinate 
planes of the ON YZ system, selected in the manner described, 
and situating the axes of the reactive flywheels parallel to the 
ixes of the very same system, gives us the minimum cross 
coupling between oscillations of the missile about its three 
axes of stabilization. 

It is easy to see that first term in Eq. 1 represents the angu- 
lar momentum which a system would have if its flywheels 
were stopped, were fastened rigidly to the missile’s body, and 
rotated with it as a single rigid unit. Let us choose OX, OY 
and OZ as the principal central axes of inertia for such a 
hypothetical body; and let the moments of inertia of this 
body about these axes be 1, Band C, respectively. The first 
term in Eq. 1 would then be written with maximum sim- 
plicity as 


@ X r)|dm+ ; {1] 


i=1 


Ss Ir X (w X r)|dm = Ape, + Bae, + Cre: 
where e,, e, and e, are unit vectors along the OX, OY and 
OZ axes. Using the expression for the projection of vectors 
w; on the axes of the ON YZ system, we write the second term 
of Eq. Lin the form 


n n n n 


=1 i=1 i=1 i=1 


Grouping according to unit vectors e,, e, and e., we finally 


obtain for K 


i=1 i=1 


(cr + e. [2] 


Let us apply the angular momentum theorem to the motion 
of the system with respect to its center of mass; since the 
the angular momentum vector is reduced, we may write 

dK/dt = M [3] 
where 


dK /dt = total derivative of the vector K with respect to 


time 


Fig. 1 Mutual orientation of the missile and reference systems 
ONYZ and ONX oY 


M = sum of moments due to all external forces acting 
on the system 


The external moments acting on the system are the moment 
of external perturbation M.x acting on the missile and the con- 
trol moment of the reaction engines M, (for the case where 
reaction engines! are used as auxiliaries in conjunction with 
reactive flywheels for stabilization of the missile), and 


M = M.x + M, 


Calculating the total derivative of vector K as the sum of the 
relative and translational velocities of the end of vector K in 
the OX YZ system, we obtain 

dK /dt = K,e, + Kye, + K.e: + lw X K] 


Here, K,. AK, and kK, are derivatives with respect to time of 
the projection of vector K on the OX, OY and OZ axes 


Ap + 


K, = By + WA; ; 


i=1 


i=1 


d Equating the corresponding projection of the vectors 


dK /dt and M along the axes of the ON YZ system, we obtain 
(we will omit summation subscripts | to » from here on) 


r( By + = Me, + Mi, 


By + r(Ap + 


p(Cr Mey Mis 


C# + + p( Ba + 


a( Ap + = 

‘If the reaction engines are introduced into the problem, 
then the mass of the system and all its other inertial parameters— 
moments of inertia, principal inertial axes, etc.—ought to be 
considered as variables; however, the expenditure of mass 
necessary for the control of the missile’s position about its 
center of-mass is small, and the on ange | in mass can n be ignored in 
the calculations. 
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We then transform the system obtained as follows 


Ap + + ait) J + + + + ai) + gr(C = B) = + 


Bg + aig + ait) Ji + + aizp + aig + + rp(A —C) = 


Cr — + + ast) + + ah + aiyg + act): + A) = Mart Mn — 


Let us clarify the meaning of the third terms in the left- 
hand portions of the foregoing equations. We will examine 
the motion of the ith wheel in the coordinate system whose 
axes are connected with the flywheel. As one of these axes, 
let us select the suspension axis of the wheel. Let the posi- 
tive direction be the same as that defined earlier; the origin 
is taken at the center of mass of the wheel. Let us then 
arbitrarily locate the other two axes in the plane perpendicular 
to the first axis. Owing to the axial symmetry of the wheel, 
the moments of inertia about these two axes will be equal to 
each other. Let us designate them as Jie (the equatorial 
moment of inertia of the ith gyro). 

Since the chosen system of axes represents a system of prin- 
cipal central axes of inertia with respect to the wheel in ques- 
tion, we are able to use the dynamic equations of Euler as 
the equations of motion. 

The projection of the total angular velocity of the flywheel 
on the first coordinate axis (axis of suspension) is determined 
by the equation 


= o + + + aur 


The projection of the total angular velocity of the flywheel 
on the other two coordinate axes will be designated as g’ and 
r’. In this notation, the first of the dynamic equations of 
Euler (the equation referring to the projections on the first 
coordinate axis) may be written as 


+ Gish + diyg + ait) = 


Here M’; is the projection of the total moment applied to the 
flywheel along the positive direction of its axis. It is obvious 
that .W’; is equal in absolute value to the axial moment at 
the shaft of the wheel and has a positive sign if the axial 
moment acting on the wheel forms a right-handed screw with 
the positive direction of the axis, or a negative sign if the axial 
moment acts in the opposite direction. The axial moment at 
the shaft of the wheel depends only on the active moment of 
the mechanism driving the wheel 1/’,; and on the frictional 
moment at the shaft of the gyro M’.; 


M'; = + 


Here ’,; and /’.; are algebraic quantities calculated exactly 
in the same way as the value for M’;. Thus 


J (@: + + + ait) = + Mes 


Multiplying this equation by a;; and summing with respect 
to i from 1 to n, we obtain the expression for the third term 
in the first equation of system 4 


+ ap + aig + ait)Ji = Da:.(M ‘ac + 


It is easy to obtain analogous expressions for the corre- 
sponding terms in the other two equations as well 


+ + + ast) Ji = + 


Replacing the active moments and the frictional moments 
M’,; and M’.; applied to the shafts of corresponding wheels, 
by moments M,; and ./.; directly opposing these and applied 
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to the body of the missile, we obtain finally 


+ + M-i) — + i 


By + aig + + rp(A — C) = Mexy + 


Cr + + J + pq(B A) MM + 
Mat + Med) — plead: + 


The last two terms in the right-hand sides of Eqs. 5 repre- 
sent projections of the gyroscopic moments applied to the 
body of the missile and arising from the induced precession 
of the wheels in the missile. The terms .J/,,, > a:, M,; and 
so forth, represent projections of the control moments pro- 
duced by the output members of the stabilization system; 
these moments depend on the angular coordinates of the mis- 
sile and their derivatives, and may also be determined by a 
special control program for controlling the attitude of the 
missile or by guidance signals from Earth. 

Let us now examine the case of stabilization of a missile 
with the aid of three reactive flywheels. We will assume that 
the rotational axes of the first, second and third flywheels are 
parallel to axes OX, OY and OZ, respectively, of a connected 
system OXYZ, selected in the body of the missile in the 
manner described. The positive directions of the rotational 
axes coincide with the positive directions of the respective 
axes of the OX YZ system. Then for the direction cosines 
we have in 


0 


0 
Asy = 0 1 
3, 0 1 
Designating further 


= We = Wy W3 = 
= Me Ma = 


Ma = Mez = M., 


and setting the axial moments of inertia of the wheels equal 
to each other (J; = J2 = J3 = J), we obtain the equations 
of motion for the stabilized missile in this case a 


Mez — + 


J)g+ (A — C)rp = Mexsy + Mey + May + 
Mey — witd + 


Mexe + + Mas + 
Me — + 


+ (B — A)pq 


The equations obtained are similar to the ordinary dynamic 
equations of Euler for the body of a missile. 
only difference between the left-hand portions of these equa- 


In fact, the 
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tions and the left-hand portions of the Euler equations is the 
presence of the correcting factor J in the first terms of the 
‘quations; by introducing the same corrections in the second 
ierms of the equation in the following manner 


(C—Bgr= 


nd so forth, it would be easy to obtain complete similarity 
f the left-hand portions of Eqs. 6 with the left-hand portions 
f the Euler equations. The difference would consist only in 
‘hat, instead of factors 1, B and C, we would everywhere en- 
‘counter the differences — J), (B — J}, (C J). The 

right-hand portions of Eqs. 6 are entirely similar to the right- 
iand portions of the Euler equations written for the body of 
he missile; they contain the projections of all moments which 
an be considered as external with respect to the body. In 
ddition to moments M., and M,, such moments for the body 
if the missile would be, obviously, the active moments and 
rictional moments at the wheel shafts transmitted to the 
ody, and also the gyroscopic moments of the wheels, the 
rojections of which w,rJ, w4J and so forth, are also con- 
ained in the right-hand portions of the derived equations. 
However, the noted similarity is only superficial, and is a 
esult of the circumstance that both the Euler equations and 
iiqs. 6 represent equations of motion of a system consisting 
if one or four rigid bodies, written with respect to a reference 
system in which the expressions for the projection of the angu- 
ar momentum on the coordinate axes have maximum sim- 
plicity. One should not forget that, unlike in the case of the 
uler equations, axes OX, OY and OZ (Eqs. 6 are written in 
terms of projections on these axes) are neither the principal 
nor the central axes of inertia of the missile’s body, and that 
neither the factors 1, B and C nor the differences (A — J), 
(B — J) and (C — J) which appear in Eqs. 6 coincide with 
the moments of inertia of the missile’s body about these axes. 


2 Conelusions From the Law of Conservation 


of Angular Momentum 


In order to refer the nonrotating system ONo¥oZ to the 
system OX YZ connected with the body of the missile, we will 
introduce three angles 3, Y and ¢ (Fig. 2). 

We will designate angle 3 as the angle between OX’— 
the projection of the axis OX on the plane XoOYo—and axis 
OX». Angle W is designated as the angle between axis OY 
and the line of intersection between planes XoO Yo and YOZ. 
In the given instance, the following succession of rotations 
of the body of the missile holds during transition from a posi- 
tion where the systems OX YZ and OXo¥oZo coincide, to the 
given position of the system OX YZ. First, a rotation about 
axis OZ) through the angle 3, then a rotation about the new 
position of axis OY» (see OY’ in Fig. 2) through the angle 
y, and then finally, a rotation through angle ¢ about axis 
OX. The table shows the direction cosines for the axes 
connected with the missile and the nonrotating coordinate 
system. 

Let us express the projections of the angular velocities of 
the body of the missile in terms of angles 3, Wand ¢ 


p =¢-—dsiny 
q = ¢ + sin ¢ cos 
r= 


Jcos gcos — 


Fig. 2 Sequence of rotations through angles J, y and y for 


transferring system O\,),Z, into YZ 


The motion of the missile with respect to its center of mass 
is expressed by Eq. 3. If no external moments are applied 
to the missile (M = 0), then from Eq. 3 we obtain an obvious 
consequence, which is known in mechanics as the law of con- 
servation of angular momentum 


K=C 


where C isa constant vector. 

Based on this law, we can draw a series of interesting conclu- 
sions, for the case of the problem of rotation of the missile 
about its center of mass. 

Let the projections of vector K on the axis of a coordinate 
system which is nonrotating in inertial space, be Cro, Cy, 
Cx». Then its projections in the OX YZ system of coordinates 
are connected with Co, Cyo, Co by the equations 
Cw = K,cos Ycos + K,(sin ¢ sin cos — 

cos ¢ sin 3) + K.Acos ¢ sin ¥ cos & + sin ¢ sin #) 
Cy = K,cos sind + K,(sin ¢ sin sin + 
cos g cos 3) + K.(cos sin sin — sin ¢ cos 


Co = —K,sin y+ K,sin gecs + K. cos ¢ cos 


Moreover, the following obvious equation is true — 


Quantities K,, A, and K, are easily obtained from Eqs. 2. 
Let us examine the motion of the missile characterized by 
the fact that its control system, with the aid of flywheels, 
keeps a certain plane (connected with the missile) continuously 
parallel to a given plane in the inertial space. Without loss 
of generality, we can assume that axes ON and OY of the 
moving coordinate system lie in a given plane parallel to the 


Axes of the 
nonrotating — 
system xX } 
Xo cos y-cos 6 


Yo cos y-sin 


—cos g:sin + sin g-sin y-cos 3 
cos ¢:cos 3 + sin g:sin y-sin 
sin ¢-cos 


Z 
cos ¢g:sin y-cos J + sin g-sin 
cos g:sin y-sin — sin ¢-cos 3 
cos ¢-cos 
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plane OXo¥o. Then, in unperturbed motion, Y = ¢ = 0, 
and Eqs. 8 take the form 


Co = K,cos 8 — K,sind 


Cyw = K,sind + K, cos 
3 Cx = 


In Eqs. 9, angle 3 may be a given function of time. For 
example, let the missile fly around some heavenly body along 
a plane orbit, and let it be required of the control system that 
the missile axis lying in the plane XOY remain pointed con- 
stantly to the center of the heavenly body while plane YOY 
itself be parallel to the plane of the orbit. Then, in one 
revolution, angle 3 will change by 27, since the body of the 
missile will complete one revolution about axis OZ. 

The first two equations of 9 show that under such condi- 
tions of motion about a heavenly body the quantities A, and 
K,, as well as angle 3 will be periodic functions of time. 


It is easy to see that a 
k, = C,, sin [8(t) + 6] 
Kk, = C,, cos [&(t) + 6] 
[10] 
Cry = + Cy’ 


6 = arctan (Cy0/Cy) 


In the unperturbed motion assumed here, angles Y = ¢ = 
0, and the angular velocities of the body of the missile, as can 
be seen from Eq. 7, will be 


p=0 q=0 


In this case, as follows from Eq. 2, projections of the angular 
momentum A, and A, which differ from zero can be con- 
nected only.with the rotation of flywheels which have attained 
certain angular velocities over a time period preceding the 
moment in question. In the given instance, this permits us 
to consider A, and A, projections of the angular momenta 
of the flywheels. Consequently, the rotational speed of the 


wheels whose axes are motionless with respect to the body 


of the missile, must, in general, vary periodically. In the 
simplest case, where the missile has three flywheels whose 
rotational axes are parallel to axes OX, OY and OZ, those 
wheels whose axes are parallel to OX and OY will have a 
periodically varying rotational speed. Using the notation 
adopted previously for the case where the projection of the 
angular momentum of the body of the missile on the OX and 
OY axis is equal to zero, we obtain 


Wad kK, Wyd y [11] 
where 


J.J, = moments of inertia of wheels set parallel to axes 
OX and OY, respectively 
w,,w, = their angular velocities 


(subscripts 7 and w are omitted, since there is no possibility 
of confusion in this case). 

From Eqs. 10 and 11 it follows that w, and w, are periodic 
functions of time. Over the time period of a single revolu- 
tion of the missile about the central body, the angular ve- 
locity of each wheel attains a maximum value, drops to zero, 
changes its sign, and attaining in absolute magnitude the 
very same maximum value, again falls to zero and then re- 
turns to its initial value. This behavior of the wheels is 
understandable if we consider that for r + 0 the angular 
velocities of wheels w, and w, produce gyroscopic moments 
which differ from zero in the first two equations of 6; these 


moments can be eliminated by suitable acceleration or de- — 
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celeration of the wheels. It appears, indeed, that damping 
of the gyroscopic moments is necessary for fulfillment of the 
requirements that p = O and q = 0. 

It is important to note that these accelerations and de- 
celerations, which are unavoidably connected with an ex- 
penditure of energy, are not at all necessary for orderly motion 
in orbit. In fact, if at the initial moment the projections of 
vector K on the plane of the orbit were equal to zero, i.e., 
Co = Cy = 0, then, as can be seen from Eqs. 10 and 11, #, = 
w, = 0 during the entire time of motion in the given orbit. 
Thus, if one is aiming for a reduction in the expenditure of 
energy, then the reactive wheels whose axes lie in the plane 
of the orbit should be slowed down (for example, by em- 
ploying moments created by reaction motors). 

Such braking of the wheels is not always possible; for this 
reason let us examine more closely the problem of expendi- 
ture of energy on the periodic speeding up of wheels. 

The kinetic energy of the wheels is expressed by the equa- 
tions 


9 J,w,? = ( ry Si [d(t) 
. + 9] 


Comparing these two equations term by term, it is easy to 
see that where J, # J, the total kinetic energy of the wheels 
is a periodic function of time. Consequently, in the general 
case, part of the mechanical energy of the wheels must 
periodically be transformed into another form of energy (for 
example, electrical energy), which must then once again 
undergo transformation into rotational energy. It is to be 
noted, that even when J, = J,, the braking of one wheel and 
the simultaneous speedup of another is practically always 
convenient to achieve by means of the transformation of the 
kinetic energy of the braked wheel into another form of 
energy and the speedup of the second wheel by some motor 
which is set into motion as a result of the application of non- 
mechanical energy. ‘Thus, in the general case, the operation 
of wheels is inescapably related to the transformation of 
energy, and this, in turn, means additional losses. We will 
now evaluate these losses. 

Let the efficiency of transformation of mechanical energy 
into some other form of energy which can then be retrans- 
formed into mechanical energy be m, and the efficiency of 
retransformation be ym». The quantity of mechanical energy 
transformed into the new form of energy when the wheel is 
braked completely will be equal to (1/2J)mC.,*, and the ir- 
reversible losses of mechanical energy will be (1/2J) X 
(1 — m)C.,?. Upon acceleration of the wheel, an expendi- 
ture of energy amounting to (1/2/72)C,,? takes place, which 
exceeds the achieved kinetic energy by (1/n2 — 1)(1/2/)C;,?. 
Thus, the total loss of energy over the time period from brak- 
ing to acceleration of the wheel will be 


{ 


The equation obtained shows that for given C,,, m and 
n2, it is possible to obtain a reduction in losses only by in- 
creasing the moment of inertia of the wheel. 

An increase in J is however connected with an increase in 
the dimensions or the mass of the wheel or a simultaneous 
increase in both. Let us therefore examine the problem of 
the optimum values of J for two wheels, with the condition 
that the total mass of these wheels will be limited. The 
total loss of energy of both wheels over a single cycle of de- 
celeration and acceleration will be 
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In view of the given construction of the wheel, its mass m can 
usually be approximated by an exponential function of its 
moment of inertia, m = aJ”, where a is a constant and n > 0. 
The condition that n > 0 is entirely natural, since with an 
increase in J the mass of the wheel also increases. 

The condition of constancy of the overall mass of the two 
wheels is expressed by the equation gs 


Jy" =b [12] 
where 6 is a constant. 7 
Let us find the minimum AF satisfying the conditions of 
Kq. 12. Using the Lagrange method, we will seek an ex- 
tremum of a function Z of two variables J, and J, 


Z 2 rales n) NOE diy — b) 


where \ is the Lagrange multiplier. After elementary com- 
putations, we obtain the condition for an extremum J, = Jy. 
It is easy to see that it is the condition for a minimum. 

Thus, to reduce losses of energy in a control system, it is 
necessary, as far as possible, to brake the wheels, and if this 
is not possible, then to increase their moments of inertia. 
In so doing, it is necessary to select equal moments of inertia 
for wheels which interact during motion in orbit. 


3 Example of Application of the Derived 
Equations to the Case of Small Deflections of the 
Missile From a Given Attitude 


As an example, let us examine the case of stabilization of a 
missile about a given, unvarying position in space. Let the 
stabilization system be required to maintain the missile’s 
position, and during this let coordinate system ON YZ always 
coincide with a certain nonrotating coordinate system 

The mutual positions of systems ON YZ and OXo) Zo will 
be characterized as before by the angles 3, Y, ¢. As for the 
orientation of the wheel axes in the body of the missile and 
the magnitudes of their axial moments of inertia, we will 
make the same assumptions we made previously in deriving 
Eqs. 6. We will neglect the magnitude of the axial moment 
of inertia of the wheels J relative to the moments of inertia 
A, Band C. We will assume further that M.. = M, = 0 
and will neglect the magnitudes of the gyroscopic moments and 
frictional moments relative to the magnitudes of the control 
moments V/4., Way, Ja: We will limit ourselves to the case 
of small deviations of the missile from the given position. 
eqs. 7 then are simplified 
q = r= J 
if we suppose further that the derivatives of the projections 
of the angular velocity /, ¢, * are of the same order of smallness 
as the projections p, g, r themselves, then the second terms of 
the left-hand sides of Eqs. 6 can be eliminated, since they are 
small terms of higher orders, and the equations take the form 
Ag = M,, By=M, 


p=¢ 


the angle in the control program. 


Let us suppose that the moment at the shaft of the reactive 
wheel mounted in the direction of axis OX is a function of 
the angular velocity it develops, and, in addition, depends 
(through some control system) on angle ¢g and its derivative ¢ 


If we can write analogous expressions for .1/7,, and W/,:, we 


obtain 
Ag = 
Ci = 


From this it follows that the oscillations of the missile 
about axes OX, OY and OZ will become independent. 

Let us consider the first of these equations in greater de- 
tail. We expand the function 1/’,, into a Taylor series near 
the point (wo, 0, 0), corresponding to an unperturbed pro- 
gram. Retaining the first members of the expansion, we 


obtain 


where the partial derivatives are evaluated at point (@,o, 0, 0). 
The motion of the missile in plane YOZ will be described by 
the equation 


e+ e+ bw, = Vas(wo, 0, 0) 
In order to eliminate 6w,, we use the first equation of 8. 
For the case of the small angles 3. neglecting terms of 
the first and higher orders of smallness—we obtain C.o = 
K, or Co = AG + Jw,. Variation in the neighborhood of 
the point of interest vields 


Consequently, the final equation for motion of the missile 
in plane YOZ will have the form 


OW”. A OM’... OM’... 
Ag+ 


The term .1/,,(@,0,0,0) on the right-hand side of the equa- 
tion represents the forcing function determining the static 
error of the system 


Fstatic = 


It is interesting to note that the dependence of the moment 
M,, on the angular velocity of the wheel for 0.1/'..,0¢ <0 
facilitates damping of the disturbance. For a sufficiently 
strong dependence of the moment J/,, on the angular ve- 
locity w,, it is not necessary to introduce the derivative of 
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account of the attitude control problem for space vehicles can 


Reviewer’s Comment 


The problem of controlling the orientation of a rocket in 
space is related to a more general problem in mechanics, i.e., 
the dynamics of rigid bodies. The usual problem in me- 
chanics prescribes the external forces and constraints and 
requires a solution for the subsequent motion. In the pres- 
ent case, the objective is to specify the desired motion and 
then to determine the control forces and moments, so that 
the desired result is obtained. A very thorough survey of 
the past work on the dynamics of rigid bodies can be found 
in the book by Lemanis and Minorsky (1). In the case 
of a spinning body, the analysis of the control by reaction 
engines is not relatively simple. An analysis of this prob- 
lem is contained in a recent note by Suddath (2). A detailed 


be found in a paper by Roberson (3). This article contains 
an extensive bibliography. An account of a NASA study of 
attitude control for space vehicles is also available (4). 


-BERNARD PAIEWONSKY 


> Aeronautical Research Associates of Princeton 
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\ OLECULAR flow conditions are considered for bodies 
4 of any shape (not necessarily convex). The problem 
is found to be that of solving an integral equation, which in 
form is the same as one in the theory of illumination; the 
problem may be solved by means of photometric measure- 
ments. A solution is given for the interior of a spherical sur- 
face, and a use of the solution is discussed. 


1. Molecular flow occurs when Knudsen’s number A = 
\// is large (A is the mean free path and / is a characteristic 
length); A > 2 gives effective molecular flow. An isolated 
simple body having a convex surface at every point presents 
no problem if the law of molecular reflection is known; 
Patterson (1)! gives a detailed treatment, with references to 
the literature. It is usual to neglect the effects of reflected 
molecules on the incident flow (the incident particles are 
assumed to have an undisturbed Maxwellian velocity dis- 
tribution). The problem becomes much more difficult. if 
there are several bodies separated by distances comparable 
to their own dimensions or if the shape is locally concave. 
Here allowance must be made for repeated reflections, which 
give rise to interactions. The parameters of the incident 
flux are among the quantities to be calculated. 

The laws of molecular reflection are important. Experi- 

Translated from J/zvestiia Akademii Nauk SSSR, Otdelenie 
Tekhnicheskikh Nauk, Mekhanika i Mashinostroenie (Bulletin 
of the Academy of Sciences USSR, Div. Tech. Sci., Mechanics 
and Machine Construction), 1960, no. 3, pp. 117-120. Translated 
by Research Information Service, New York. 

1 Numbers in parentheses indicate References at end of paper. 
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ment [see (1), p. 168] shows that 95-98% of the incident 
molecules are reflected diffusely from ordinary polished sur- 
faces. I shall neglect the few particles that show a correla- 
tion between the angles of incidence and reflection. The re- 
flected molecules have a Boltzmann energy distribution whose 
temperature may be close to (but perhaps not equal to) the 
temperature of the surface. 


2. ‘The number of particles striking an element dop per 
unit time surrounding a point P at a surface with an interior 
normal vp is the sum of the number incident directly and of 
the number incident after reflection. The first component 
depends only on the direction of vp and is 


N\(Pidop = + Vr B(1 + erf@)|dop [2.1] 


where 
n = density (in particles per cm?) 
T = temperature of the flow ~*~ 


R = gas constant 


The second component is to be deduced by adopting the fol- 
lowing device. The particles are reflected diffusely, so the 
flux originating at any surface is equivalent to a distribution 
of sources having a density g (particles per em’). Lambert’s 
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law is obeyed; the reflected molecules behave as though they 
come from some hypothetical gas on the far side of the sur- 
face. The number of impacts on the surface from particles 
whose directions of motion lie within a solid angle dQ whose 
axis is inclined at an angle ¢ to the normal is proportional to 
cos gdQ. Thus [q cos gdQda/z] of the particles reflected by 
the element do enter that solid angle (Fig. 1). This second 
component is now to be found as 

N(P)dop 


log [2.2] 


q(Q)argare 
dop da 
S(P) 


TP pg” 


in which 


vp-(rp — Ye) 


apg > 
Irp 
Teg = Irp 
The integral is taken over the part S(P) of the surface 


that is visible from P (Fig. 1). 

The main difficulty lies in deducing the distribution of 
sources; the equation is easily derived from the condition 
of conservation for a small volume 


Tr pg” 


= Ni(P) + NAP) = Ni(P) + [2.3] 


S(P) 
which is a linear integral equation whose solution is essen- 
tially the solution to the whole problem. 

The temperature of the body and the aerodynamic forces 
can be calculated if the distribution is known. I shall assume 
that the accommodation coefficient for the temperature is 
unity (normal values range from 0.9 to 1), so a reflected mole- 
cule has the temperature of the point from which it was re- 
flected. Patterson (1) gives the contributions of the re- 
flected molecules to the energy and momentum fluxes (pp. 
159-178). By analogy with Eq. 2.2, the energy flux incident 
on do» is found as 


e(T 9)q(Q) apeaer 
dor — [2.4] 
S(P) TT 
where 
Tg = surface temperature at point Q q 
e(T) = mean molecular energy at temperature T , 4 


Kach reflected molecule carries, on the average, energy e(7')p, 
so the net energy reaching the body via dap is 


on e(Tr) lq(Q) apeder 


(2.5] 
TT pa? 


dog 


e(P)dop = dor 


S(P) 


[2.6] 


(rp — re) 
dop dag 
S(P) 


TT pa? 


in which 7(7') is the mean absolute momentum of a molecule 
at temperature 7. The reflected molecules make a contribu- 
tion 


1i(T p)N2(P)vp [2.7] 


because they are reflected diffusely. The additional momen- 


tum flux is then 


q(Qap 
S(P) 


TT pa? 


io(P)dop = 
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The incident flux of Eq. 2.1 brings in a momentum 5b , 


which produces at P an extra pressure p: and an extra tangen-_ 
tial force t2, which are respectively 


pAP) = i(P)-vp 2(P) = in(P) — po(P)vp [2.9] 


3. Eq. 2.3 has the same form as the equation describing the 
illumination in a space having non-absorbing walls which re- 
flect in accordance with Lambert's law. Here the source den- 
sity is the intensity of illumination produced by unit surface el- 
ement, and the incident flux is that produced by the external 
sources. The analogy can be used to produce optical analogs 
in order to solve Eq. 2.3; the scale of the analog is without sig- 
nificance. Measurements are made of the intensity produced 
by external sources corresponding to N,, which may be multi- 
plied by any factor. Green’s function for Eq. 2.3 can be repro- 
duced in the usual way. The solution for large Mach num- 
bers (Newtonian flow) is particularly easy to find, because 
the model need only be placed in a parallel beam of light. 

It is difficult to solve Eq. 2.3 in general; some exact solu- 
tions are given in textbooks on light, but they are of little 
value here. The interaction between two bodies is governed 
by the ratio of the areas to the square of the separation; if 
that ratio is small, successive approximation methods give 
the result quickly. 


4. Consider as an example the repeated reflection within 
a concave spherical surface of radius R; here part of the 
integrand takes the very simple form 


4R? 
The integrals may be simplified further if the Mach number 
is large; then 


1 


and has the solution © 


P) nos’ (4 4] 
( 4k? S 
- 

where 

S = area of the part of the spherical surface exposed to 

the flow 
S’ = projection of that area on the flow direction. (It is 


supposed that the flow strikes the entire surface.) 


We can use Eq. 4.4 to consider a spherical molecular flow 
diffuser, a device that may be of value for air intakes for use 
at high altitudes. The diffuser has the shape shown in Fig. 2. 


1631 


= 
= 
- 
’ 
7 4 
Fig. 1 7 
| — i 
and the equation takes the form ss 
Ne 


. 
Reviewer Ss Com me 


nt 


m4 
The author presents a theoretical analysis of free molecular 
flows with incident molecules coming from both the free 
stream and from diffuse reflection from a surface. Perhaps 
the most interesting part of the paper is the suggested use of 
the mathematical analogy between aerodynamic flow and 


Wehave 
— 2xR%cos 


g(P) = nv cosg + 


cos? 6 
nv (cos 5) [4.5] 


and from the small hole there emerge 


cos? — cos 
2 — cos 6 


rR?nv — cos? 6+ [4.6] 


molecules; then the emergent flux is increased by one half 
of the incident flux. 
—Received November 20, 1959 


Reference 
1 Patterson, G. N., ‘‘Molecular Flow of Gases,” 1956, 
= * 


~~ 
= 
illumination. This analogy permits solutions to aerodyna- 


mic problems to be obtained from measurements of light 
intensity. The 50% increase in exit momentum mentioned 
in Section 4 is obtained as @ approaches zero. 


—F rank W. Barry 
Hamilton Standard Div. 
United Aircraft Corp. 


Luminosity of Matter Confined Within 
a Finite Volume and Having Arbitrary 
Absorption and Emission Bands savsoy 


ECONDARY processes of absorption and emission which 

take place in volumes of finite dimensions influence the 
intensity, yield, shape of the lines, the duration, and other 
characteristics of the radiation emitted from the volume. It 
is very important to take this effect into account in astro- — 
physics (1),! in the investigation of gas discharges, or in the © 
studies of the luminescent properties of substances with 


Translated from Akademiia Nauk BSSR, Doklady, 1959, vol. 
3, no. 12, pp. 479-483. Translated by A. Werner, New York. 
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strongly overlapping absorption and emission bands, ete. 
The following equation (2,3) underlies the investigation of the 
effect of secondary processes 


| 
fy, vOK + (em 7) U1 
l 4 
K(p) = = inp? €(w) y(w)b(w)e~ [2] 
+a 
T) = ¥(w)e(w)dw [3] 
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Here, y is the ratio of density of emitted particles to the 
density of all particles; 7 is the ratio of time ¢ to the duration 
of luminescence of the elementary volume 79; w is the fre- 
quency, expressed in terms of the half width of the emission 
band computed from its maximum [w = (vy — v)/Av]; y(w) 
denotes quantum yield of elementary volume; €(w) and 6(w) 
the shapes of emission and absorption bands. Also 


e(w)dw = | 


Volume V’ and radius vector p are incorporated in Eq. 1 as 
dimensionless magnitudes: p, derivative of radius vector r 
over absorption coefficient at the maximum of the emission 
band ko; V’, derivative of volume V over ko* [total absorption 
coefficient equals kob(w) |. Eq. 1 is written without taking the 
forced emission into account and is valid when the emission 
band is independent of the frequency of stimulating light. 

The free term yoF(p, 7) characterizes the possible particle 
densities at point p, stimulated without the participation of 
secondary absorption and emission processes. The analytical 
form of this term depends on the condition of stimulation. 

Iq. 1 for a constant y(w) was formulated by Biberman 
(2). Agranovich (3) assumes the value of y(w) dependent on 
the frequency of stimulating light; however, applicability of 
Equation [1] is subject to the following restrictions: 

1 The law of radiation damping of the elementary volume 
is identically exponential] for all emission band frequencies. 

2 The dependence of the quantum yield on the stimulat- 
ing light frequency is determined by inactive absorption (ex- 
tinction processes of the first kind), and consequently the 
duration of luminosity is not dependent on the spectrum com- 
position of light. 

Many solutions of a particular kind for Iq. 1 for the coinci- 
dence of frequencies of the emitted and absorbed quanta 
(resonant emission) are given in literature (1,4). Solution 1, 
in a general case, is extremely complicated owing to the 
existence of a frequency integral. An attempt is made in the 
present work to reduce Eq. 1 to a better known equation for 
resonant radiation. 

As a rule, the total density of radiating particles y(p) is of 
the same order of magnitude with the density, which is de- 
termined (4) by the magnitude of the free term yoF(p, 7). 
Consequently, the error of solution caused by the approximate 
kernel, is of the same order of magnitude as the error of 
kernel.? Therefore, if it is possible to substitute the kernel of 
the integral Eq. 1 with a simpler kernel, incorporated in the 
integral equation for resonant radiation, the specified problem 
will be satisfied. 

Comparison of Eq. 1 with the corresponding equation of the 
propagation of resonant radiation (4) shows that they differ 
only in the kernel, which in the latter case is of the form 


b(0) = 1 [4] 


Kres(p) = 


[5] 


Comparison of Eq. 2 with Eq. 5 shows that the kernels 
differ by the value of f(p). In the resonant case (fres) has an 
exponential form. In the case of Eq. 2 it has the form of a sum 
(integral) of exponentials. Kernels K and AKy,s rapidly de- 
crease with increasing argument p and, consequently, solution 
of Eq. 1 is dependent essentially on the behavior of kernels 
only for small values of the argument [for sufficiently uniform 
y(p)]. This circumstance indicates that substitution of f(p) 
with a simpler function fre; must be such that they do not 
differ too much for small values of p. 

Expanding f(p) and fres(p) into series in steps of p, and ad- 


2 It can be shown that for a stationary system, if 
vo p — poly(p)dV' 


then Ay «< y (AK is the absolute error of kernel, Ay the densities 
of radiating particles). 
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justing the first terms of the expansion to each other, we have 
for f(p) the approximate expression which incorporates only 
one constant Ker; 

S(p) =~ filo) = (6] 


where 


[7] 


ker = f €(w)b(w)dw 


Magnitude kerr is the correction factor for the absorption co- 
efficient. For e(w) = 6(w) or for b(w) = 1, kerr is equal to 1. 

It is clear that Eq. 6 gives values close to the exact ones only 
for small values of the argument. It is possible to find the 
formula for f(p) which gives a better approximation when two 
constants are introduced. Comparing the two first terms of 
the expansion f(p) and f,..(p), we have for f(p) an expression 
which is determined with the aid of two constants 


f(p) ~ folp) = ak’ [3] 
where 
ak’ ets kets = J €(w)b2(w)dw [9] 


Here k’.«, as before, is the correction factor for the absorption 
coefficient, a is the correction factor for the yield yo. As for 
kerr, the magnitude of k’er, and @ equal 1 when the emission 
band e(w) is much narrower than the absorption band 6(w). 

The foregoing expressions for the multiplicative factor f(p) 
of the kernel] of the integral equation give a good approxima- 
tion to the exact values for small p. In order to evaluate the 
nature of the approximation for p, several concrete calcula- 
tions for the following y(w), e(w) and b(w) were made 


e(w) = b(w) = [10] 
Yo VT 
¥(w) 1 
=> ) 11 


Magnitude A in Eq. 10 characterizes (for constant y) in half 
widths the distance between maxima of the emission and 
absorption bands. The values 0 and 0.5 have been selected 
for the calculation of A. The calculated data are compiled in 
Table 1. 

Exact values of f(p) were calculated with the aid of quadra- 
ture formulas of the Gaussian type (5). If the absorption 
and emission bands have a dispersional form (Iq. 11) a 
simple equation is derived for f(p) 


where Jo) and J; are Bessel functions of the imaginary argu- 
ment. The values f(p) computed in accordance with this 
formula are also given in Table 1. This table shows that for 
p< 1,fiand fo only differ slightly from the exact values. The 
values fo(p), as expected, approximate f(p) better than fi(p), 
even for p > 1, whereas fo(p) < f(p). For example, fo(p) is 
smaller than f(p), on the average, for p = 1.0 by 1.5-4%; 
for p = 1.5 by 5-10%; for p = 2.0 by 10-20%. For the same 
p, the difference between f,(p) and f(p) is greater, whereas for 
a rather considerable range of variation of p (to p = 3-4), fi(p) 
is greater than f(p). 

As already pointed out, the kernel of the integral equation 
decreases considerably more rapidly than f(p), and, conse- 
quently, the error in the value of the kernel for large values 
of the argument must not be strongly reflected in the solution 
of the equation. 

In Fig. 1 is shown the behavior of the kernel of the integral 
equation for a plane parallel layer, when the shape of the 
absorption and emission bands are given by Eq. 10. The 


[12] 
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Table 1 
——A=0— A=0.5 ——  ———-Scattered from———— 

p f(e) file) file) Sop) 

0.0 0.707 0.707 0.707 0.801 0.801 0.801 0.500 0.500 0.500 

0.2 0.601 0.614 0.601 0.656 0.683 0.655 0.431 0.452 0.430 

0.4 0.512 0.533 0.510 0.540 0.582 0.536 0.372 0.409 0.370 

0.6 0.437 0.462 0.433 0.445 0.495 0.438 0.323 0.370 0.319 

0.8 0.373 0.402 0.368 0.369 0.422 0.358 0.280 0.335 0.274 

1.0 0.319 0.349 0.313 0.307 0.360 0.293 0.244 0 303 0 236 

1.5 0.219 0.245 0.208 0.198 0.241 0.177 0.176 0.236 0.162 

2.0 0.152 0.172 0.138 0.132 0.161 0.107 0.129 0.184 0.112 

3.0 0.077 0.085 0.061 0.064 0.072 0.039 0.074 0.112 0.053 

4.0 0.042 0.042 0.027 0.036 0.032 0.014 0.047 0.068 0.025 

Kr . with a kernel of the equation for resonant radiation by in- 
F | troducing two constants, dependent on the shape of the emis- 
06 | sion and absorption bands. In this case, inclusion of two 
\ \ constants (in accordance with formula 9) must give satisfac- 
\ A-0 A=05 tory results: 
wk y 1 For comparatively small volumes (the optical thickness 
\ at the maximum of the emission band is of the order of magni- 
a2 ‘ k tude of two to three units). 

‘Y 2 Inthe case when stimulation of particles is concentrated 

i Pil on a comparatively small portion of the volume. 
0 10. 20 30. ~OP 0 1.6 20 «3.0 p If the stimulated particles are distributed more or less uni- 


solid curve corresponds to to the kernel 


€(w)b(w) [b(w) p [13] 


The crosses and circles correspond to approximate values of 


kernels according to equations ao 
x Ki(p) = 
O Ke(p) = (ak'ets/2) Ei(k' essp) 


where E;(p) is the integral exponential function. — 

The figure shows that the approximate values (crosses and 
circles) agree well with the exact curve. The agreement 
is better, the smaller the argument. Values K:2(p), for 
example, in the interval of p from 0 to 2, differ from the exact 
values by one or two units in the second decimal place. The 
values of K(p) vary from @ to 0.04 in the same interval. 

The problem, then, of determining the density of radiating 
particles of a substance with arbitrary absorption and emis- 
sion bands confined in a volume of finite dimensions can be re- 
duced to an identical problem for resonant radiation (coinci- 
dence of frequencies of emitted and absorbed quanta). Only 
the kernel of the corresponding integral equation is replaced 


formly over a large volume, better results are probably ob- 
tained by introducing only one constant kes; determined by 
formula 11. This can be expected on the basis of the follow- 
ing representations. For small p, fi < f, Ai < A; and for 
large p, conversely, fi > f, Ki > AK. Thus the error AK is 
such that the integral of AK over an infinitely large volume 
equals zero. For a uniform distribution of the density of ra- 
diating particles (y = constant), the error Ay, dependent on 
such a kernel, is minimal. (See footnote 2.) This implies 
that in certain instances it is expedient to include a’ contained 
within the following limits instead of magnitude a 


kerr << a’ <1 [14] 


In conclusion, I should like to express my appreciation to 
8. I. Stepanov, Member of the Academy of Sciences BSSR, for 
his invaluable advice. 
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Reviewer’s Comment 


In this paper, it is shown that the problem of the secondary 
absorption and emission processes in a medium with arbitrary 
spectral bands confined in a finite volume may be reduced to 
an identical problem for resonance radiation. On replacing the 
kernel of the corresponding transport equation by the kernel 
of the diffusion equation of resonance radiation and introduc- 
ing two constants which depend on the shape of the spectral 
bands, the transport integral equation appropriate to the 
numerical analysis is derived. 

In what follows, it is of interest to mention some articles 
concerning the resonance radiation not cited here. Assum- 
ing that the escape of resonance quanta can be treated 
statistically as a decay of atoms in the excited state in a man- 
ner similar to that used in the diffusion process, Milne (6) 
first formulated and solved approximately the equation of 
transfer for the imprisoned resonance radiation. 

Holstein (3.4), who presented a rigorous method for 
handling the imprisonment by means of the integral equation, 
investigated the rate of emission of the gases by determining 
the cessation of the optical excitation for a variety of spectral 
shapes in plane and cylindrical geometries. With the aid of 
the method of discrete ordinates, Chandrasekhar (1) and 
Edmonds (2) have obtained the solutions of the diffusion 
equation of imprisoned radiation in the first and the second 
approximations, respectively. Based on the integral method 

of Holstein, Walsh (8) has computed approximately the den- 


sity and the imprisonment lifetime of resonance radiation in 
a gaseous discharge between parallel plates. The particular 
case with Doppler broadening was investigated. 

With the aid of the diffusion treatment of resonance radia- 
tion based on the two-fluid theory, Little (5) showed that a 
linear dependence of the decay time on pressure is in good 
agreement with experiment in natural mercury, using the 
measured values of lifetimes and cross sections and assuming 
Doppler broadening. Recently, Veklenko (7) also treated 
the diffusion problem of the resonance radiation. Applying 
the Ambarzumian transform, i.e., the point-to-plane trans- 
form, to the transient transport equation for the resonance 
radiation, a one-dimensional integro-differential equation is 
obtained. With the aid of the Fourier transform, the source 
function is given. Then, using the inverse Fourier and Am- 
barzumian transforms provides a complete solution of the 
problem for an arbitrary shape of spectral lines. 


UENO 


Rand Corp. 


1 Chandrasekhar, S., Radiative Transfer, Clarendon Press, Oxford, 1950. 

2 Edmonds, F. N., Jr., Phys. Rev., 1950, vol. 78, p. 424. 

3 Holstein, T., Phys. Rev., 1947, vol. 72, p. 1212. 

4 Holstein, T., Phys. Rev., 1951, vol. 83, p. 1159. 

5 Little, P. F., in Proc. Fourth Int. Conf. on Ionization Phenomena in 
Gases, N. R., Nilsson, Ed., North Holland Pub. Co., 1959, p. 202-205. 

6 Milne, E. A., Jr., J. Math. Soc. London, 1926, vol. 1, p. 40. 

7 Veklenko, V. A., Zh. Eksper. Teoret. Fiz. (Soviet Phys.—JETP), 1959, 
vol. 36, p. 138. 

8 Walsh, P. J., Phys. Rev., 1957, vol. 107, p. 338. 


A STUDY is made of the existence and stability of a 


steady-state plasma torus; the velocity vector v and 


the magnetic vector H are assumed to be collinear 
H = = xV4rpv 


in which % is a scalar constant and p is density. Shafranoy 
(1)! has detected such formations, but he did not examine the 
solution in detail; he supposed that such a ring can be stable 
only if the velocity of the gas in the ring is larger than the 
corresponding Alfven speed (v 2 HV i.e., if It 
is found here that such a ring can be stable for « > 1 subject 
to certain restrictions on the perturbations. Further, the 
pressure within the ring can exceed the external pressure very 
greatly, so the temperature can be high. Such a ring may be 
a model for ball lightning; some of the effects observed with 
such lightning (motion in closed spaces, magnetization of 
metal objects, might be explainable in terms of such a model 
(2). 


Translated from Jzvestiia Akademii Nauk SSSR, Otdelenie 
Tekhnicheskikh Nauk, Mekhanika i Mashinostroenie (Bulletin of 
the Academy of Sciences USSR, Div. Tech. Sci., Mechanics and 
Machine Construction), 1960, no. 4, pp. 7-13. Translated by Re- 
search Information Service, New York. 

‘ Numbers in parentheses indicate References at end of paper. 
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Magneto-Vortical Rings 


YU. P. LADIKOY 


The precise significance of the model is determined by the 


lifetimes of such rings and by the possibility of observing 


them. 


1 Steady-State Motion of a Magneto-Vortical 
Ring 


1 We assume that the interior of the ring is filled with an 
incompressible liquid having infinite conductivity. The 
azimuthally directed vortexes and currents cause the flow and 
the magnetic field in meridional planes. The equations of 
magnetohydrodynamics then become 


v? 
(1 — «2 (rot v X v) = —grad (? + 5) 


which is integrable for the axially symmetric case (3) if Q/r = 
f(y). This condition on the current function WY gives us an 
equation in polar coordinates r, 0, z 
tloy doy 


{1 2] 
or? Or + Oz ri) 


of i 
div v = 0 rot. vy = 2 


and for the surface of the torus 


v 


constant 


which are related to r, 0, z by 


asin & 
cosh + cos 


a sinh 


~ cosh + cos = 


[1.4] 


The coordinate surfaces are then: For a torus 


7 = constant (r — acoth n)? + 2? = (a/sin n)? 
for a sphere 


for a plane 


constant a cot &)? + r? = (a/sin &)? 


(z 


constant 
The square of a line element takes the form 


ds* = [a?/(cosh + cos £)*] X 
(dé? + sinh? dé? + dn?) 


[1.5] 


The radius ¢ of the cross section is small relative to a, the 
radius of the ring (Fig. 1); then the y for points near the 
surface on both sides are very large, so x = exp(—7) is very 
small. 

Hicks (4) has solved Eq. 1.2 subject to Eq. 1.3 and c/a « 1 
on the assumption that f(y) = 1 = constant; the approxi- 
mate form is 


+B+Qx?- Ox? cos E+... | [1.6] 


—Ta/4m) 


1636.0 


[1.3] 


In what follows toroid coordinates § 6, 7 shall be used — 


Now H = 


so 


2a*b? 


H, = 


(2x :) 


xV 4rp x? sin & 


The integral of Eq. 1.1 is 


p 


y+ (1.10] 


WE, n) 


c — (v?/2 


and is related to T by 
= velocity circulation i in the rnal flow vds = ff Qde = 
B = constant 
The cross section is not assumed to be an exact circle; the 4rra*h? ath? 
equation of the bounding surface is assumed to be 
Substituting into Eq. 1.10, we have 
The condition Y = constant gives us that 8; = 9b/4; the eu? 9 
other 6, are smaller quantities. The projections of the (x x? cos :) 
velocity are related to the current function by ; as 
where 
ron : u = Q/a*bis the principal value of the velocity at the sur- 
face of the ring 
pe 2) of\? po = pressure at the central line for x = 0 | 
or oz 
The total pressure at the surface is 
p+ = mt (1 (1+ beos [1.12] 
= a’? 4 
=> 
1" s5g3* g, 2 We assign to the whole mass of fluid a velocity V, and 
consider the external flow about a torus having at infinity a 
Q | velocity V and a given [ at the surface of the ring. A mag- 
‘oie Jah? (2x — % x* cos ) ea netic field (caused by the current in the ring) exists in the ex- 
sed ternal space, but this field is not related to the motion of the 


medium. 
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Hick’s solution for this V and T takes the following approximate form for points near the ring 


y= — 4x cos &) + (1 xeo 8 Jal (in : +1)]- 
x X 


2 2 ‘OS és 
In b 4) a 4 xX (in h 3) {1.1 3 | 


sO 
J 10 (in 1 el The first formulas in Eqs. 1.18 and 1.19 show that large pres- 
: ay | 20) x 2 sures (and temperatures) can occur with the ring if the mag- 
' 4 netic energy exceeds the kinetic energy. The ring will move 
= Eg + 2Q sin in a sense opposed to the sense of the circulation at points 
— <a° x lying on the axis of symmetry (an incident flow will compress 


( 4 ' the ring). In the converse case, the ring can be a region of 
= J, Ea — 49 (in 5) | reduced pressure, in which case it will move in the other sense. 
ax \ 2q) x - The ring will not move (V will be zero) if 


4 lp\'? 


The external m: ignetic field satisfies the same equations and » 
boundary conditions, so, putting V = 0 and = (in 4 
which J is current) in Eq. 1.14, we have a p (in _ ) p yi 
pi b 4 
/ \ 
x V 4 l 
= 1 + E (in >) | COS Then the pressure at the center is 
7 7 Po = Po + (pi — p)(u?/2) 
i, =~ iin - = 2) will exceed the pressure at infinity if > p. 
2a x : 
x 2 


parameters not accompanied by a change in the ring shape 
For x = b{1 + (96/4) cos &] the total external pressure on 
the ring is ee 


H? 4 l 


The normal components of the field and velocity are zero at 


the surface of the ring, so we have the condition ‘ 
eae 7 will produce only a change in the speed of the forward motion; 
ip + H?/8r; = 0 [1.17] it will not produce instability. Therefore we shall examine 
(Here and in what follows denotes the difference in the the surince, whore equation in the 
values of A across the surface of discontinuity.) The co- turbed state is 
efficients to cos n&é (n = 0.1) in Eqs. 1.12 and 1.16 give us ‘ 
( = b(1 + 2b cos &) [2.1] 
Po = po — (pr +s + iil ‘Ue = P= — \s Basset (4), who examined the out-of-plane oscillations 
of a vortex ring, we shall assume that any perturbation is of 
a ; the second order of smallness relative to 6 = ¢/2a, i.e., that 
aye. 2 b in the equation for the perturbed state is 
20 4 lp l : 
V = [1 — x? In 1 + x? — 22 
K + x + + nO + at) [2.2] 
: x2 Sa 1+ x? 1 7 in which @ is of the order of b. This is equivalent to the as- 
In sumption that the unperturbed cross section is the 
where xX = exp (—y) = b n = contant 3| 
u = speed of the circulation at the surface Then the unperturbed velocity, field and pressure are: In the 
c= 2ab = mean radius of the cross section internal region 
v = p,/p, ratio of the densities 
for = we have ox = ve° = = H;° = 0 


Po = Po — pur(1 — x?) @ 
& 
V = ub (in (1 — 


u? 
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in the external region Te have for the internal region that 
—i(A + m)v, + 203 — = —(imb/x) I(x) 
—(A + — 2xnh,) = —2bnIl(x) 
+ — 24, + + + imhs) = OII’(x 
(A + m)hiy = moj 3) 
+ 2xinvs — v'sx 0 
The perturbed values are : imh, + 2xinhs — h’sx — hs = 0 
v=v°+v BH =H° +h p = po + p& 
This system gives Bessel’s equation for v. 
in which the perturbations v, h and p are high orders of 
small quantities relative to v°, H° and p°; these perturba- » dv2 dv 2 
tions must (after linearization) satisfy these equations: For J dy? +Y 
the internal region 
in which 
+ grad (v°-v) — v° X rotv — v X rot v° + y= 47? 
r= [A+ m(1 — + m)? — 12m?) 
The solution bounded at x = 0 gives us that 
Oh/dt = rot [v°-h] + rot [v-H°] 
vo = Adn(y) [2.15] 
(2 6) ~ whence we have from Eq. 2.8 that 7 
and for the external region A = [2bnB(A + m)(L — [yo m(yo) + 
ov 1 yo = 2bnV4r? — 1 
All other unknowns are expressed as functions of v2 
rot v = 0 divh = 0 roth = 0 
2Qnx(1 — 4r2) 
The junction between the two solutions at the surface of the si 


° 
perturbed torus is governed by the conditions Yd m(y) + 2mrJd,.(y) 


oF HY A+m 2nx(1 — 4r*) 
+v-gradF =0, {H,} = 0, 
ot Saf IT(x) = 
2bn(X + m)(1 ) 
F(é, 0, x, t) = 6[1 + Bexp (imE + + tot)] — [29] 
The total pressure from within at the surface is 
= 
Pts =), + Sa +a + oP (imé + + iot) = 
2 — m*x?\(1 — 47? 2) 
1 — x?) x 6 t) [2.18 
The projection of the velocity in the external region is to be 
i found by solving Neumann’s problem (the boundary condi- 
We put tions of Eq. 2.8 are used). The solutions are 
= [vg(x)ig + v0(x)ie + Vq(x)i,] exp + nO + ot) mB(A + m) Kn(2xn) Kyn(2xn) 
y= — = +m) 
the(xdig + he + (x)iehy(x)in] exp i(mE + nO + ot) 
+ m)K’,,(2xn) Bm? Kn(2xn) 
pt = pr(x) exp + nb + ot) 2xnk",,(2bn) 
lransferring to dimensionless quantities by means of he = K’,,(2bn) K’..(2bn) 


a to the accuracy implied by the assumptions; here K Kn(2 xn) 
are MaeDonald functions of mth order. From Eq. 2.7 we 
Pp: = pu*ll(x) o = —(u/2ab)r have 
[2.11] I(x) = B(A + [2.20] 
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The total pressure from outside at the surface is 


H? H° J re) H°he\ 
pt+ (» + + Eee 4 =) exp t (m&, + 


u 


9 u? 9 9 9 9 
Po t+ pr? — px? + + m)*? — px?m3] 


Then from Eqs. 2.18 and 2.21 we have for \ the equation 
(A + m)? — x*m?](1 — 47?) J m(Yo) 
Yo! 'm(Yo) + m(Yo) 


2bnK'm(2bn) +1—y=0 [2.22] 


The motion is stable if the roots do not have positive imagi- 
nary parts. Let us assume that 2bn < 1, i.e., that the per- 
turbations are not of very high frequency. 

First let us examine the case m = 0 (pinch disturbance); 


(4 — A*)So(yo)/ Yo ’o(yo) + vA? In bn + — 1 = O 


[2.23] 


Let us suppose that A] = 1; then 


and Kq. 2.23 gives us 


= 2(v — 1)/(1 — 2b2n2v In bn) [2.24] 


The motion is stable (A is real) if vy = pi/p > 1; if A is very 
small, yo ~ 4bn A, and Eq. 2.23 becomes 


4J (Yo) Yod ‘a(Yo) = 0 [2.25] 


which has only real roots if y 2 1, or two conjugate purely 
imaginary roots if y < 1 (5). 


Kn(2xn) et(mé+nd+ ot) 
K’,,(2bn) 2bnf 


[2.21] 


Now we shall turn to the cases m 2 1, of which only 


yo = 2bnV4r? — 1] K An «1 


is considered. Here Eq. 2.22 becomes 


+ v) + 2A[m(1 + v) — 1] + 
mim — +» — 2x?) = [2.26] 


Stability is indicated by the condition 


2x?m(m — 1)\(1 + vy) + — 1) + 1>0 [2.27] 
Iq. 2.26 is simply the equation for the out-of-plane vibra- 
tions of a vortex ring (4) if x = 0. Thus the ring can be 
stable in every case if the density within the ring does not ex- 
ceed the density outside it. 

I am indebted to L. I. Sedov for guidance in this work. 


—Received February 20, 1960 


References 


1 Shafranov, V. D., ‘‘Magnetic Vortex Rings,’ Zh. Eksper. Teoret. Fiz., 
1957, vol. 33, no. 3 (9). 

2 Stekol’nikov, I. 8., Physics of Lightning and Lightning Protection, 
Moscow-Leningrad, 1943. 

3. Lamb, H., Hydrodynamics, Gostekhizdat, 1947, p. 165. 

4 Basset, A., A Treatise on Hydrodynamics, Vol. 2, Cambridge, 1888, 
chap. XIV. 

5 Lebedev, N. N., Special Functions and Their Uses, GITTL, 1953. 


Reviewer’s Comment 


This interesting paper by Ladikov appears to be one of 
several Russian efforts to understand a relatively simple 
axially symmetric plasma torous configuration, moving 
steadily in free space. This configuration lends mathematical 
simplification but retains the possibility of some interesting 
physical behavior. The ring current is azimuthal, and it is 
assumed that the ring plasma velocity is parallel to the mag- 
netic field. The entire ring moves along the axis of symmetry. 
Thus, the plasma behaves somewhat like a spinning smoke ring 
in free space. A single fluid, continuum MHD approach is 
used throughout. 

I know of no similar work in Western journals. The prob- 
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lem is reminiscent of plasma toroids experimentally produced 
by Bostick (Phys. Rev., 1956, vol. 104, no. 2, pp. 292-299), al- 
though in these experiments the toroid motion was along the 
plane of the ring. 

Although there is a wealth of analysis and experiments on 
confined plasma toroids resulting from controlled thermo- 
nuclear efforts, the free space plasma toroid seems to have 
been overlooked. Ladikov mentions some lightning phe- 
nomena which behave in a manner similar to the model 
analyzed. 

—Rosert A. GRoss 
Department of Mechanical Engineering 
Columbia University 
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AGNETIC measurements on the second cosmic rocket 
I had the following aims: To discover if a magnetic field 
of dipole type) exists on the moon, and if it exists, to find its 
magnitude; to obtain new data on the magnetic field in the 
outer corpuscular zone (radiation belt) of Earth. Previous 
information on the latter was obtained by the measurements 
on the first cosmic rocket (1).! 

Solutions to both these problems are of fundamental im- 
portance for explaining the magnetism of the Earth and cos- 
mic bodies, and for explaining magnetic disturbances and 
their relation to solar activity, polar aurorae and cosmic rays. 

The rocket container sent to the moon on Sept. 12, 1959 
contained three independent single-component magnetome- 
ters with saturable-core detectors of even harmonics. The 
detectors measured the components of the magnetic field along 
each of three mutually perpendicular axes Y, Y, Z associated 
with the rocket container. 

For measurements of weak fields, the only kind we might 
expect close to the moon, the sensitivity of the magnetometers 
in the second rocket was increased by approximately four 
times over that of the magnetometers in the first rocket. 
This involved a reduction in the measurement range of the 
magnetometers, and this was the reason why measurement of 
Karth’s magnetic field by the second rocket did not begin 
until it was approximately 18,000 km from Earth’s center; 
on the first rocket, measurement of the field began at a dis- 
tance of approximately 14,700 km. 

The magnetometers used were relative instruments. The 
values of magnetometer zeros were determined first of all from 
the results of measurements at great distances (outside Earth’s 
field). The degree of stability of the zeros of magnetometers 
X, Y and Z is illustrated by the graphs in Fig. 1. Each of the 
points represents the mean reading during 20 min of flight. 

As the graphs show, the zeros of channels Z and Y remained 
constant to within +20 y close to Earth. The zero of channel 
X varied a little. For the part of the trajectory near Earth, it 
may be said that the zero shifted at a constant rate of around 
7 y per hr. On the part of the trajectory outside Earth’s mag- 
netic field, the mean error in determining the zeros of the 
magnetometer channels was +30 y. 

In addition, it was possible to check the zeros of the 
magnetometer channels in parts of the trajectory within 
EKarth’s magnetic field from the measured values of the field 
components and the total vector 7 caleulated on the basis of 
these measurements. 

Let us consider the rotation of a detector set at a constant 
angle y to the horizon and rotating round a vertical axis with 
angular velocity ¢. The detector is at a point where the 
vector Ty is inclined at angle @ to the horizon. The detector 
readings will vary periodically - 

Hact.. = To(sin a sin + cos @ cos sin 

Translated from Artificial Earth Satellites, 1960, no. 5, pp. 

16-23. Translated by Research Information Service, New 


York. 
‘ Numbers in parentheses indicate References at end of paper. 
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The mean value is given by the first term of this expression. 
The readings of another detector lying in the same plane as 
the first, but at an angle of 90 deg to it, will be 


Hact.2 = To(sin a@ sin Y + cos @ cos yp Cos ¢) 


ic., they will also be periodic functions of ¢ with the same 
mean value, but for each position there will be a phase shift of 
between Hact, and Haet. ». 

The difference in levels relative to which the readings of two 
such detectors vary (i.e., the variation in the mean values of 
Hace, and Hacer, 2) indicates the zero shift in one of the chan- 
nels. In the coordinate system associated with vector 7) the 
readings on the detectors of a three-component magnetometer 
will be 


7 


Z = T) cos 0 N = Ty sin sin ¢ 
Y = 7) sin cos ¢ 
where 
6 = angle between the axis of detector Z and vector 7% 


¢ = angle of inclination of the axis of detector VY to the 
plane of the meridian 


It is obvious that the square of the modulus T) (7)? = XY? + 
2 + Z?*) does not depend on the orientation. However, if 
the zeros in channels Z, VY, ¥ shift through values a, 6, c, re- 
spectively, then the expression for the square of the total vee- 
tor will contain two additional terms 


T? = (a2 + +c?) + 
27,(a cos @ + b sin @ sin g + ¢ sin 4 cos ¢) 


aw 


- 30 
40 = 
ot 
50 
x 0 
100 150 2U0 


Fig. 1 Readings of magnetometers 7, ) and \ at 75-200 x 10° 
km trom Earth, showing stability of zeros 
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The second term in this expression reduces to the total shift 
of the experimental curve and does not depend on the orienta- 
tion. 

The third term depends on the magnitude of the field and 
on the orientation of the container. It may be found on the 
experimental curve of the total vector of the field 7 from the 
modulation of the amplitude with a period equal to the period 
of rotation of the container. From the results it is possible to 
determine the magnitude and sign of the zero shifts of the de- 
tectors which are affected by the container rotation. 

Thus, the container rotation provides an additional oppor- 
tunity to exclude the effects of the zero shifts, and this, to a 
large extent, overcomes the limitations imposed by the rela- 
tive nature of the apparatus used 

The errors due to nonsimultaneity of interrogation of the 
magnetometer channels attain magnitudes of about 30 y only 
in the initial region of the measurements. The error in de- 
termining the total vector was approximately 50 y. 


Vleasurements of Earth’s Magnetic Field 


The results of the field measurements are shown in Fig. 2. 
The values of 7 calculated from the measured values of VY, ¥, 
Z are shown here by points ‘The continuous curve 2 gives the 
values of Earth’s total magnetic field intensity, calculated 
from the results of a spherical harmonic analysis of magnetic 
field measuremen‘s at Earth’s surface. 

From analysis of the results of measurements in the region 
near Earth and the differences between the measured and 
calculated values of the total field intensity, we can draw the 
following conclusions. 

The measured and caleulated values of the total com- 
ponent of the field agree fairly well with one another at dis- 
tances of more than 23,000 km from Earth’s center. At 
shorter distances, we find discrepancies similar to those ob- 
served in measurements on the first rocket, but of much 
smal er magnitude. 

The intensity of corpuscular radiation in the outer radiation 
zone was measured by American scientists with the Pioneer 
II] and TV rockets (2,3!, and by Soviet scientists with the 
first and second cosmic rockets (4-6) 
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Magnitude of total vector of Earth’s magnetic field as function of distance from center. 
calculated dependence for dipole field; 3: difference AJ’ = Tog. — Tea. between meas- 
ured and calculated intensity, 4: 


1: Meas- 


smoothed curve for difference \7’ 
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Fig. 3 Radiation intensity (Fig. 3a) and difference (Fig. 3b) 
AT = Treas — Teac in region of outer radiation belt in relation to 
distance from Earth’s center. 1,2: Radiation measurements 
on rocket Pioneer III (ascending and descending legs); 3: meas- 
urements on rocket Pioneer IV; 4,5: measurements of total 
ionization on first and second cosmic rockets (from scintillation 
counter data); 6,7; differences \7’ from magnetometer meas- 
urements on first and second cosmic rockets 
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Table 1 
Decrease in 
=| Distance from mean diurnal 
 Earth’s center value of hori- 
to zone maximum, zontal compo- 
Rocket 103 km nent AH,7 
Ist Cosmic (USSR) 25 
Pioneer III (U. 8.) 22-24 70 
2nd Cosmic (USSR) 18 78 7 


When the rocket Pioneer III returned to Earth it had made 
two sets of measurements in the radiation belt: One on its 
way out from the Earth and one on its return. The results of 
the measurements (2-5) are shown in Fig. 3a. Curves 1, 2 
and 3 show the intensity of corpuscular radiation in relation 
to the distance from Earth’s center, as recorded in (2 and 3) 
(the radiation intensity is characterized by the number of 
particles incident on 1 cm? in 1 sterad per 1 sec). Curves 4 
and 5 are plotted from the data of (4 and 5) (the radiation in- 
tensity is characterized by the total energy evolution in the 
erystal of the scintillation counter). The ascending and 
descending legs of the trajectory of Pioneer III intersected the 
corpuscular zone at different distances from the geomagnetic 
equator, which explains the difference in the position of the 
radiation maximum recorded on ascent and descent of 
Pioneer III. 

The magnetic field in the outer radiation zone has been 
on the first and second cosmic rockets. 
3b show the differences between the 
values of the total intensity of Earth’s magnetic field, 
measured on the first and second cosmic rockets, and the 
theoretically calculated values from the results of spherical 
harmonic analysis. 

An examination of the available data shows that the di- 
mensions of the outer zone, the intensity of radiation in it, 
and the distance of the region of maximum radiation intensity 
from Earth varies with time (from measurement to measure- 
ment)(2,5,6).. These changes are correlated with magnetic 
storms and polar aurorae, indicating that the outer corpuscu- 
lar region is formed by solar corpuscles trapped by the mag- 
netic field. The analysis reveals a definite relationship be- 
tween the position of radiation maxima and magnetic activity. 

Table 1 shows the distances, taken from Fig. 3a, of the 
corpuscular radiation maxima from Earth’s center and the 
diminution in the mean diurnal values of the horizontal com- 
ponent:(AH). As Table 1 reveals, the greatest radiation in- 
tensity occurs closer to Earth, the more violent the magnetic 
storm preceding the formation of the corpuscular region. 


measured twice 
Curves 6 and 7 on Fig. 
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variation of the mean diurnal values of the horizontal 
Comparing the graphs in 
Figs. 3 and 4, Dolginov suggested that the radiation maxi- 
mum is probably located farther from Earth, the greater the 


The 
component is shown in Fig. 4. 


number of calm days preceding the experiment. The follow- 
ing model of the effect seems possible, Under the action of the 
corpuscular stream, the zone is pushed toward Earth. When 
the main phases of the magnetic storm are completed, the 
zone is repelled by the inhomogeneous magnetic field of Earth 
into the region of smaller fields and gradients (like the 
ionized gases of a burner in the inhomogeneous field of an 
electromagnet). Thus, we can postulate that the outer zone 
of corpuscular radiation possesses diamagnetic properties. 
This last hypothesis does not contradict the fact of the non- 
coincidence of the cosmic radiation maximum and the maxi- 
mum depression of the magnetic field, a feature discovered in 
the flight of the first cosmic rocket. 

A comparison of the intensity variations of the corpuscular 
radiation and magnetic field with distance, as recorded during 
the flight of the first cosmic rocket, shows that the region 
where magnetic effects associated with the outer radiation 
zone were observed was located several thousands of kilome- 
ters closer to Earth than the maximum intensity of corpuscu- 
lar radiation. Ifa similar relationship had occurred during the 
flight of the second cosmic rocket and, what is more, if this 
time the region of magnetic effects had been located closer to 
the maximum of the corpuscular radiation intensity, such 
magnetic effects could still not have been recorded, since they 
lay in that region where the field intensity was high and was 
outside the measurement range of the magnetometer. 

The pattern noted suggests that the outer corpuscular re- 
gion acts as a special intermediate reservoir of solar corpuscles, 
which fills up during magnetic disturbances and gradually 
loses them in the period between disturbances. This, in fact, 
explains why aurorae are observed at high latitudes every 
day, and not only in days of magnetie disturbances. 

In the first communication on the experimental discovery 
of the magnetic field of the outer corpuscular region (1), it 
was suggested that one of the most probable causes of the 
magnetism of the radiation zone might be electric drift- 
currents, created by the drift of charged particles in Earth’s 
magnetic field. The possibility of a drift of charged particles 
in Earth’s magnetic field was confirmed by observations of 
high altitude atomic explosions conducted in the U. S. (7). 

As in modern theories of the origin of magnetic storms and 
aurorae, we can assume that the outer corpuscular zone is 
formed by the entry of solar corpuscular streams, consisting 
of neutral particles and different numbers of positive and 
negatively charged particles (protons, positive ions and 
electrons), into Earth’s magnetic field. When such streams 
flow around Earth, some of the charged particles will be 
trapped by the magnetic field and become part of the outer 
corpuscular zone of Earth. 
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Fig. 4 Time variation of mean diurnal values of H from ground station measure- 


ments during flights of cosmic rockets. 
P-4, Pioneer IV; 2-CR, second cosmic rocket 


Launching times: P-3, Pioneer III; 
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n Charged particles in the radiation zone complete numerous 
passages from one hemisphere to the other, spiraling around 
the field lines, and at the same time drifting in a direction per- 
pendicular to the field and its gradient. Positive particles 
move toward the west, and negative particles toward the east. 
It can be shown (8) that, owing to the small radius of curva- 
ture of the electron tracks, their contribution to the magnetic 
field of drift-currents in the corpuscular zone will be insig- 
nificant. Hence, the magnetic effects recorded by mag- 
netometers on rockets must be attributed to drift-currents of 
protons and positive ions. 

The possibility of the existence of protons in the outer 
corpuscular region is indicated also by observations of 
aurorae, the spectrum of which reveals hydrogen lines. At 
present there are no other experimental data on the presence 
of protons with energies of several tens of kiloelectron-volts, 
but indirect evidence of the presence of such protons in the 
corpuscular zone can be obtained from the results of magnetic 
measurements in this region and from the recording of hydro- 
gen lines in auroral spectra. 

From the data of the magnetic measurements on the first 
cosmic rocket Pushkov put forward the hypothesis that the 
two regions of magnetic field reduction, recorded by the 
magnetometers on this rocket (one at a distance of 23,000 km, 
and the other at a distance of 25,000 km, from Earth’s center), 
were associated with two maxima of proton intensities, one of 
which was formed during the great magnetic storm of Dee. 
3-6, 1958, and the other during the two moderate magnetic 
storms observed during the period Dee. 12-20, 1958. 

An analysis of the variations of the field components 
showed that during flight, the container precessed in such a 
way that the magnetometers Y and Y rotated round the Z 
axis, and the Z axis itself rotated round another axis, the 
axis of precession, the spatial orientation of which remained 
constant. The angle of nutation (the angle between the axis 
of precession and the Z axis) was about 89 deg. The con- 
tainer’s natural rotation period, manifested in the periodical 
changes of the XY and Y components, was approximately 900 
sec, and the period of precession, manifested in the periodic 
variation of the Z component, was approximately 86 sec. 

The values of the XY, ¥ and Z components varied periodi- 
cally in Earth’s magnetic field, owing to the rotation of the 
container and the magnetometer detectors rigidly attached 
to it. This rotation continued after the container had left 
Earth’s magnetic field, and it could provide an additional eri- 
terion for a decision as to the existence of a field close to the 
moon and in interplanetary space. If a field exists in these 
regions, we should observe the same periodic variations of the 
field components as we observe close to Earth. 


Measurements of the Moon’s Magnetic Field 


The magnetometers on the second cosmic rocket were also 
intended for detecting the moon’s magnetic field. They 
functioned normally right up until the container struck the 
moon. The last measurement was made approximately 50 km 
from the moon’s surface. The measurements did not reveal 
a magnetic field near the moon. 

From an analysis of the precision of measurement of Earth’s 


Table 2 
Distance Distance Distance 
from from from 
moon's moon's moon's 
center, center, center, 
km T,Y km km Ty 
5150 30 4000 70 2800 50 
5100 80 3800 70 2700 110 
4950 50 3700 30 2600 40 
$850 50 3400 90 2290 40 
4700 20 3300 40 2160 170 
4550 50 3200 70 2055 60 
4400 70 3100 70 1935 50 
4300 40 3000 60 1795 120 
Note: The values are rounded off to tens of gammas. 
The most reliable values, gaged from the quality of the 
telemetric traces, are given in bold face. 


magnetic field and the accuracy of the telemetric records of 
the measurements, we can conclude that if the moon had a 
magnetic field of intensity more than 50-100 y at the moon’s 
surface, this field would have been detected. 

Table 2 gives the measured values of the total magnetic 
field intensities at different distances from the moon's center. 

As Table 2 reveals, the measured values of the field lie 
within the limits of measurement errors, and do not show any 
noticeable increase of field with approach to the moon. We 
know that if the moon possessed a general field like Earth’s, 
similar to that of a homogeneously magnetized sphere, then, 
with approach to the moon, there would have been an in- 
crease in field inversely proportional to the cube of the dis- 
tance from the moon's center. In view of the data given, we 
can conclude that either the moon possesses no magnetic 
field at all, or its field, even at the surface, is so small that 
it lies wholly within the limits of measurement error. The 
greatest possible value of the mean intensity of magnetization 
of the moon is about 0.0002 cgs units. In any case, We can 
state that the intensity of the lunar magnetic field at the 
moon's surface is no more than 1/400 of the magnetic field 
intensity at Earth's surface (the mean value of the moon's 
magnetism is not more than 0.257% of the mean value of 
Earth’s magnetism). 
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Reviewers’ Comments 


This paper is highly significant. The Russian scientists 
are to be congratulated for a pioneering effort and for making 
a prompt report on their work available to scientists all over 
the world. The measurements themselves deal with one of 
the fundamental physical parameters of space science and at 
least set limits on the magnitude of magnetic fields in space 
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and near the moon, estimates of which have been varied and 
highly controversial. 

Although the instrumentation used by the Russians in the 
measurement of weak magnetic fields was crude by present 
technology, the instrumentation operated as part of a complex 
device which achieved a spectacular success. The Russians 
very honestly reported on the limitations of their instrumen- 
tation, and the fact that the experiment was simple probably 
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insured success. Probably the most significant criticism in 
regard to the instrumentation is that the range of zero in- 
stability in their experiment is approximately equal to the 
range of actual field variations observed by Explorer X, which 
contained a sophisticated magnetometer, in traversing the 
same range of distance from Earth. 

The paper reports on both measurements in space and 
measurements in the vicinity of the moon. These will be 
commented on separately. 

A significant result of the Russian experiment was that it 
measured a low field value near 23,000 km. These data are 
controversial. The Russians attribute this phenomenon to an 
actual field disturbance due to particles in the outer belt. 
However, such a measurement can also result from an in- 
accurate representation of Earth’s field in their spherical 
harmonic analysis or trajectory errors in which the spacecraft 
was actually at a slightly greater altitude or lower latitude 
than computed. This combination could easily lead to large 
discrepancies. If one accepts the Russian interpretation, it 
implies a highly accurate and sophisticated tracking system. 

On the other hand, it is possible that the Russian result is 
entirely correct. At 23,000 km, the spacecraft was appar- 
ently located near the inner edge of the outer Van Allen belt, 
or the outer edge of the slot, depending upon its latitude in 
geomagnetic coordinates. There is recent supporting evidence 
from a NASA Javelin sounding rocket, which carried a 
rubidium vapor magnetometer and was fired at Wallops 
Island, and there was a diamagnetic type of field disturbance 
of around 200 y in the horn of the outer belt near its inner 
edge at the time of this flight. Vangu: tard ITI also showed dis- 


turbances, mainly next to the slot. This satellite, however, 
did not reach the outer belt. Finally, preliminary interpreta- 
tions of Explorer X data show that the field is depressed in 
intensity over the entire range of the outer part of the inner 
belt, the slot, and the inner part of the outer belt in the region 
traversed over South Africa. If this is not found to be the 
result of an error in trajectory computation, there will be ad- 
ditional reason to believe that the radiation belts are mag- 
netically active in the vicinity of the slot. The result of these 
indications is that the Russian data may be valid. Future 
experiments should resolve this problem, and NASA mag- 
netometers are scheduled for flight on S-3 spacecraft in the 
near future that should cover the region of interest. 

Concerning the moon’s magnetic field, we accept the Rus- 
sian measurements as showing the lunar field to be less than 
50 y. However, this is not too surprising, since the moon’s 
field may be quite weak and influenced greatly by the mag- 
netic field of the sun. Perturbations of the lunar field by the 
solar ‘‘wind’’ would be expected. In fact, the magnetic field of 
the moon, if present at all, may be blown out into a teardrop 
shape by solar influence. 

In making their measurement of the lunar magnetic field, 
the Russians have pioneered and have set an upper limit which 
will be of considerable value to space scientists in understand- 
ing the nature of the lunar body. 


—James P. HEPPNER and 
JoHN W. TOWNSEND JR. 

NASA Goddard 
Space Flight Center 
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Propagation and Interaction of 
Compression and Rarefaction Waves 
in Elastic-Plastic Media 


CchAnunnindins (1,2)! was the first to examine the propagation of compression and unloading waves 


LIAKHOV and 


G. M 


in elastic-plastic media. He obtained an analytical solution for loads which instantaneously in- 
crease and then decrease, according to the power law, when the relationship ¢ = f(€) was linear. 
The solution to a series of problems connected with the propagation and interaction of waves was 
provided by Shapiro (3), Lebedev (4), and Biderman (5) employing a method of plotting grid charac- 
Lee (6), Kaliski and Osiecki (7), and others have examined the propagation and interaction 
of waves where there is an infinitely large unloading velocity. B. A. Olisov, C. Lempson, N. V. 
Zvolinsk, S. S. Grigoryan and others haye examined the application of elastic-plastic medium 


This article will examine the propagation, reflection and refraction of fixed plane compression 
and rarefaction waves in elastic-plastic media, employing the Lagrangian mass-time coordinate 
system, when the law for the compressibility and unloading of media is linear. This made it possi- 
ble, on the basis of (8), to obtain a formula defining the propagation and interaction of waves. 


Translated from /zvestiia Akademii Nauk SSSR, Otdel. Tekh. Nauk, Mekhanika i Mashinostroenie (Bull. 
Acad. Sci. USSR, Div. Tech. Sci., Mechanics and Machine Construction), 1960, no. 3, pp. 99-108. Translated 
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1 Compression and Unloading of Media 


KT US examine a medium in which compression and un- 

loading to some stress o, occurs elastically, i.e., in ac- 
cordance with a single law; when o > o,, it occurs according 
to various laws (Fig. 1). When the load in the medium is 
removed (if ¢, > o,), residual deformation ¢€ is preserved; its 
magnitude depends on the maximum stress ¢, which is 
achieved during compression. A secondary load takes place 
according to the unloading law to the stress reached with the 
first compression. 

The o = f(e) relationship can be made linear, both during 
loading and unloading of the medium, but with different 
straight lines. 

In the future, the system of magnitudes p, p, used in gas- 
dynamies, will be used. In addition 

— =. - 
Vo p 


where p is the medium’s density. 

If the o = f(e) relationship were linear, then the connection 
of pressure p and specific volume V would also be linear. Let 
us assume that when loading 


Pp — po = —Ag(V — Vo) when p<ps_ [1.1] 
p— p= — when [1.2] 
p— pe = —AX(V — Vi) when p> [1.3] 
where 
Ao = Cop 
‘Ay = 
Co, (1 = propagation speeds of elastic and elastic-plastic 


deformations, respectively 
po, p1 = media densities, when p = po and p = p, 
During unloading of a medium 
p— pr = —Ar(V — V,) [1.4] 
ps, V, correspond to a state to which compression and un- 
loading occur plastically. Quantities p,, V- correspond to a 
state at which unloading begins. Secondary loading occurs 


according to the law of Eq. 1.4, where p < p;, and, where 
pressures are greater, according to a corresponding loading law. 


2 Propagation of Plane Compression Waves 

The solution is effected with the aid of the method stated 
(8); the Lagrangian coordinate system h,t is employed, where 
h is the mass of the substance contained within the examined 
and the initial cross section, and t¢ is the time 


x(h,t) 
h= pdx 
x(0,t) 


The main gasdynamic equations in coordinates /,t have the 
form 
ou . Op 0 Ou 
ot oh oh ot 


According to Eqs. 1.1 to 1.4, there is a linear connection be- 
tween the pressure and the volume, i.e. 


p= —-AV+B [2.2] 


When loading, magnitudes 4 and B are constant within 
each section of the compression diagram; when unloading, 
1 = Ap (constant magnitude), and B = B(h). 

The magnitude B does not depend on ¢; therefore, from 
iiqs. 2.1 and 2.2, we obtain 


Ou Op =0 Ou 1 Op =0 23) 


ot oh Oh A? ot 
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[2.1] 


Fig. 1 “Loading; "unloading and secondary loading 


and the system reduces itself to a homogeneous wave equation 


Oh? 


whose solution is 


p= Fy(h At) + + At) 
u = (1/A)[Fi(h — At) — + 


Discontinuity conditions at coordinates h,t have the form 


= hy?(Vy V2) le UW = hy (Vy = V2) [2.6] 
where h, is the rate of the discontinuity propagation. 

If both states—before discontinuity and during it—lie on a 
single straight line of the linear compression diagram, i.e. 


then the first of the conditions in Eq. 2.6 determines the dis- 
continuity rate 


hy = +A (237) 
the second provides 
po + = py AW [2.8] 


Thus, in passing through a discontinuity (weak or strong) 
one of the expressions p + Au did not experience a jump. 

Let the pressure in the initial cross section of medium h = 0 
grow instantly to p»(px > Pm > ps) and then drop, according 
to the given law p = p(t). Then, two fronts would propagate 
along the medium (Fig. 2). The first, corresponding to a 
pressure jump from po = 0 to p,, moves at velocity Ao; the 
second (which we designate S), corresponding to a pressure 
jump from p, to a maximum pressure (in a given cross sec- 
tion) moves at velocity -A;. 

In region 1 (elastic wave) the solution takes the form 


— In region 2, the medium is unloaded. We seek the flow in 


egion 2, in the form of Eq. 2.5, where A = Ao. 
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0 drop according to the 


Let the pressure in section h = 
linear law 


Pp = — t/8) [2.10] 


where @ is the operating time of the wave. 

Let us obtain functions F; and Fin region 2 in a linear form 
ay by(h Act) 
dy + + Aot) 


Fy(h Apt) 
> Apt) 


boundary conditions 
= —t when h = 0 
Awe = py — Ayu when fA = 


provide four linear equations for determining the four con- 
stants a,, b}, dg and by. a result, we obtain a solution in 
region 2 in the form 


Ds (Ao 


Pm Ao h Ay? A,? t 
Ao (3 ;) 


The particle velocity and pressure, in each section, drop in 
the course of time. At each moment of time, the maximum 


values of these magnitudes lie on the line h = Ait. In ad- 
dition 
1, Ps 1? 6 


At the moment of time ¢ = @, when the pressure in section 
h = 0 drops to value p = 0, a new wave front will proceed at 
velocity Ao. Making calculations, analogous to those carried 
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(Ao = Ai)? 


6 


(Ao — Pmt 


The pressure in each section does not change with time; 
particle velocity is identical in the entire region 3, at each 
moment of time. 

The wave front, reflecting from S, begins to move in the 
opposite direction. A region 4 emerges, in which the solution 
has the form 


= lo 
— + Ae?) (Ag — Ay)? pl 
Ao + Ay)? 0 (Ao + A)? 
Ag = aA, (Ap A,)? 
u=- 2 
(Ao + Pmt (Ao A,)?p, 
Ao + A;)? 6 (Ao + Ay) 
On line h = .A,tin region 4, we have : 
Ag Ay A,)? Poh 
Ag + Ay + Ad) 6 
1, — A, 4 (Ago — 
+ AL) 


Reflecting from plane /, = 0, the wave front will again move 
toward an increase in /. Under repeated reflections, there 


emerge regions 5, 6,7... 
| 


The solution in region 5 is 
+ A)? 6 

(40 — Dn 


The solution in region 6 is 


(Ao — Ay) + A;?) (Ao — Dnl 


(Ao + Ai) 


A;)? 
1 


Ay(Ao + + An)?” 


At the point of intersection (ts,43) of the boundaries be- 
tween regions 6 and 7 and lines h = .A,t, we have 


(Ao — A,)* 


(Ao + Ay)? Ao 1, 
Im 


(Ay — A,)3 
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gut, we obtain a solt 
‘ee 
p™ 
7 
4 “= — — 
Ayo Ay Ag 
; 
i 
Fi 
a 
2 
10° + | ‘led h l 
P= Pn {lt 
=. { 1 0 
: 


Sunilar computations can be made as well for sueceeding 
regions. However, in a series of cases, this has no practical 
significance, inasmuch as the pressure in region 7 is already 
small. Actually, if 4; = 3Ao, then p = 3gpm whenh = hs; 
if Ay = then p = psp, when h = hs. 

The operating time of waves, as follows from the calcula- 
tions given, grows with separation from the initial cross sec- 
tion. If the operating time 6 of pressure p > p, is considered, 
then when the distance is increased there is a growth in @ from 
the beginning; then, when the maximum is reached, it is re- 
duced. The time and place for the maximum to be reached 
depend on the value of p used. 

Fig. 3 shows pressure changes, with time, in some sections 
of the media. 

We have examined the case when the pressure in the initial 
section h = 0 was set in the form of Eq. 2.10. In this sec- 
tion, let the pressure be given as 

p = pr(l — [2.20] 
Let us represent this function in the form 


p=a+att+...+ a,l” 
We seek functions /; and F» in region 2 in the form of 
polynomials of the power n 


ao + — Aol) + anh — Apt)” 
Bo = Bi(h + slo) + + Aot)" 


>| 
| 


Conducting appropriate calculations, we obtain 


—ay( Ag = A,)**} 
Ao*[(Ao + Ay)AT! — (Ao — 1) 


[2.22] 


Having determined all the a; and 8;, we can find the pres- 
sure and the velocity of particles in accordance with Eq. 2.5. 
Thus, a solution in region 2 was obtained. Solutions to sub- 
sequent regions may be found similarly. 


3 Reflection of Compression Waves 


When the compression wave in an elastic-plastic medium is 
reflected from boundary 7 with another elastic-plastic or 
clastic medium, reflected compression or rarefaction waves 
emerge, from a free surface or obstacle, which move along the 
disturbed medium. 

If the state in front of and on the wave front lies on one line 
of the compression diagram, then, as was shown, the front 
velocity ish; + A. 

The conditions on the front are 


pot Aue = + Au if 


— Aw =p — Au if 


If the front velocity coincides with the velocity of the weak 
disturbances behind the front, then functions F.(h + At) = 
3(p — Au), or Fi(h — Act) = 3(p — Au) behind the front re- 
main the same as on the front and in the region before the 
front. 

If the states before and behind the wave front lie on different 
straight lines of the compression diagram, then the pressure 
discontinuity is broken down on two fronts. Behind the first 
front, which is the discontinuity moving with the velocity 
of elastic deformations, a pressure growth occurs in each par- 
ticle to the value reached with the first compression. On the 
second front, appearing as a strong discontinuity which is 
propagated with the velocity of elastic-plastic deformations, 
a pressure jump takes place to the maximum value. 
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-1(h-4,0) 


2(h=2A,0) 


5 10 t 


Fig. 3 


Let us consider the reflection of compression waves, deter- 
mined by Eqs. 2.9 and 2.13, from boundary 7’ with another 
elastic-plastic medium having a greater compressibility. 

The coordinate of the medium boundary is h*. 

All magnitudes, belonging to the second medium, will 
henceforth be designated by a prime. Let us assume that Ao < 
Ao, A’, < A, and p’, > pz. 

At the moment t* = h*/Ao, the elastic wave front reaches 7’. 
The elastic compression wave begins to propagate along the 
second medium and the rarefaction wave along the first 
medium (regions 1 and 3 on Fig. 4). The equation for the 
compression wave front is 


h = + A’), [3.2] 


The equation for the rarefaction wave front (boundary of 
regions | and 3) is 


h = —Aot + 2h* [3.3] 

This line is a continuity, since all states are propagated with 
an identical velocity where there is a linear relationship p = 
p(V). 

The flow in regions 1’ and 3 is determined from the follow- 
ing conditions: It is equal on the boundary of media p and u, 
in both regions; on the compression wave front p — A’ou = 0; 
on the boundary of regions 3 and 1, by virtue of Eq. 2.8, p + 
Agu = ps + Agus = 2p,. 

Hence, the solution in regions 1’ and 3 have the form 


9 
9 - 
‘Ap + 
The relationship of the pressure in reflected wave 3 and in 
incident wave | (reflection coefficient) is 
When the rarefaction wave front encounters the elastic- 
plastic wave front S, regions 4 and 5 form. In region 4, p 
p,, and the function 2 is the same as in region 3. 7 
The solution in region 4 is 


p= DP, “= [(3A0 Ao(-4 + A’) |p. [3.6] 
The particle velocity in region 4 is greater than in regions | 
and 3. 
The expression p + .1,;4 is maintained on the boundary of 


regions 2 and 5, while the expression p — Au is maintained 
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on the boundary of regions 5 and 4. Hence, the solution in 


region 5 is 


2A;(Ao 


+ A;? ml 

(Ay + + Av) 2474, 80 
Dm + AoA,? — 


c= — Ds 


Ay (Ao + A1)(Ao + 
Pmh Ao? + A? Dat 


The pressure and velocity in regions 5 and 2 differ only by 
the constant terms. The pressure is greater and the particle 
velocity smaller in region 2 than it is in region 5. This is 
manifested by the influence of elastic wave reflection from T 
on the flow behind the elastic-plastic wave front. Jumps in 
pressure and velocity occur on the boundary of 2 and 5. 

When there is a reflection of the elastic compression wave, 
determined by the equations in 3.6, from 7’ a new rarefaction 
wave emerges which proceeds in the opposite direction and 
once more reflects from the elastic-plastic wave front S. Wave 
reflection would continue until S approaches 7 . When there 
are numerous reflections, a series of regions form 6, 7, 2’. . ., 
which adjoin front S or boundary 7. Their solutions are 
found as are the solutions in regions 3, 4 and 1’. 

At the moment front S approaches boundary 7’, elastic and 
elastic-plastic compression waves emerge in the second 
medium (regions C and B in Fig. 4), and a rarefaction wave in 
the first medium (region A). It is assumed that the pressure 
on the medium boundary exceeds p’,. 

The solution in region C would be 

p= p's u = p’./A’o [3.8] 

The equation for the rarefaction wave front (discontinuity) 

has the form 


h = —Aot + [(Ao + Ai)/Ai|h* 
We seek the solution in region B in the form 
and in the region A in the form 


The first boundary condition is: When h = h*, the pres- 
sure and the particle velocity in both media are equal 


+ = F, + Fo (F, — [3.9] 


The second boundary condition is: On the elastic-plastic 
compression wave front in the second medium 


+ = = = D's A’yu, [3.10] 


The third boundary condition is: Function F\(h — Aot) 
passes through the rarefaction wave front in the first medium 
without any changes. Inasmuch as F,(h — Aot) does not 
change at all discontinuities which move with velocity Ao, 
F\(h — Aot) must be the same in region A as in region 2. 
Thus, a; and 6; are known. 

Equating the coefficients in Kqs. 3.9 and 3.10 with identical 
powers of t, we obtain another six equations for determining 
the six coefficients de, be, a1, a2, 8; and By. The solutions in 
regions A and B are, correspondingly 


p=F,+F;, u = (1/Ao)(Fi — F2) 
= — 2) 


| 


The pressure in regions A and B drops in the course of time. 

When the rarefaction waves (regions 5, 8, etc.) reflect from 
section h = 0, compression waves form which, in their turn, 
reflect from section h = h* in the form of rarefaction waves. 
In all of the formed regions, the solutions can be obtained 
similarly. 

If the pressure on boundary T does not exceed p’, when the 
elastic-plastic wave front S approaches, then the elastic- 
plastic wave B does not form. As before, function F\(h — 
Aot) in region A is the same as in region 2. Let us designate 
®,; and %, as arbitrary functions which determine u and p in 
region C. It follows, from performance curve conditions, 
bounding region C in the second medium, that 


= — = 0 


In order to determine the four constants in two functions 
= dg + bo(h + Aot) ®, = ao + Bol(h — 


we have two equations, expressing the equality of u and p at 
boundary 7’. Having equated the coefficients where there 
are identical powers of t, we obtain four equations. Thus, the 
problem was solved completely. In Fig. 5, the relationship 
p = p(h) is shown for the moments of time, designated as ¢, 
and ¢, on Fig. 4. 

When h = h*, let the examined medium border on the 
other elastic-plastic medium, in which A’, > Ao, A’; > A; and 
o> Dis 

When the elastic wave front reaches 7, compression waves 
begin to propagate from the boundary on both sides (in the 


first medium, with a velocity of A; and in the second medium 


Fig. 6 
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with a velocity of A’). The pressure and particle velocity in of the boundary for regions 4 and 4’. Then, the second 


regions 3 and 1’ (Fig. 6) are equal and constant and are de- discontinuity ceases to exist and region 5 is formed, bound by 
termined from the conditions on these wave fronts the two characteristics in which medium unloading takes 
; ~ place. The flow in region 5 is determined from the conditions 

p— =O when = + (1! that here the functions F2(h + Aot) is the same as in 4’, 

Ag whereas the function F'\(h — Aol) is the same as in 2. 

A \ When front S reaches the boundary of regions 5 and 4’, a 

p+ Aw = p, (: + ) when 4 = —Ayé + h (: o ) region 6 forms in which the function f(A — Aot) is the same 
Ao as in region 2, whereas the expression p — A,u is maintained 

Hence, the solutions in regions 3 and I’ are on the boundary with region 3. The solution in region 6 is 

At point [+ h (Ao? + 
Ag LA: AoA 
Ag + At Ao + Ai 

the reflected wave front encounters the front S. Two dis- 
velocity Ao, a pressure growth occurs to the e slastic limit for At the moment of time t = A*/A,, the front S reaches the 
each particle. On the second discontinuity, moving at boundaries of the media 7. Further, an elastic and elastic- 
velocity A,, a pressure jump occurs for the elastic limit. The plastic compression wave emerges in the second medium, 
flow in region 4 is determined from the conditions that the ex- while two reflected compression waves emerge in the first 
pression p + Apu is maintained on the boundary with region 2 medium with front discontinuities and moving at velocities of 
while on the boundary with region 4’ the pressure reaches the — Ay and — A, (Fig. 6). The pressure behind the first (weak) 
elastic limit for each particle discontinuity grows to the elastic limit of each particle. Func- 
tion /'\(h — Apt) in region 7 is the same as in 2, whereas the 


2 
P= Pm[l — (lo? — Ar?)h/2-A0°A 8] pressure on the boundary with region A is equal to the elastic 


limit. The solution in region 7 has the form 


(A, 19) 4,6 (A, Ay + Ao)(. A, + A" 
Pm Ao Ao)h* (A, + (A, + Ao) (AY? — Ao’ 2AoA Ps Ay 2Ai(: A’ Ao) ] 
— — 2AoAi(Ar — Ai (Ay + Ao) + Ai) 
= The flow in regions A and B is determined from the con- 
Hence, the solution in region 4 is ditions: The expression p -+ Agu is maintained on the 
(Ap + Ay)2n* boundary with region 7; the pressure and particle velocity 
= Dm on boundary 7’, in both media, are equal, and the expression 
AoAi(Ao — Ai) a p — A‘,u is maintained on the elastic-plastic wave front in 
(Ag SAGAS = the second medium. 
— AO Where the first medium borders with an incompressible 
- medium, the solution in Gave A has the form 
2),* 
‘ 
(Ap + A,)(Ao? + — Ao — Ay Ds 24 Ao 2A (Ao +. 4,) _ A, A)? 2p, 
| (Ao + Ai\(h — h*) 
The medium is unloaded in region 4’; expression p — Aju pm + A,2)9 
is maintained on the boundary with region 3, and the exe [3.17] 
pression p + A,uon the boundary with region 4. Hence, the 
solution in region 4’ is The configuration of the reflected wave system depends 
(Ag? — Ash on the values of magnitudes 6, h*, Ao, Ai, pm and p,. For ex- 
= Ms ample, if p, and @ are relatively large, then the boundary of 
p= Pal at 
: 4. ia Ai(Ao? + AY regions 4’ and 5 intersects directly with the boundary of 
do + Ah (Ao + ADA] — Add 
au p The expressions for the elastic limit of particles in region 2 
(Ao* + o + Ai) 4’ and 6 are different. Therefore, discontinuities arise in the 
D { (Ao? — A12)*h* density, fixed in the coordinate system h,t and a series of new 
| : fields emerge, which were not noted in Fig. 6. Solutions in 
a AeAr\Ad’ + Ar’) regions 7 and A, determined by Eqs. 3.16 and 3.17, were ob- 
(Ao + Ai) Ash (Ae Ae tained for particles in direct proximity to the walls, passing 
(Ag? + 241.0 + Ay) Ds through the discontinuities on the boundary of regions 3 and 6. 
13 14 If p, is substantially less than p,,, the interaction of the 


neglected. Then, region 7 borders on 
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Let the pressure jump become equal to zero at some point plastic wave can be 
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region 2, and the solutions in regions 7 and .{ are determined 
by expressions 3.16 and 3.17, respectively, in which the terms 
containing ps are absent. 

In the foregoing, the cases considered were of the propaga- 
tion and interaction of waves when the pressure had not 
reached the value of p, (Fig. 1); above that point the load- 
ing, unloading and secondary loading lines of the linear dia- 
gram p(V) coincide. If pp > p, in section h = 0 then, follow- 
ing the elastic wave where p = p, and us = ps/Ao, a shock 
wave would be propagated; its front velocity changes from 
V (pm — ps)/(Vs — Vm) when t = Oto V — ps)/(Vs — Vx) 
= A, at that moment of time when the pressure on the front 
declines to py. The propagation of a similar wave and its re- 
action with a fixed obstacle was examined earlier (8). 

The authors would like to express their gratitude to B. A. 
Olisov, Kh. A. Rakhmatulin, and L. I. Sedov for their con- 
sideration of this work. 
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Reviewer’s Comment 


This paper deals with the propagation of dilatational waves 
in an elasto-plastic solid. The nature of the propagation, the 
effect of waves of unloading, and the interaction between com- 
pression and unloading waves are considered, and the reflection 
of plastic waves at the interface between two elasto-plastic 
solids is also treated. The reviewer believes the main weak- 
ness of the paper to lie in the ill-defined physical properties of 
the medium which are postulated. Thus it is assumed that 
the material is elastic up to a given stress and above this stress 
two different linear relations apply, one for loading and the 
other for unloading. The stresses considered, however, are 
hydrostatic stresses, and the relations between pressure and 
volume appropriate to fluid dynamics are employed. This 
is justifiable for solids at pressures very large compared with 
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the shear strength of the material, but plastic deformations 
are then not relevant, since the material behaves in a fluid 
manner for both loading and unloading. The references given 
to other work in this field are rather inadequate, and for this, 
reference should be made to the review paper by Cristescu (1) 
to the book by Goldsmith (2) and to a forthcoming article by 
Craggs (3), 

HerBert Kosky 
Applied Mathematics Div. 

Brown University 
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N AU TOMATIC control system often contains an ampli- 
fier in which the plate circuits are fed with an alternating 
voltage. The mean current in the load is proportional to the 
input voltage, whereas the direction of the current depends on 


the phase; the amplifier is phase-sensitive. Such amplifiers 
are much used, particularly in servo electric drives (2,3) ;! 


Translated from /zvestiia Akademii Nauk Azerbaidzhanskoi 
SSR, Sertia Fiziko-Matematicheskikh Nauk (Bull. Azerbaijan 
Acad. Sci., Physico-mathematical and Technical Sciences Series), 
1959, no. 5, pp. 127-134. Translated by Research Information 
Service, New York. 

' Numbers in parentheses indicate References at end of paper. 
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Transistor Phase-Sensitive Amplifier 


2, North Holland Publishing Co., 1961. 
CH. S. AGALAROV 


they may be made cae and more re Jiable by substituting 
transistors for the vacuum tubes. In some cases such an 
amplifier is required to amplify an a-c signal whose frequency 
is that of the line and also a d-c (slowly varying) signal. The 
sense of the load current must bear a definite relation to the 
phase of the a-e signal and to the sign of the d-c signal, 
as in a servo whose error detector produces an a-c signal and 
which has a d-e operated effector stabilized by d-c feedback. 

Fig. | shows a circuit satisfying these requirements. The 
diodes D, to Ds; prevent the current from flowing in both 
directions in the collector circuits. The load receives a current 
whose direction is controlled by the phase of the a-c signal 
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. Fig. 2. Theoretical circuit of the phase-sensitive amplifier used in an automatic control system 
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and by the sign of the d-c signal, and whose magnitude is de- 
termined by the algebraic sum of the amplitudes. That cur- 
rent is governed by the difference of the reverse currents in the 
diodes, not by the initial reverse currents of the collector 
junctions, when the signals are zero. A proper choice of com- 
ponent values insures that the collector circuit does not be- 
come open while the emitter circuit still carries a current. 
The circuit provides a large power gain. 


06 07 OB U, 


I, ma 


Fig. 3 Output current of the amplifier (Fig. 2) as a function of 


a-c signal 


Calculation of Phase-Sensitive Amplifiers 


The large signal parameters may be used here; Konev (3) 
gives the basic relationships needed for this purpose. A 
specified load current can be provided for if we know Ry», the 
input impedance of the transistor in the common-base con- 
nection, and az, the total static current gain, which is the 
ratio of the total change of collector current (from J.» to 1.) to 


the total change of emitter current (from zero tol.) 
ar = (I. — 1a)/l. 


It is usual to specify the load resistance and the required 
difference in the mean currents, J; — Jz. The line voltage 
required in one section of the transformer winding in order 
to provide that difference in a full-wave circuit is 


Urine 1.110, — + (Ra + R.)/a + {1] 


where 
R. = external resistance in the emitter circuit 
U. = collector-to-base voltage at maximum collector cur- 


rent 


The current is known, as is the reverse voltage on the collector, 
namely = Uiine V 2, so the transistor needed can be selected. 
The maximum power dissipated at the collector of a transistor 
is 


P = 0.16U/Ru [2] 


and must not exceed the limit set for that kind of transistor. 
The signal voltage needed to produce the specified load cur- 
rent is 


U, = 0.7071 (Ri + R.)/ar [3] 
in which /,, is the amplitude of the collector current. The 
input impedance of the stage is we 

Ry = (Ry + ar) [4] 


Fig. 4 Output current of the amplifier (Fig. 2) as a function of d-c signal 
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and the power gain is 


= 0.8laz*Ry/(1 ar) (Ri R.) 


[5] 
These equations apply for a load consisting of an inductor 
whose resistance is 2, and of a capacitance sufficient to remove 
a large fraction of the ripple from the load current. 
= 


Experimental Results 


The four-stage amplifier (Fig. 2) was designed for use in an 
automatic load control for a drill used to make oil boreholes. 
he input stage uses a P6V transistor in the common-collector 


‘onnection. That connection provides an input impedance 
of about 60 kQ, which is fairly high for a transistor a-c ampli- 
fier; a high value was needed to match the output im- 
pedance of the transducer. The second stage uses a P6G 
transistor and is a voltage amplifying stage. The third stage 
‘ontains two P2B transistors and acts as a power preamplifier; 
t also provides the sensitivity to the sign of the d-c signal. 
rhe d-c negative feedback signal is taken to the junction of 
the two secondaries on transformer 7’r.1, which couples the 
second stage to the third. The output stage is a phase-sensi- 
tive power amplifier containing two P4G transistors and is fed 
rom the previous stage via transformer 7r.2.. The amplifier 
eeds the control windings of an amplidyne, which have a re- 
sistance of 18.5 Q and are shunted by capacitors. 

Fig. 3, curve a, shows J; — J,asa function of U’,~, the a-c in- 

yut signal, when the phase of that signal is 0 and 180 deg rela- 
tive to the line voltage. The zero drift at the load is small, 
vhich eliminates any steady-state error in the servo. The out- 
f-balance current corresponding to zero input is 0.1 ma, which ; 
s negligible relative to the nominal control current of the 
implidyne (220 ma, which is produced by an input of 0.13 v). 
The response remains linear up to nearly four times the 
1ominal current (880 ma, which is produced as the maximum 
‘urrent by inputs of 0.6 to 65 v). This insures that the 
transient response of the servo is good. The mutual conduc- 
tance of the amplifier is 1.3 to 1.45 ma per mv; the maximum 
power supplied to the load is 14 w. 

In Fig. 3, curves b and ec show J; — /2 as a function of U’.~ for 
d-e signals UL of land3 v. These signals were applied in the 
sense such that they reduced J; — Js, as would occur in the 
presence of negative feedback. Fig. 4 shows J, — J. as a func- 
tion of U_ for U.~ of 0, 0.3, and 0.4 v. 

This amplifier has been tested in the drill-control system 
under various conditions; it provided the required per- 
formance. 


Se 
Thermal Stabilization 
_ 


Compensation and stabilization are the two methods of 
providing required performance in the presence of temperature 
fluctuations. Compensation for temperature changes is pro- 
duced by including in the circuit of the emitter or base an 
clement whose resistance varies with temperature in such a 
way that the gain of the stage remains constant. Suitable 
elements are thermistors, germanium diodes and_ special 
ceramic resistors having negative temperature coefficients. 
The disadvantages here are that circuits utilizing temperature 
sensitive elements are difficult to analyze analytically and are 
difficult to adjust. 

Stabilization gives satisfactory results for temperature 
rises of 40 to 50 C above normal. Shea’s system is the one 
inost often used; certain specified resistances are included in 
the circuit in order to stabilize the working point (4). Nega- 
tive a-c feedback also gives good results; heavy negative 


a 


1S 


wide temperature range. 


overall gain is reduced. 


for 55 C. 


current is unaffected; 


significant for the purpose. 


eedback insures that the effective parameters are dependent Summary 
in those of the feedback circuit rather than on those of the 
ransistors. This method of stabilization also provides the l 


»dvantage that no difficulty will be encountered if a transistor 
weed be replaced. For example, 100% series negative feed- 
ack may be introduced into a common-collector stage, in 
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Fig. 5 Response curves of amplifier for several temperatures 


which case the voltage gain of the stage becomes about | and 
is almost independent of the transistor parameters. 
put impedance of the stage, subject to the condition R > Ri, 


R; = — @) 


i.e., it depends very little on the parameters of the transistor. 
Tests show that this is so (1); R; is almost constant over a 


Further advantages of thermal stabilization are that it is 
simple to achieve, it permits analytical description, and makes 
transistors interchangeable; its main disadvantage is that the 


All stages in the amplifier were stabilized by negative a-c 
feedback; the working point of the second stage was stabilized 
by providing the bias to the base through a divider. 
input stage had a constant input resistance by virtue of its 
common-collector circuit. Fig. 5 shows J; — Jz = f(U.~) for 
three temperatures. The initial slope (up to U.~ = 0.3 v) in- 
creases with temperature between +10 and +55 C; the in- 
crease (relative to the slope at 10 C) is 7% for 30 C and 28% 
These changes are permissible in an automatic 
control system. The characteristics are the same for all tem- 
peratures for signals larger than 0.3 v. 
the out-of-balance error increases 
slightly with temperature (1.4 ma at 55 C) but remains in- 


A description is given of a phase-sensitive amplifier 
that can be used in any automatic control system in which 
the mean values of two rectified currents are controlled as 
a function of low power a-c and d-e (slowly varying) signals. 
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2 The amplifier has very little zero drift; its slope is 
about 1.4 ma per my, and the maximum power delivered to the 
load is 14 w. 

3 The amplifier is fitted with means of stabilization against 
temperature changes and operates normally to 55 C. 

I am indebted to M. G. Eskin for advice and discussions 
on this work. 
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Reviewer’s Comment 


The output stage is a full-wave power phase detector, quite 
efficient by the necessary use of rectified a-c on the collectors. 
The driver stage makes use of half-wave rectified a-c applied 
to the collectors and the signal is smoothed out into a full sine 
wave by the associated capacitor and transformer primary. 
D-c feedback from possibly a tach generator is mixed with the 
input a-c signal by insertion in the driver base circuit. 

The design of the power stages, although quite an efficient 
method of amplifying at a-c and converting to d-c by phase 
detection, is not new to American scientists, but does indicate 
a Russian application of this circuit using transistors. The 
method of mixing a d-c signal into the base biasing circuit has 
also been used before and is a convenient and simple way to 
mix where extreme linearity is not required. Under circum- 
stances where extreme linearity is required, resistive input 
mixing of the two signals at d-c is usually required. 


The schematics seem correct as printed. It would have 
been better for the author to have spent more time in the ex- 
planation of the curves in Figs. 3, 4, and 5, than in the de- 
sign and stabilization of the first two amplifier stages which 
are quite straightforward in design. 


—WILLIAM QO. Brooks 
Technology Laboratories, Inc. 
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BULLETIN OF THE ACADEMY OF SCIENCES 
USSR, GEOPHYSICS SERIES (Izvestiia Akademii 
Vauk SSSR, Seriia Geofizicheskaia). Published 
by American Geophysical Union, Washington, 


D.C. 


No. 9, Sept. 1959. 


Propagation of Elastic Waves in Layered Media. II, G. 8. Pod’- 
vapol’ski, pp. 913-919. 

Abstract: The paper continues the analysis of an arbitrary 
elementary wave begun elsewhere. An exact expression is ad- 
duced for an axisvmmetric case and a genera! method is given for 
the determination of approximate expressions for the individual 
rapidly changing parts of the displacement field—the seismic 
waves 


Intensity of Earth’s Magnetic Field in the Historical and Geo- 
logical Past, I°. Thellier and O. Thellier, pp. 929-949. 

Abstract: The authors suggest here a method for investigat- 
ing the remanent magnetization of ancient vases and_ bricks. 
These measurements were used for determining the intensity of 
Karth’s magnetic field of times past during which the investi- 
gated objects acquired their magnetization. A consistent and 
steady decrease of Earth’s magnetic moment during the past two 
thousand vears becomes evident; in some regions this decrease 
reaches 65°¢. The methods proposed by the authors take into 
consideration all external influences during the magnetization 
process, and provide in particular, a way for eliminating the 
effect of repeated heating, the so-called ‘temperature cleaning” 
method. 


Problem of an Effective Recombination Coefficient in the 
Ionosphere, A. Bagaryatsky, pp. 966-968. 

Abstract: The article discusses the influence exerted by the 
reaction of an ionomolecular exchange on the magnitude of the 
effective recombination coefficient in the ionosphere. It is 
shown that the concentration [Nz] represents a parameter which 
determines possible changes in the character of the recombination 
with altitude. 


Application of a Formula of Ostrogradsky to the Approximate 
Determination of the Mass of a Body, V. A. Kazinski, pp. 
1003-1004. 


Thermodyramics of Earth’s Mantle, B. N. Zharkov, pp. 10C5- 
1009. 


Influence of Parameters of the Free-Molecular Flow of Raref.ed 
Gas on the Instrument Readings During the Rocket Measure- 
ments of Density of the Higher Layers of the Atmosphere, 
N. G. Belekhova, pp. 1012-1016. 

Introduction: In this paper the relative movement of the con- 
tainer with the instruments in the free-molecular flow of rarefied 
gas, is discussed. Such a container furnished with various instru- 
ments, including the magnetic and ionizing manometers for 
measuring the density of the higher lavers of the atmosphere, 
is a hollow evlinder, the hollow of which is connected with the 
outside medium by means of a number of ‘‘windows,’’ uniformly 
situated along its side surface. In the case discussed, these 
‘windows”’ are narrow slots parallel to the axis of the container. 
Obviously such a construction of the container influences the 
manometer reading to a greater degree than the construction with 
wide windows, described elsewhere. 

Conclusions: 

1 In the ease of a free-molecular flow within the hollow of the 


NOVEMBER 1961 


Information con- 


container of the discussed construction, a state of balance between 
the number of flying in and flying out molecules is established 
quite rapidly (in less than 0.05 see.). 

2 Concentration of the inside molecules within the container 
can be computed with sufficient accuracy by formula 17, and with 
the aid of barometric formula 15 it can be simply related with the 
concentration of molecules .V of the surrounding medium. 

3. The number of molecules getting into the opening of the 
manometer, placed inside of the container of described type, in a 
unit of time, is practically determined by concentration of the 
inside molecules in the hollow of the container and is computed by 


Some Particular Solutions of the Problem of the Propagation of 
a Free Tidal Wave in a Channel of Varying Depth, T. Ya. 
Sekerzh-Zen’kovich, pp. 1041-1044. 

Abstract: The propagation of a free tidal wave in an infinite 
channel of constant width is investigated. It is assumed that 
the depth of the channel varies linearly from one bank to the 
other, the depth being zero at one bank and having some finite 
value at the other. The Coriolis force is taken into account in 
the solution. Particular solutions are obtained in the form of 
Tchebycheff-Laguerre polynomials, and these solutions are 
analyzed. 


formula 27. 


No. 10, Oct. 1959. 


Problem of the Speed of Spreading of Small Cosmic Material 
From the Upper Layers of the Atmosphere Downward, A. A. 
Dmitriev, pp. 1045-1047. 

Abstract: The spreading of very small particles of cosmic origin 
downward, by means of the turbulent mixing, is discussed. A 
scheme with the coefficient of turbulence of kinematic viscosity, 
increasing lineally with height, is used. 


Approximation of Magnetic Anomalies and Reductions by 
Interpolation Polynomials, V. A. Kazinski, pp. 1067-1068. 

Introduction: The polynomial method which has already been 
applied to the reduction and interpretation of subsurface and 
surface gravitation measurements can also be instrumental in 
the testing of models of the elements of 2 magnetic field. 


BULLETIN OF THE ACADEMY OF SCIENCES 
USSR, PHYSICAL SERIES (lIsvestiia Akademii 
Nauk SSSR, Seriia Fizicheskaia). Published by 
Columbia Technical Translations, New York. 


Vol. 23, no. 8, 1959. 


Formation of a High Frequency Discharge at Low Pressures, 
G.N. Zastenker and G. 8. Solntsev, pp. 925-930. 

Introduction: In most investigations of electric discharges the 
fundamental characteristic of initiation of the discharge (break- 
down) is assumed to be the ignition voltage. From the stand- 
point of elucidating the elementary processes involved in the 
development of discharges, knowledge of this parameter alone is 
insufficient. This explains a number of attempts to measure 
another parameter characterizing the onset of discharges, namely, 
the discharge formation time. Investigation of the discharge 
formation time involves observation of the time variation of the 
various parameters characterizing the state of the discharge gap. 

The purpose of the present work was to investigate the forma- 
tion of high frequency (3.3 me) discharges in argon at pressures 
from 0.4 to 15 mm Hg. 

Conclusions: 

1 The developed experimental procedure is suitable for in- 
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vestigating the time variation of different parameters of a high 
frequency discharge (the voltage on the gap, the discharge current, 
and the luminous output of the gas) in the process of its formation. 

2 The formation time of an hf discharge at low pressure (in 
the range of our experimental conditions) varies from 300 to 10 
microsec with variation of the over-voltage from 5 to 60%. 

3 The pressure dependence of the discharge formation time is 
similar to the pressure dependence of the ignition (breakdown) 
voltage. 

4 Theoretical calculation of the duration of the initial stage of 
discharge development, during which the role of space charge is 
minimal, yields satisfactory agreement with the experimental 
data for pressures exceeding 9 mm Hg. 


Measurement of the Charge Concentration in Plasma by Means 
of Super-High Frequency Probes, 8S. M. Levitskii and I. P. 
Shashurin, pp. 938-941. 

Introduction: Within the last decade significant advances have 
been made in the investigation of gas discharges, owing to the 
introduction of a new microwave method, namely, the resonant 
cavity method. Without detracting from the recognized merits 
of this method, we want to mention some of its disadvantages. 

1 To employ the resonant cavity method, the design of the 
discharge tube must be adapted to the design of the cavity, which 
may not always be possible or convenient. 

2 Inits simplest variant the method allows only of measuring 
relatively low electron concentrations (~10%em~*). The method 
yields only the average charge concentration for the volume in 
which the plasma interacts with the microwave field. 

3. Use of the interferometric method which is based on free 
transmission of radio waves through the plasma, largely eliminates 
the first two shortcomings, but in this case, too, the electron con- 
centration obtained is an average value. 

For local measurement of plasma parameters one can hardly 
propose any method not based on the introduction of some type of 
probe into the plasma. In order to retain the advantages of the 
microwave methods and at the same time obtain local values of 
charge concentration, the introduced probe must be an _ shf 
oscillating system with a locally concentrated electric field 
penetrating to a limited extent into the surrounding plasma. 
By immersing such a system into the discharge and observing 
its behavior, one can judge of the charge concentration in its 
vicinity. 

By way of shf probes one can use such oscillating systems as 
various types of symmetrical and asymmetrical dipoles, a two- 
conductor (twin) line, ete. One should, naturally, give pref- 
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erence to systems in which the electric field is more concentrated 
and which are characterized by low emission. 
In our work we used an asymmetrical dipole and a twin line. 


Effect of Disappearance of Negative Ions on Their Concentration 
in the Discharge Column, M. V. Konyukov, pp. 960-963. 

Introduction: The negative component of the plasma in the 
positive column of a discharge in electronegative gases (or in 
media with an appreciable admixture of such gases) consists of 
electrons and negative ions. The latter are apparently re- 
sponsible for some of the distinctive characteristics of such dis- 
charges (anomalous gradients, the formation and subsequent be- 
havior of striations, facilitated contraction of the column, etc.), 
and hence it may be of interest to calculate their concentration. 
The purpose of the present work was to evaluate the influence 
of the decay of negative ions on their concentration. 


Elementary Processes in the Ionization Zone of Corona Dis- 
charges at Atmospheric Pressure, G. N. Aleksandrov, pp. 
976-981. 


Broadening and Shift of Spectral Lines in Gas Discharge Plasma, 
S. L. Mandel’shtam and M. A. Mazing, pp. 1005-1008. 

Introduction: Measurement of the spectral line widths is one 
of the fundamental methods of determining the concentration of 
charged particles in plasma. Investigation of line broadening 
is also of independent physical interest inasmuch as in a plasma 
the radiating atoms are subjected to external influences that can- 
not be duplicated under other conditions: The atoms are under 
the influence of strong, inhomogeneous, rapidly varying fields. 
Theoretical analysis of the interaction of a radiating atom with 
the ambient fields leads to certain relationships between the line 
width y and the charged particle concentration \, i.e., the rela- 
tionships that must be used in experimental determination of \’. 

In the present work, we carried out an experimental investiga- 
tion of broadening under the influence of charged particles of lines 
exhibiting a quadratie Stark effect. 


Measurements of Electric Fields in a High Frequency Low 
Pressure Discharge by Means of an Electron Beam, G. S. 
Solntsev, A. G. Porokhin and N. M. Chistyakova, pp. 1014 
1018. 

Introduction: The techniques of measuring electric fields in 
discharges by observing the deflection of an electron beam passing 
through the discharge has been used in a number of investiga- 
tions of d-c glow discharges at low pressures. In a low pressure 
high frequency (hf) discharge the electric field at each point may 
be represented in the form of a superposition of two fields: 
the hf field E~, produced by the hf oscillator supplying the 
discharge, and a d-c field E- due to distribution of space charges 
in the discharge gap. By shooting an electron beam through 
the hf discharge one can simultaneously determine F~ and E- 
and investigate the distribution of the two fields along the 
discharge axis. 

Conclusions: 

1 The electron beam procedure, which we employed for in- 
vestigating a high frequency, low pressure FE discharge, allows 
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Fig.l. Discharge tube and electron 

beam probe: EG - electron gun, e) & 
C5 - electrodes, DJ - double joint, 
HFG - hf generator. 
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of simultaneous measurement of the hf field and the space charge 
field at different points in the discharge gap. The procedure has 
certain limitations and to obtain meaningful results for hf dis- 
charges certain specific requirements must be met. 

2 In the present work there have been obtained for the first 
time data on the distribution of the electric field, potential and 
space charge in a high frequency discharge in argon at a pressure 
of 10-2 mm Hg for different instants during the hf voltage cvcle. 
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New Data on a Slow Neutron Detector, T. V. Timofeeva and 
S. P. Khormushko, pp. 1275-1277. 


Physical-Chemical Nature of Silver Activator in Zinc Sulfide 
Phosphors, A. A. Cherepnev, pp. 1321-1327. 

Conclusions: 

1 The formation of luminescent ZnS:Ag is apparently con- 
nected with segregation of elemental silver, which is distributed 
through the system in dispersed form (as was earlier suggested 
to the case for copper in ZnS:Cu phosphors). Inasmuch as the 
silver in ZnS is initially present in the form of some compound, 
the process of its precipitation is accompanied by the separation 
of an anionic component; thus there are simultaneously formed 
luminescence centers and traps, characterized by blue emission 
and the low temperature thermostimulation peak. 

2 The process of formation of solid ZnS phosphor proceeds 
primarily in the gaseous phase: The fluxes NaCl and H;BO; 
are comparatively volatile. This determines the condition of 
building up of luminescence centers (emission, trapping and 
quenching centers). Here there occur the effects of precipitation 
in the solid phase, diffusion of the activators, and formation of 
their active states through ‘‘volatile’’ components as well as other 
physical-chemical processes. Owing to the presence of the solid 
system in all the technological stages, the reactions occurring 
in the formation of ZnS phosphors are topochemical (localized). 
Owing to the limited time of the heating operation, equilibrium 
states are not always attained. 

3 Although the final states of the activators silver and copper 
are identical (dispersed elemental particles), the two activators 
exhibit differences as well as similarities. The similarity is 
evinced in the rapid falling off of the absorption curve from the 
maximum value and in weakening of the zine band in the lumines- 
cence spectrum owing to displacement of the zine by the activator 
metal. As regards dispersion and coagulation, however, silver 
behaves differently: Silver is segregated more easily than copper 
and is more stable. The difference between the peaks of the 
blue luminescence bands of zine and silver is difficult to dis- 
tinguish. 

4 The specific characteristics of the two fluxes, NaCl and 
H;BO;, used in our experiments result in differences between 
the processes of formation of the phosphor crystals and between 
the behavior of activators (Cu, Ag and others) in the ZnS phos- 
phor system. Asa flux, NaCl promotes dispersion, chlorination 
and desulfurization; H;BO,; favors vitrification and elimination 
of ZnO by solution and formation of borates. Intensification of 


the phosphor formation process in the presence of H;BO; with 
increasing temperature occurs primarily owing to dissociation of 
the activator sulfide; in the presence of NaCl, as compared with 
H;BOs, owing to the reactivity of the chloride such processes occur 


at lower temperatures. 
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5 The absence of luminescence when ZnS:Ag is heated in an 


atmosphere of hydrogen sulfide apparently occurs owing to 


formation of nonluminescent silver sulfide; in the case of pre- 
formation active luminescent centers appear at 1100 deg. The 
appearance of luminescence is due to the presence of oxidized 
products. These results are connected with the temperature 
dependence of the chemical equilibrium in the systems under 
consideration. 

6 Particularly in view of the work of Kynev on ZnS:Cu 
phosphors, one cannot exclude the possibility that the luminescent 
agent in ZnS:Ag is silver sulfide in colloidal form. On the basis 
of our experimental data and the considerations set forth, how- 
ever, we are inclined to believe that the active agent under the 
conditions of our experiments is dispersed elemental silver. 


Suitability of the Thiosulfate Procedure for Preparation of Zinc 
Sulfide for Synthesis of Phosphors, M. I. Tombak, O. V. Popova, 
A. F. Komar and A. A. Bundel’, pp. 1349-1355. 

Introduction: At present zine sulfide for the production of 
phosphors is prepared by precipitation of the ZnS from a water 
solution of zine sulfate with hydrogen sulfide. 

Guntz proposed a number of other procedures, not involving 
the use of H.S, for preparing zine sulfide; one of these—the 
thiosulfate procedure—was investigated in detail by Grillot. 
Investigation of the thiosulfate procedure is of considerable in- 
terest and value because the phosphors prepared from the product 
differ as regards properties from phosphors prepared from ZnS 
produced by the conventional procedure. Grillot explains this 
by the fact that the reaction product in the thiosulfate process 
is zine disulfide rather than the usual monosulfide. The disulfide 
sulfur is uniformly distributed through the bulk of the product 
and is eliminated only at high temperature, which, in the opinion 
of Grillot, creates conditions favorable for complete reduction of 
all the oxygen compounds present in the batch. 

Thus the distinctive characteristics of ZnS phosphors obtained 
by heating amorphous zine sulfide prepared by the thiosulfate 
procedure are attributed by Grillot to complete absence of oxygen 
in the final phosphor, in contrast to phosphors obtained from the 
ordinary sulfide even after their deoxidation by H2S. 

Inasmuch as some of the points in Grillot’s report appeared 
obscure to us, we undertook an investigation of the reaction of 
sodium thiosulfate with zine sulfate and a study of the phosphors 
prepared of the resultant ‘‘zine disulfide.’’ 


Luminescence of Glasses Containing Cerium, G. O. Karapetyan, 
pp. 1368-1371. 

Introduction: Investigation of glasses activated by cerium is of 
considerable theoretical and practical interest inasmuch as these 
glasses exhibit an exceptional resistance to coloring by ionizing 
radiation. At the same time cerium is a luminescence activator 
and its behavior as such in glasses merits attention. Yet it is 
not feasible in general to extend to glass with cerium the data on 
the absorption and luminescence of allied systems, i.e., crystals 
and solutions containing Ce ions. 

Investigation of the absorption and luminescence spectra of 
glasses containing Ce showed that the element may be present in 
either the trivalent or the tetravalent state. Each valence state 
is characterized by its own distinctive absorption spectrum. 
Change in the type of glass (phosphate, silicate or borate) af- 
fects the level of oxidation of the cerium and also leads to quanti- 
tative changes in the absorption and luminescence spectra. 

It is known that only glasses containing cerium in the trivalent 
state luminesce. Transition of part of the cerium to the tetra- 
valent state results in a substantial decrease in the intensity of 
luminescence. Some authors, including Ginther, report that with 
increase in the cerium concentration its absorption spectrum in 
crystals undergoes great modifications connected with the forma- 
tion of aggregations or clumps of Ce ions. 

The present communication is concerned with a number of 
concentration dependences in glasses activated by cerium. The 
glasses used in the experiment were characterized by high purity: 
The iron content did not exceed 0.005%. Carbon or ammonium 
salts were introduced into the initial batch in order to produce 
reducing conditions. The glasses were founded in quartz cruci- 
bles in an oil fired furnace. 

The absorption spectra were measured on a Beckman spectro- 
photometer using samples 0.1 to 0.5-mm thick and then nor- 
malized to a sample thickness of 1 mm. The absorption spectra 
of cerium were obtained in the form of differential curves: 
The ordinates of the absorption spectrum of the glass without 
cerium were subtracted from the corresponding ordinates of the 
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absorption curve for ‘the glass with cerium. The absolute 
luminescence efficiencies were measured by the sphere procedure 
on the photoelectric setup of V. A. Arkhangel’skaya. The total 
experimental error did not exceed 10% of the measured values. 


Improvement of the Color Rendition of ‘“‘Fluorescent’’ Lamps, 
R.A. Nilender, pp. 1377-1382. 


Vol. 23, no. 12, 1959. 


7 
Nuclear Reactions of Multiply Charged Ions With Carbon and 


Oxygen and Their Influence in Investigation of Coulomb Excita- 
tion, D. G. Alkhazov, A. P. Grinberg, G. M. Gusinskii and I. 
Kh. Lemberg, pp. 1451-1457. 

Introduction: Investigation of Coulomb excitation of high 
lying levels requires increasing the energy of the bombarding 
particles. The use of protons or a-particles in such experiments 
leads to a sharp increase of the y-background (particularly in 
the case of bombardment of low mass number nuclei) that arises 
due to reactions involving the formation of a compound nucleus. 
In principle, the intensity of the background y-lines must de- 
crease greatly if heavy ions are employed as the projectiles. 

However, even in our first attempt to use 25 Mev nitrogen ions 
for excitation of high energy nuclear levels in Sn isotopes we 
encountered a background of y-lines, the intensity of which was 
much greater than the expected intensity of y-rays emitted as a 
result of Coulomb excitation. 

For the purpose of elucidating the origin and character of this 
background y-radiation we investigated the y-spectra of over 40 
different elements, their compounds and isotopes, arising as a 
result of bombardment by C!2, N'4, O'%, Ne?? and Ne? ions. 
The y-radiation was detected by means of a scintillation spec- 
trometer consisting of 40-mm diameter, 40-mm thick Nal( Tl) 
ervstal coupled to an FEU-11 photomultiplier and a 50-channel 
pulse height analyzer. In these experiments we used the same 
“close geometry’’ as in our earlier work on the 4-ra\s emitted as 
a result of Coulomb excitation of nuclear levels. As before, the 
distance from the target to the face of the crystal was 2.7 mm. 
The y-radiations were studied in the region from 0.1 to 2 Mev. 
The spectrometer was calibrated in energy by means of ‘‘stand- 
ard”’ radioactive y-sources. 

By way of bombarding particles, we used the following ions 
accelerated in a cyclotron: 13.6 Mev 11 to 40 Mev 
18.1 Mev O"8, 23.1 to 27.8 Mev Ne” and 25.6 Mev Ne??. 


Inelastic Neutron Scattering by Odd Nuclei, A. M. Korolev, 
pp. 1470-1474. 

Introduction: Excitation of collective nuclear surface levels 
through direct interaction between the nucleus and an inelastically 
scattered neutron is ordinarily described in terms of the Bohr- 
Mottelson model. It is possible, however, that at resonant 
energies the neutron may remain for some time close to the nuclear 
surface, and then a metastable state, which is not a compound 
nucleus state, may form. We have included this mechanism in 
treating inelastic neutron scattering by spherical nuclei by 
Heitler’s theory of damping. 

In the present work we use the methods of our previous analysis 
to study inelastic neutron scattering by nuclei with half-integral 
spin. It was found that because of its complexity the problem 
ean be solved only for weakly deformed nuclei and only taking 
into account excitation of the first two collective levels. 


Magnetic and Quadrupole Moments of Weakly Deformed Nuclei, 
B. D. Konstantinov and A. M. Korolev, pp. 1480-1485. 

Introduction. Bohr’s generalized model of the nucleus is now 
commonly used to explain the experimentally observed magnetic 
and quadrupole moments of nuclei. 

Several authors have used this model to treat these moments 
while taking account of the interaction of an odd nucleon with 
the nuclear surface. Investigations have been published both for 
strong and weak coupling. Nevertheless, there remain some 
questions, for instance, as to the validity of the adiabatic ap- 
proximation in finding the energy levels and the contribution of 
two-phonon states to the magnetic and quadrupole moments. 
Although Choudhury has extended the investigation to three- 
phonon states with intermediate coupling, the analytic depend- 
ence of these quantities on the ground state quantum numbers of 
the nucleus has yet to be given. 

The method used in the present work is essentially the same as 
that of Kerman and Choudhury. We shall investigate the 
magnetic and quadrupole moments of nuclei by approaching 
intermediate coupling from the direction of weak coupling and 
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making use of the results of our preceding work, in which we 
obtained the wave functions and energy levels of odd nuclei. 
Nonadiabatic terms and two-phonon states are included in the 
calculation. Further, the influence of the odd ‘‘extra’’ nucleon 
on the magnitudes of the nuclear moments will be considered. 


Some Characteristics of New Photomultiplier Tubes, A. G. 
Berkovskii, I. Ya. Breido, O. S. Korol’kova and L. G. Leiteizen, 
pp. 1498-1500. 

Introduction: There have been developed two new photo- 
multipliers for use in scintillation spectrometers. The FEU- 
35 photomultiplier has a 25-mm diameter cathode, a maximum 
diameter of 34 mm and an overall length of 108 mm. The tube 
has a 10-pin key-type base and no side leads. 

The other photomultiplier has a cathode diameter of 38 mm, 
a maximum diameter of 48 mm and an overall length of 190 
mm. Like the FEU-35 it has no side leads. 


INDUSTRIAL LABORATORY (Zavodskaia Labo- 
ratoriia). Published by Instrument Society of 


America, Pittsburgh. 


Vol. 25, no. 9, Sept. 1959. 


Determination of Diffusion Coefficients in Polycrystals According 
to Concentration Curves, V. T. Borisov, V. M. Golikov and G. V. 
Scherbedinskii, pp. 1115-1118. 

Introduction: The present article describes a method of de- 
termining the volume D and the boundary D, diffusion coefficients 
according to integral y-radioactivity curves for the residue. This 
method can be applied in those cases where y-radioactive isotopes 
are used in experiments and the concentration on the surface of 
specimens under investigation remains constant during diffusion 
annealing. 


Determination of the Diffusion Coefficient by the Method of 
Residual Gamma Activity, V. Z. Krasil’shchik, I. L. Svetlov and 
M. B. Bronfin, pp. 1118-1121. 


Optical Method of Determining the Total Volume of Foreign 
Matter in Solid Bodies, A. F. Naumov, pp. 1143-1145. 


Method for Determining the Saturated Vapor Pressure of 
Nonvolatile Substances, G. I. Novikov, A. V. Suvorov and A. K. 
Baev, pp. 1148-1150. 

Introduction: In our method the substance being tested is 
evaporated in a closed vessel in the atmosphere of a chemically 
inert gas, which acts as an elastic medium in transmitting the 
vapor pressure to a U-tube manometer placed outside the high 
temperature region. During this process, the vapor of the sub- 
stance being investigated is forced to condense after entering the 
low temperature part of the apparatus. However, the rate of 
condensation is so slow that the equilibrium of evaporation will 
not be disturbed, provided that the transfer of the substance 
resulting from the evaporation in the hot part, and condensation 
in the cold part of the device, is eliminated. 


Determination of the Elastic Moduli of Heat Resistant Materials 
by the Dynamic Method at High Temperatures, A. I. Kovalev 
and I. I. Vishnevskii, pp. 1154-1159. 

Introduction: We developed a method for the high temperature 
measurement of dynamic elastic moduli of heat resistant ma- 
terials, which permitted the simultaneous determination of the 
normal elastic modulus and of the shear modulus at each tem- 
perature point. This method has been used earlier by one of the 
authors of the present article in measurements of elastic properties 
of metallic materials. 


Apparatus for Determining the Thermal Expansion of Films and 
Fibers, A. V. Sidorovich and E. V. Kuvshinskii, pp. 1179-1181. 

Introduction: We have developed an apparatus designed for a 
quantitative linear dilatometric study of polymers in a 107? to 
10-> mm Hg vacuum under dynamic rigidly controlled tem- 
perature-time conditions over temperature ranges from 20-250 
deg. This apparatus makes it possible to heat and to cool samples 
during measurement at constant rates varying from 0.5 deg/hr 
to2deg/min; it maintains a constant temperature during anneal- 
ing to an accuracy of +0.2 deg for several days, etc. While the 
samples are studied dilatometrically, they can also be weighed 
(in order to check and control the invariability of their contents 
while the experiment is being carried out). 
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Automatic Measurement ot the Density of Flowing Liquids, 
B. Z. Votlokhin, pp. 1184-1186. 


Vol. 25, no. 10, Oct. 1959. 


Volumetric Determination of Microgram Quantities of Rare- 
Earth Elements, V. I. Kuznetsov and L. A. Okhanova, pp. 1218 
1221. 

Introduction: After the separation of rare-Earth elements, 
by chromatography or some other method, the absolute quantities 
available may be very small, of the order of a few micrograms. 
It is difficult to determine such small quantities by the normal 
methods. Complexometric methods for the determination of 
rare-Earth elements are described in the literature. For this 
purpose, instead of trilon B, we have used citric acid and arsenazo 
indicator. The method developed for 10-20 y and upward 
of these elements, gave a mean error of 1-2% relative. 


Application of Polarography to the Identification of Plastics, 
V.D. Bezuglyiand V. N. Dmitrieva, pp. 1235-1239. 

Introduction: A number of methods of identifying plasties 
and other high-molecular compounds are described in the litera- 
ture. However, they do not make it possible to determine the 
nature of these materials sufficiently rapidly and accurately. 
We set out to develop a polarographic method for use in deter- 
mining different plastics quantitatively. Many of the monomers 
and also other destruction products, obtained by dry distillation 
of plastics, are reduced on a mercury dropping electrode and are 
‘haracterized by definite half-wave potentials. Some of the 
depolymerization products are not reduced directly, but they may 
he converted into nitro-, nitroso- and halogen derivatives which 
are polarographically active and have definite polarographic 
characteristics. The values of the half-wave potentials of the 
dry distillation products and their derivatives for definite forms 
of plastie may be compared with the characteristics obtained by 
polarographic examination of analysis samples. 


Use of Hard X-Ray Radiation for Determining Residual Stresses, 
M. Ya. Fuks and L. I. Gladkikh, pp. 1248-1252. 


Determination of the Gas Content in Liquids From the Change 
of Electrical Cor-ductance, Kk. Perepelkin, pp. 1265-1268. 

Introduction: Wuring the processes of removal of gas from 
liquids: Flotation, electrolysis, ete., it is necessary to know the 
amount of gas dispersed or dissolved in the liquid phase. A 
number of different methods has been suggested for determining 
the gas contained in liquids; they are, however, not universally 
applicable or not entirely satisfactory, while some of them exclude 
the possibility of remote or automatie measurement. 

For determining the content of gases dispersed in a liquid, a 
method can be used which is based on the dependence of the 
properties of the dispersion system on the content of the dispersed 
phase (for example the electrical conductance), which has been 
considered earlier. 

The known methods of determining the relationship between 
the resistivity of the dispersion systems and the content of the 
electrically nonconductive inclusions for measuring the gas 
content in the liquid were found to be unsuitable for the intended 
purpose since the exact value of the electrical conductivity of the 
medium must be known and also its constancy in time, which in 
some cases cannot be realized. 

We therefore developed a new method, the dilatometric con- 
duction method, which is based on the relationship between the 
gas content in the liquid (and, consequently its electrical resist- 
ance) and the absolute pressure in the system. 


Dynamic Method of Studying the Degassing Process, Bb. I. 
Gorfinkel? and Yu. A. Arkhipov, pp. 1268-1271. 

Introduction: A knowledge of the mechanism of the gassing 
processes is indispensable for solving many problems connected 
with the use of a vacuum. The nature and the speed of such 
processes determine the requirements imposed on the equip- 
ment used, and make possible the correct choice of parameters 
for the production processes. 

The emission of gases by various bodies in vacuum was studied 
by anumber of investigators. However, the necessity of studying 
fast, often highly unsteady processes calls for the use of tech- 
niques different from the usual static methods. The dynamic 
method was found to be the best in the study of unsteady proc- 
esses. 

This paper describes a dynamic method of investigation of the 
iggregate emission of gas. 
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Rapid Method of Calculating the Damping Factor From Oscil- 
lograms, R. 1. Garber and Yu. G. Miller, pp. 1291-1292. 


Initial Amount of Dust in Work With Models of Dust Collecting 
Devices, V. Ic. Maslov and Yu. L. Marshak, pp. 1317-1318. 

Introduction: A simple and reliable method of studying the 
separation of aerosols in models of various dust collecting devices 
has been developed at the All-Union Institute of Calories. Aec- 
cording to this method, the catching surface is covered with a 
viscous liquid (for instance, vaseline). The amount of dust 
settled on the surface can be determined by various physico- 
chemical methods. 


Method of Sampling Liquid Hydrocarbons, A. G. Cheglikov, p. 
1333. 

Introduction: In transferring samples into low pressure gas 
meters, liquid gas, which is throttled after it is let out from the 
storage tank, is evaporated in a small heat exchanger. In this, 
a portion of heavy hydrocarbons can remain in liquid form. In 
further heating, the hydrocarbons evaporate and then partially 
condense in the gas meter. 


Diagram of the sampler, 1) In- 

let valve; 2) packing gland; 

3) inlet valve stem with the shut- 
ing cone; 4) cover; 5) inlet fit- 
ting; 6) frame; 17) check open- 
ing; 8) check valve, 


We developed a device for the sampling of liquid hydrocarbons 
directly from the storage tank. The sample can be kept in the 
sampler for a long time and can be easily transported to the 
place of analysis. In filling the gas meter, the liquid hydrocar- 
bons evaporate, and heavy hydrocarbons which do not evaporate 
at room temperature, remain in the sampler. They are measured 
and taken into account in the processing of analysis data. 


JOURNAL OF APPLIED MATHEMATICS AND 
MECHANICS (Prikladnaia Matematika i Mekha- 
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York, in conjunction with American Society of 
Mechanical Engineers. 


Vol. 24, no. 3, 1960. 


Problems of the General Theory of Plasticity, A. A. Il’iushin, 
pp. 587-603. 

1. The isotropy postulate. The general mathematical theory 
of plasticity is being developed primarily for solid bodies, whose 
material in the undeformed state is isotropic or quasi-isotropic 
(polyerystalline), obeys Hooke’s law in the elastic region and in 
which the formation of plastic deformations is characterized by a 
plasticity condition which coincides with sufficient accuracy with 
the Huber-Mises condition (as, for example, the Tresca and other 
conditions, which replace the Mises surface by polygones which 
are close to it, etc.). For the sake of brevity such bodies will 
be called isotropic in the initial state. 

In approaching the problem of stress-strain relations, which 
determine the basis of the theory, the various plasticity conditions 
of the indicated type may be called approximate representations 
of the Mises condition in all those cases in which the ensuing 
consequences do not differ essentially. 


Stability in the Critical Case of a Zero Root for Systems With 
Time Lag, 8. N. Shimanov, pp. 653-668. 

Abstract: A practical method is given for solving the problem 
on the undisturbed motion of a system with time lag for the 
critical case when one of the roots is zero. 


Equations of Motion for Systems With Non-Ideal Constraints, 
G. K. Pozharitskii, pp. 669-676. 

Abstract: The mechanics of systems subject to constraints 
is based on the assumption that the equations of constraints 
are satisfied exactly during the whole course of motion for arbi- 
trary forces acting on the system and for arbitrary initial con- 
ditions consistent with the constraints. It is clear, however, 
that in an arbitrary mechanical system, the model of which is a 
system with constraints (for example, ideal and holonomic), 
the reactions of constraints arise in consequence of violating the 
latter. In the majority of problems known in the literature, the 
reactions arise in consequence of elastic deformations of bodies 
belonging to the system and those exterior to the latter. 

Since the disturbances of the reactions which arise during the 
motion are very small, their neglect does not cause any essential 
differences between the theory and experiment. This model 
together with a certain hypothesis concerning the properties of 
reactions (their being ideal) permits us to obtain in quite a general 
form the equations of motion for systems with ideal holonomic 
as well as nonholonomic constraints, and also the equations of 
motion for holonomic systems with friction, provided that a 
certain law of friction holds. 

However, due to a frequent appearance of automatically con- 
trolled systems, it seems to us interesting to consider systems 
with automatic devices, the action of which can be interpreted 
as that of a constraint, not belonging to any of the fore-mentioned 
types. 


Explosions Above Surfaces of Liquids, V. M. Kisler, pp. 724- 
734. 

Abstract: The problem of explosions above surfaces of liquids 
is studied under the following assumptions: 

(a) The effect of the external explosion on the motion of the 
free surface of the liquid can be simulated by an appropriate un- 
steady pressure distribution over a time-varying area. 

(b) The motion of the liquid is linearized, which can be justi- 
fied by the large difference between densities of liquid and gas. 

(c) The liquid is considered as incompressible. This as- 
sumption increases in validity as the ratio of the speed of sound 
in the real liquid to the propagation speed of the shock wave over 
the surface increases. 

These assumptions make it possible to reduce the problem at 
hand to a problem of infinitely small surface wave disturbances 
over a heavy incompressible ideal liquid. Apparently, this 
conceptual setting was first used by Lamb in connection with 
problems of long surface waves. In modern times, wave motion 
from this point of view has been studied in some detail by 
Finkelstein. A similar approach was also used by Voit, Cher- 
kesov, and others. 

From the results for explosions above heavy liquids, it is easy 
to derive the motion of free surfaces of weightless fluids by letting 
the gravitational constant g approach zero. This condition 
corresponds to the initial effects of the explosion when the pres- 
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sure forces dominate the gravitational forces. 


Calculation of Flow Past Axisymmetric Bodies With Detached 
Shock Waves Using an Electronic Computing Machine, O. M. 
Belotserkovskii, pp. 745-755. 

Abstract: In solving the title problem, the majority of authors 
have assumed the shape and position of the shock wave, after 
which the inverse problem was solved; a detailed review of 
existing methods is given in the paper. The method of Dorodnit- 
syn permits the direct problem to be solved with the necessary 
degree of accuracy, in its exact formulation. 

Such calculations have been carried out at the Computing 
Center of the Academy of Sciences USSR, on the electronic 
computing machine BESM-1. 

The problem is posed of treating the calculation by a method 
that would be equally suitable for handling both plane and 
axisymmetric bodies of various shapes (smooth, with corners, 
combinations) with detached shock waves, for different values of 
the adiabatic index «(x > 1) and of the free-stream Mach number 
(1 < Me < o@). The plane problem was considered by the 
author, and calculation of flows at 1J~2 = 1 was carried out by 
Chushkin. The calculation scheme and also results of computa- 
tions for certain simple shapes (ellipsoids, spheres and disks) are 
given. 


Flow of an Electrically Conducting Fluid in Tubes of Arbitrary 
Cross Section in the Presence of a Magnetic Field, 8. A. Regirer, 
pp. 790-793. 

Abstract: One of the problems of magnetohydrodynamics is 
the study of the flow of a conducting viscous fluid in a tube located 
in a magnetic field. Its solutions are known for plane and round 
cylindrical tubes. The general problem treated here of the 
steady flow in an infinitely long tube of arbitrary cross section 
reduces to the successive solution of two linear boundary value 
problems. A result of a similar nature was obtained earlier for 
a limited class of externally applied magnetic fields. 


Almost-Periodic Solutions of Nonlinear Systems of Differential 
Equations, V. Kh. Kharasakhal, pp. 838-842. 


Stability of a Gyroscope Having a Vertical Axis of the Outer 
Ring With Dry Friction in the Gimbal Axes Taken Into Account, 
V.V. Krementulo, pp. 843-849. 

Abstract: Rumiantsev has investigated the influence of viscous 
friction on the stability of vertical rotation of a gyroscope on 
gimbals. Other authors have investigated the vertical rotation 
of a gyroscope taking into account dry friction in the suspension. 
The latter problem is further investigated in this paper. 


JOURNAL OF PHYSICAL CHEMISTRY (Zhurnal 
Fizicheskoi Khimii). Published by The Chemical 


Society, London. 

- * 
Vol. 34, no. 5, May 1960. 


Thermodynamics of Equilibria in Similar Reactions, V. A. 
Kireev, pp. 449-455. 

Abstract: The number of chemical compounds used in the 
various branches of the national economy increases so rapidly 
that it is impossible to keep pace with requirements for thermo- 
dynamic data. It is therefore important to be able to calculate 
the equilibrium of one reaction from data for another similar 
reaction. The accuracy of such results is, in general, comparable 
with that of the usual experimental determinations. This in- 
creases the utilization of the existing data and enables the 
equilibrium to be determined for experimentally difficult con- 
ditions such as very high temperatures. 

We shall consider three of the numerous thermodynamic 
methods of determining chemical equilibrium, differing in the 
initial assumptions. 

Comparison of Two Reactions at the Same Temperature 

The second method of calculation is based on the equation 


AH; — AHy AS; — ASy 
RT R 


(11) 


In Ky = In Ky + 


which was obtained by the author from 
Hy ASr 
R 


InK,y = 
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and an analogous expression for In Ky. When Eq. 11 is expressed 
in terms of activities it becomes a rigorous thermodynamic rela- 
tion. 

Comparison of Two Reactions With the Same Values of Equilib- 
rium Constants 

Summary: Basic thermodynamic relationships, which cor- 
relate the parameters of chemical equilibrium in similar (mono- 
typical) chemical reactions, have been given and the more im- 
portant analytical methods of calculating the chemical equilibria 
have been considered. It has been shown that the accuracy of 
the results obtained by such methods is practically the same as 
that of the direct experimental determinations. Approximate 
results may be obtained from very meager initial data, or for the 
less similar (monotypical) reactions. 


Burning Velocity of Ozone-Oxygen Mixtures, N. A. Slavinskaya, 
V. E. Kazakevich, S. A. Kamenetskaya, V. M. Cherednichenko 
and 8. Ya. Pshezhetskii, pp. 462-464. 

Introduction: Kamenetskaya and Pshezhetskii have shown that 
the critical ignition conditions for ozone and its mixtures with 
oxygen correspond to the kinetics of its slow decomposition. 
Values of the lower inflammability limit and the ignition param- 
eters, calculated from Semenov’s theory of thermal ignition 
and from kinetic data were in satisfactory agreement with experi- 
ment. 


€ 
= valve 
voltmeter 
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Fig.1. Circuit for measuring flame-propagation velocity: 

1) photocells; 2) pulse amplifiers; 3) pulse-shaping stages; 
4) charging valve; 5) voltage regulators; 6) valve voltmeter; 
7) milliammeter. 


It is of interest to establish whether there is a similar relation 
between the kinetics of the slow reaction and the burning velocity 
of ozone. Until recently there has been only limited informa- 
tion in the literature about the burning velocity of ozone, and 
we have therefore measured the flame propagation rate for various 
mixtures of ozone and oxygen. The measurements were made in 
a horizontal tube, using a photoelectric method because of the 
low actinic value of the flame. 

Summary: The burning velocity of ozone in ozone-oxygen 
mixtures in a horizontal tube has been measured by a photo- 
electric method. Calculations based on the thermal theory of 
flame propagation have shown that the burning velocity of ozone 
conforms to the kineties of its thermal decomposition. 


Kinetics and Mechanism of Methane Oxidation. II. Kinetics 
of the Accumulation of Intermediate Products, L. V. Karmilova, 
N.S. Enikolopyan and A. B. Nalbandyan, pp. 470-473. 

Introduction: The present paper deals with the kinetics of ac- 
cumulation of formaldehyde and hydrogen peroxide during the 
oxidation of methane in a quartz vessel, treated with HF, and 
operating during the course of a year, with a wide variation of ex- 
perimental conditions 

Summary: 

1 Formation of hydrogen peroxide as well as formaldehyde 
has been shown to occur during the oxidation of methane at 
423-513 deg in a quartz vessel treated with HF. 

2 Formaldehyde is a primary intermediate in the reaction, 
and peroxide is formed by oxidation of the formaldehyde. 

3 The maximum vield of formaldehyde is a linear function 
of the initial pressure of the reaction mixture and of methane, 
and is independent of the oxygen concentration over a wide 
range of pressure. 

4 The activation energy for the formation of formaldehyde 
has been determined (Ecu.0 & 7.8 kcal). The relative yield 
of hydrogen peroxide has been shown to fall with rise in tem- 
perature. 


Aerosols Formed in Reactions, 8. Bogdanov, 
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_E.A. Malysheva, pp.541-548. 


Summary: 

1 An optical siethion has been used to study the formation 
of aerosols during the exposure of gaseous organic substances to 
fast electrons under various conditions, and their subsequent 
changes. 

2 The particle sizes of aerosols formed from methane are 
a few tenths of a micron; the particles carry both positive and 
negative changes (1-11 elementary charges); the majority of 
the charged particles (66%) having 1-3 elementary changes. 
The number of charged particles depends on the degree of dilu- 
tion of the aerosol. 

3 A study has been made of the variation of the weight con- 
centration and the yield G with the amount of energy absorbed 
for aerosols obtained from methane. The vield for ethylene is 
32 times as great as the maximum yield reached with methane. 


Inhibition and the Structure of Inhibitors. III. Inhibited Poly- 
merization of Vinyl Acetate, Z. A. Sinitsyna and Kh. S. Bagdasar’- 
yan, pp. 529-532. 

Introduction: In a previous paper the inhibition of polymeriza- 
tion was used to investigate the reactivity of different aromatic 
compounds toward the methacrylate radical. The inhibited poly- 
merization of vinyl acetate has been investigated in the present 
work. 

Summary: 

1 Each inhibitor molecule terminates two kinetic chains: 
chain regeneration does not occur, and each kinetic chain leads 
to the formation of a single polymer molecule. 

2 Electron-donor and electron-acceptor substituents intensify 
the reactivity of the benzene ring toward the polyvinyl-acetate 
radical. 

3. The value of k:/kp is about 103 times as great for nitro- 
benzene as for the other benzene derivatives. The effect of 
substituents in nitrobenzene is governed by the Hammett rule, 
with p = 0.76. 

4 The reactivity of the polyviny] acetate radical exceeds that 
of the polymethacrylate radical by a factor of 103-104. 

5 It has been shown by the inhibited-polymerization method 
that k,/ko'/? = 0.275 for vinyl acetate at 50 deg. 


Apparatus for Determination of Temperature at Different 
Heights in the Jet of Flame From Gunpowder, P. F. Pokhil, 
i ” Mal’ tsev and L. N. Gal’perin, pp. 539-541. 


Fig.1. Block diagram of apparatus for measuring temperatures 

at different heights of jet of flame from gunpowder: 1) constant- 
pressure bomb; 2) incandescent tungsten lamp; 3) xenon lamp; 

4) optical system; 5) photomultiplier bank; 6) MPO-2 recording 
oscillograph; 7) power pack; 8)-10) FEU-42 photomultipliers; 

F, red filter; F, blue filter; F neutral platinum filter; B biprism; 
Ob objective. 


* Transliterated from the Russian—significance not clear 
(Ed. of Translation). 


Introduction: The apparatus described was designed for measur- 
ing temperatures between 1300 and 3000 deg at different heights 
in the jet of flame from gunpowder at different pressures. The 
temperature was investigated by means of the blue:red ratio 
with a flame emitting a continuous spectrum in the visible region. 
The apparatus enables the absorption to be measured at the height 
of the jet of flame (at definite wave lengths), thus giving the true 
temperature of combustion. 


Simple Apparatus for Production of Pressures up to 200 Atmos- 
pheres, S. I. Sklyarenko, I. V. canis L. B. Belyaeva and 
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Diagram of apparatus for the production of specified pressures 
up to 200 atm (dimensions in mm). 


Introduction: In many kinds of physico-chemical investiga- 
tions it is often necessary to produce high pressures in various 
types of apparatus. Most frequently pressures not exceeding 
100-200 atm are adequate. 

The production of high pressures usually requires quite com- 
plicated apparatus, which may not be available in a nonspecialist 
laboratory. In the absence of such special apparatus, a simple 
device may be employed for producing pressures up to 200 atm 
(and if necessary somewhat higher) in a closed vessel of small 
volume. 


Calculation of Rotational Partition Function for Molecules in a 
Liquid, A. M. Ivseev, p. 548. 

For complex molecules possessing rotational degrees of free- 
dom, the self-consistent field method! allows the asymmetry of 
the force field of the molecule to be taken into account on the 
assumption that the effective field of the other molecules is iso- 
tropic. This circumstance simplifies calculation of the partition 
function. For this calculation the rotation of the molecules in 
the liquid should not be considered free, but hindered, as in intra- 
molecular rotation, in accordance with the hypothesis put forward 
by Evseev and Lebedev.? 

To find the energy spectrum of rotation of the molecule it is 
necessary to solve the Schrédinger equation in a general form 

h? 


for potential energy V. In order to allow for hindered rotation 
we can separate that part of the molecules for which the wave 
equation has the form 


0? 8x7] 


E=90 
h? 


with the condition E > V. 
freely rotating. 
For the rest of the molecules the potential ene ergy 


pressed by 


where ¢ is the angle of rotation of the molecule. 


These molecules 
will perform torsional oscillations or librations. The frequency of 
ose “ills ation is 


Obviously these molecules will be 


can be ex- 


9 J max 
<7 Y2/ 
he fraction of the molecules possessing free rotation is? .. 


/kT 


iEvseev, A. M., Zhur. Fiz. Khim., 1957, vol. 31, p. 2414. 
*Evseev, A. M. and Lebedev, V. P., zbid., 1953, vol. 27, p. 106 
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the fraction performing torsional oscillations is oe 
The energy of the system of interacting molecules can be divided 
into the following parts: 


—€0/kT 


N 
1 Transitional energy of all molecules, yp’, 2mj. 
1 


2 Potential energy of interaction of all molecules, calculated 
by the self-consistent field method, U = You : 

3 Rotational energy of molecules, = Ne~©/*T with 
energy spectrum from e9 to ©. 

4 Energy of torsional oscillation of V2 molecules, 
e—©/kT) with energies of oscillation from 0 to eo. 

Consequently, the Hamiltonian of the system of .V molecules 


= 


can be written in the following way 


N 2 
H; Pi 


9 
2m; 


N 
= 4; + Ejrot > €ilibr 


i=] i=] 


Hence the partition function is equal to 


Vi € Ne 
where | 
‘rot KT = Q, 
0 
j< 


PHYSICS OF METALS AND METALLOGRAPHY 
(Fizika Metallov i Metallovedenie). Published 
by Pergamon Institute, New York. 


Vol. 8, no. 6, 1959. 


Electrical Conductivity of Ferromagnetic Metals in a Radio 
Frequency Field, P. S. Zyrvanov, T. G. Izyumova and G. V. 
Skrotskii, pp. 1-5. 

Introduction: We know that ferromagnetic metals have addi- 
tional specific resistance compared with nonferromagnetic metals. 
This additional resistance is usually explained by dispersion of 
conductance electrons on thermal fluctuations in the degree of 
magnetization. 

In metals which are in conditions of ferromagnetic resonance 
the character of the fluctuations in the degree of magnetization 
can be materially altered. 

The specific resistance p of a metallic ferromagnetic can be 
considered to consist in the following parts 


Piattice + Pym + Pym [1] 


Le eee where Plattice and Pym are, respectively, the resistances caused by 


dispersion of conductance electrons on lattice variations thermo- 
dynamically in equilibrium (phonons) and on variations in 
magnetization (ferromagnons), and py is the additional resistance 
saused by a change in magnetization in a radio frequency field. 
The temperature relation and order of magnitude of plattice are 
well known; pym was calculated by Turov for low temperatures 
on a spin wave model. The temperature relation of pym agrees 
qualitatively with experience. Here an attempt is made to 
construct a qualitative theory for the increase in resistance of 
ferromagnetics in a radio frequeney field. 


Development of Metal Physics and Physics of Metallurgy in 
the Chinese Peoples’ Republic During Ten Years, Chou Chih- 
Hung, Ko Tsun and Tsin Lin-Chao, pp. 16-23. 


Electric Conductivity of Iron at Elevated Temperatures, V. 5S. 
Gumenyuk and V. V. Lebedev, pp. 38-41. 

Abstract: The paper gives data on the study of the temperature 
_ dependence of the electric conductivity of high purity iron, which 
was obtained by vacuum distillation. The curve plotting the 
electric resistivity as a function of temperature shows clearly the 
points of magnetic and phase transition. Comparison is made 
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with the measured values of the electric resistance for armco 
iron and with the results of other authors. A description is given 
of a high temperature (up to 2500C) furnace. 


Concentration Inequalities Experimentally Detected in Solid 
Bodies, V. 1. Arkharov, pp. 51-55. 

Conclusion: 

1 The phenomena of the appearance of concentration non- 
uniformity in allovs investigated by Zav’ialov and Bruk are dif- 
ferent in their nature from the phenomena of internal adsorption. 

2 The phenomena studied by Zavialov and Bruk can be 
explained on the basis of simple thermodynamic conditions in- 
fluencing the form of the constitutional diagram of the alloy; 
however, these explanations are unsuitable for the phenomena of 
internal adsorption because of the completely different nature of 
the latter. 

3. The thickness of the zones of intererystalline internal 
adsorption estimated from the experimental data obtained by 
lifferent methods has a value of 102-108 A; direct auto-radio- 
craphie detection of these zones is extremely difficult and pos- 
sible only in particularly favorable conditions in exceptional 
‘cases. The zones detected in the work of Zav’ialov and Bruk 
have a thickness which exceeds the transverse size of the zones of 
intererystalline internal adsorption. 

4 The remarks of Zav’ialov and Bruk relating to some of the 
ispects of the phenomena of internal adsorption, apart from the 
‘uct that they (the remarks) are made on the basis of incorrect 
dentification of internal adsorption with the phenomena of 
stratification in solid solutions, contain also a series of misunder- 
standings caused by a much too simplified conception of the 
mechanism of internal adsorption. 

5 Despite the shortcomings noted in our article, we consider 
that the work of Zav’ialov and Bruk contains a description of 
phenomena which are important for metallography and correct 
‘xplanations of these phenomena on the basis of the laws repre- 
sented by the constitutional diagrams of alloys. 


Failure of a Martensite Crystal Under the Simultaneous Action 
of External Tensile and Compressive Stresses, A. Mashin, pp. 
95-100. 

Abstract: The present article deals with changes in a martensite 
crystal which is actually a martensite monocrystal, under the 
simultaneous action of external tensile and compressive stresses. 
It is found that the martensite crystals in failing, first form micro- 
cracks, break up, then fracture completely. The interpretation 
of these changes was based on an analysis of stresses acting on the 
martensite crystal during its loading. The crushing of the ervs- 
tals is explained by the fact that the compressive stresses on one 
side of the erystal are not compensated in these cases. The 
problem considered here is solved in connection with problems 
of propagation of elastie and plastic waves in an austenitic- 
martensitic medium. 


Effect of the Original Structure and Temperature of Deforma- 
tion on Phase Transformations During Plastic Formation, 8B. A. 
\paev and Yu. A. Sysuev, pp. 101-106. 

Introduction: Several papers express an opinion on the dif- 
ferences in processes caused by plastic deformation, when equilib- 
rium and nonequilibrium alloys are subjected to plastic deforma- 
tion. It is usual to consider that, in the first case, the deforma- 
tion causes fragmentation of the phases, change in the block struc- 
ture and stressed state of the alloy. In the second case, as well 
as these processes, on reaching a certain degree of deformation 
ore important changes can take place in the alloy, characterized 
by the formation of new phases, which can differ from the original 
phase both in their crystallographic structure and in the chemical 
‘omposition. 

In the light of a number of papers, this classification is more 
than arbitrary, since the results of these papers indicate the 
presence of phase transformations during deformation of the 
alloys in the equilibrium and states close to the equilibrium. 

The opinion held by some authors that the deformation of 
equilibrium systems cannot cause the formation of new phases 
as led, in our opinion, to incorrect conclusions. Taking as a 
starting premise the fact that plastic deformation in equilibrium 
systems only causes those processes mentioned, and observing 
he similar change of the physical properties during the heating 
fa plastically deformed alloy and also the alloy in the hardened 
tate, these authors consider it sufficient proof that during tem- 
ering in the range of the “‘third transformation’’ processes occur 
vhich are caused only by the change in the block structure of the 
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a phase, of its stressed state, by destruction of the coherent bonds 
and by the change in shape and dimension of the carbide particles. 

If the equilibrium state of the alloy is characterized only by the 
ratio and chemical composition of the phases, then, on the basis 
of the existing point of view, the conclusion would be drawn that 
there should be no essential difference in the processes arising 
during the deformation of tempered steel, when the structure of 
the latter is characterized for example, by lamellar and granular 
cementite. In our opinion, this point of view should be made 
more specific, i.e., there should be a more complete determination 
of the equilibrium state of the system, taking into account the 
equilibrium of the form of the phases. 

The present paper deals partly with problems similar to those 
already published and also studies the effect of the deformation 
temperature and the original structure on the character of phase 
transformations and the character of hardening during deforma- 


tion. 


Vol. 9, no. 1, 1960. 


Determination of the Number of Independent Parameters of 
Short-Range Order in Multicomponent Solid Solutions, A. N. 
Men’, pp. 16-19. 

Introduction: When investigating the temperature and con- 
centration dependence of correlation parameters e in different 
coordination spheres for a binary solid solution by the methods 
of the statistical theory of ordering, the problem has to be solved 
for an arbitrary extreme, as the parameters ¢ are not independent. 
In the case of n-component solid solutions with a complex lattice, 
the number of independent parameters of short-range order has 
to be determined in each special case. 

In the present note an attempt is made to construct a general 
scheme for calculating the number of independent parameters 
of short-range order ¢ in any coordination sphere for n-component 
solid solutions with a complex lattice, allowing for the dependence 
of o on the distribution of the atoms with respect to the sub- 
lattices. 


Mechanism of the Formation of Liiders Lines, I.. 8. Ferlikh. pp. 
19-52. 

Abstract: The mechanism of the formation of Liiders lines 
the wrinkles on the surface of soft steel loaded in tension—has not 
so far been established. It is proposed to explain this mechanism 
on the basis of the loss of the rigidity of the thin surface lavers 
under compression. It has been shown that in soft steel test- 
pieces in the process of tension, the conditions necessary for this 
phenomenon do actually oecur. 

Conclusions: 

1 The process of the formation of Liiders lines is made up of a 
number of successive stages. 

2 Liiders lines—these are wrinkles—are the result of the loss 
of stability in the thin surface layer under the action of residual 
compressive stresses. 


Structure of Electrolytically Deposited Metals From the Data of 
Structural Investigations, K. M. Gorbunova, V. V. Bondar’, 
V. P. Moiseev, O. 8S. Popova, I. M. Polukarov and A. A. Sutia- 
gina, pp. 62-67. 

Introduction: As is generally known, the structure of electro- 
depositions determines a number of technical characteristics of 
metallie coatings. With their help it is possible to give desired 
properties (decorative, electric magnetic, mechanical, ete.) to 
the products made from very different metals. The interest 
in the coatings of this type will rise still more when polymeric 
materials commence to be used for construction purposes. 

A wide range of permissible fluctuations in the conditions of 
electro-deposition processing makes it possible to change the 
structure and the properties of coatings of the same metal. At 
the same time the mechanism by which the conditions of electrol- 
ysis affect the strueture and properties of coatings remains un- 
known to a considerable degree and requires further investigations. 


Investigation of Internal Friction in Metal-Ceramic Bodies. 
II. The Ternary System Cu-Ni-Fe, B. Ia. Pines and Den Ge Sen, 
pp. 72-76. 

Abstract: In studying internal friction in metal-ceramic speci- 
mens from a ternary mixture of powders of Cu-Ni-Fe, six peaks 
(maxima) of internal friction were observed, corresponding to 
relaxation effects at the contacts of like and unlike grains. 

Conclusions: 

1 An experimental study of internal friction has been carried 
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out for metal-ceramic specimens from a Cu-Ni-Fe ternary mixture 
of powders, by a torsional oscillation method. 

2 Onthe temperature internal friction curve of specimens from 
the Cu-Ni-Fe ternary mixture of powder, six maxima (peaks) are 
observed, which are evidently due to diffusion transitions of 
atoms along grain boundaries. Of these, three maxima corres- 
pond to processes at the contacts of like grains, and the other three 
to contacts of unlike grains. The position (temperature) of the 
maxima coincides with those found (separately) for the three 
binary systems. 

3 The energy of activation of ‘internal friction processes’’ 
has been found from frequency displacements and the background 
of the internal friction curve. 

4 The values of energy of activation #, corresponding to 
internal friction maxima, agree within the limits of uncertainty 
with those found earlier for the same maxima in the case of 
specimens of the pure components or of binary mixtures of the 


powders. 
5 The energy of activation EH’ found from the ‘‘background”’ 
of the temperature-internal friction curve is appreciably less 


than the values Z determined from the frequency displacement 
of peaks. If the presence of the background corresponds to 
fluctuating transitions of atoms then there may be taking part 
only those atoms in highly distorted region of the crystal (for 
instance on the core of edge dislocations). 

6 All the values of energy of activation for specimens from a 
ternary mixture of powders depend on the duration of preliminary 
annealing; on increasing the time the values of energy of activa- 
tion Z increase, and those for EL’ are somewhat decreased. 


| Reltecoe of Physical Explanation of the Slope of Stress 


Relaxation Curves, G. N. Kolesnikov and A. I. Moiseyev, pp. 
84-86. 


Investigation of Grain Boundary Displacement During Creep, 
V. M. Rozenberg and I. A. Epshtein, pp. 106-113. 

Abstract: A study is made of grain boundary displacement 
during creep in aluminium. It is shown that the extent of grain 
_ boundary displacement is independent of the boundaries’ orienta- 
tion with respect to the elongation axis. Displacements of grain 
boundaries make no contribution to the value of overall elonga- 
tion but they do increase the plasticity of the material. 


RADIO ENGINEERING (Radiotekhnika). Pub- 
lished by Pergamon Institute, New York, in 
conjunction with Massachusetts Institute of 
Technology. 


Vol. 15, no. 7, 1960. 


Dependence of the Voltage Standing Wave Ratio and Losses 
ina Balanced Aerial Switch on the Distance Between the Magne- 
tron and the Dischargers, B. E. Rubinshtein, pp. 24-32. 

Abstract: A study is made of the operation of a balanced aerial 

_ switch with discharge protection for the receiver under reception 

conditions. It is shown that the losses under reception condi- 

- tions are greatly dependent on the distance between the magne- 
tron and dischargers if the dischargers produce a shift in phase. 


Problem of ‘‘Jumps’’ in Electrical Circuits, S. A. Drobov, pp. 
47-60. 
Abstract: A method is elaborated for determining the moments 
at which voltage and current ‘‘jumps’’ arise in circuits with non- 
Tn two-poles and three-poles. It is shown that the solution 
does not require a complete set of differential equations for the 
system being studied. It is sufficient to know certain parameters. 
Several examples are given in the conclusion of circuits widely 
in radio engineering. 


_ Establishing the Frequency at the Output of an Ideal Narrow 
Band Filter With Frequency Phase Modulation, L. I. Iaroslavskii 
and B. I. Iakhinson, pp. 72-83. 

_ Abstract: This article is devoted to a study of transient proc- 
esses, or, more accurately, the establishment of the instantaneous 
frequency and instantaneous amplitude at the output of a narrow 
band system with a frequency phase modulated input signal, 
_ i.e., when sharp and large changes occur in its frequency and 
phase (jumps). Purely frequency and purely phase modulation 
are regarded as special cases of frequency phase modulation. 


_ Investigation Into Transient Processes, Caused by the ‘‘React- 
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ance’? of the Supply Sources, in Class B Amplifiers, A. T. 
Balanov, pp. 112-126. 

Abstract: The results are given of a theoretical and experi- 
mental investigation into periodical recurrent transient processes, 
caused by the reactance of the supply sources in reference to 
large class B-type amplifiers. The effect is shown of feedback 
on the magnitude of nonlinear distortions caused by these tran- 
sient processes. Criteria are given for the design of output filter 
capacitors, rectifiers and blocking capacitors in automatic bias 
circuits. 


Vol. 15, no. 8, 1960. 


New Theory of Active Four-Poles and Its 
tributed Amplifiers, I’. V. Zeliakh, pp. 18-35. 

Abstract: The author introduces new parameters for ative 
four-poles called characteristic (or ‘‘image’’) voltages and cur- 
rents. He formulates the theory of cascade connection for 
matched active four-poles with symmetrical bodies and gives 
formulas for calculating the currents and voltages at the ends of 
the circuit under any load. The theory can be used to analyze 
distributed amplifiers. Formulas are produced for the amplifica- 
tion factors of amplifiers. Account is taken of ‘unmatched 
conditions’’ at both ends of the grid and anode circuits. 


Application to Dis- 


Contribution to the Analysis of Harmonic Frequency Dividers, 
I. Kh. Rizkin, pp. 47-59. 

Abstract: It is shown that the analysis of harmonic frequency 
dividers described by differential equations of a higher order than 
two can in some cases be reduced to the analysis of an equivalent 
divider described by a second-order equation. The method of 
devising the equivalent system is indicated for two common types 
of system. 


Improving the Frequency Phase Characteristics of Selective 
Feedback Amplifiers for Excessive De-Tuning, V. L. Zmudikov, 
pp. 75-80. 

Abstract: In selective, frequency-dependent, feedback ampli- 
fiers, amplification greatly off-tune from the quasi-resonance 
frequency approaches unity and not zero, unlike resonance 
amplifiers. This brings about a deterioration in the selectivity 
and interference resistance of feedback amplifiers. Two methods 
of eliminating this residual amplification by means of subtracting 
circuits are described. 


Interpolation of Functions by Exponential Polynomials and Its 
Application to the Synthesis of Electrical Circuits With Respect 
to the Characteristics, N.S. Kochanov, pp. 87-95. 

Abstract: A study is made of the makeup of exponential poly- 
nomials, with respect to the values of a time function and its 
first derivative at a number of equidistant points of interpolation. 
The use of this method is shown for the synthesis of a four-pole 
forming a sine-squared pulse. 


Spectrum Analysis in Amplitude Phase Modulation, L. FE. 
Kliagin, pp. 96-106. 

Abstract: Exact theoretical formulas are produced which make 
it possible to evaluate the spectra produced in the Kahn-type 
transmitter by amplitude phase modulation. It is shown that 
oscillations with optimum amplitude phase modulation cannot 
be produced and that the second sideband can only be partially 
suppressed; in this respect, the system has no advantages over 
the quadrature modulation system. 


RADIO ENGINEERING AND ELECTRONICS 
(Radiotekhnika i Elektronika). Published by 
Pergamon Institute, New York, in conjunction 
with Massachusetts Institute of Technology. 


Vol. 4, no. 12, 1959. 


Theory of Detecting the Useful Signal in the Background of 
Gaussian Noise and an Arbitrary Number of Interfering Signals 
With Random Amplitudes and Initial Phases, Ia. D. Shirman, 
pp. 9-21. 

Abstract: The detection of a ‘“‘composite’’ signal, when the use- 
ful high frequency signal w4(t) (is known completely, including 
the initial phase or the amplitudes and the initial phase) is not 
only superimposed on a white Gaussian noise but also on m high 
frequency signals ugi(t) (i = 1, 2, ..., m) with a Rayleigh-type 
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amplitude distribution and with an equiprobable initial phase is 
examined. It is shown that the addition of interfering signals 
changes the pulse characteristic of the optimum filter v(t) = 
Aotm4i(to — 8) and the energy of the threshold signals, but it 
does not change the character of the optimum process and the 
form of the detection curve. Recurrent formulas are given for 
the determination of r»4:(t), and the energy utilization factor 
for composite detection is derived. 


Some Methods of Detecting Fluctuating Signals in Two-Channel 
Systems, Iu. B. Chernyak, pp. 22-35. 

Abstract: The problems of detecting signals in a two-channel 
receiving system with linear composition of the amplitudes, 
with maximum amplitude selection and also with detection ap- 
proaching the quadratic composition of amplitudes are examined. 
The comparison of these methods with the quadratic composi- 
tion and with the multiplication of the amplitudes is carried out. 


Accuracy of Measuring the Parameters of Oscillations Distorted 
by Low Intensity Gaussian Noise, B. A. Dubinskii, pp. 36-46. 

Abstract: We examine the problem of determining the ac- 
curacy of the parameters of oscillations during ideal reception in 
KXotel’nikov’s sense, in the presence of Gaussian low intensity 
noise. The derivation is given of expressions for the matrixes 
of the second moments of the parameter errors, which are the 
generalizations of the corresponding expressions given elsewhere 
to the case when the partial differentials of the oscillations with 
respect to the parameters are non-orthogonal functions. As an 
example of using these expressions, the matrix of the second 
moments of the errors of measuring the parameters of a sinusoidal 
pulse signal is calculated. 


Dependence of Irregular Ionospheric Radio Wave Refraction on 
the Zenith Angle, Iu. L. Kokurin, pp. 47-53. 

Abstract: The dependence of the vertical and horizontal ir- 
regular refraction of radio waves (A = 4 m) in the ionosphere on 
the zenith angle is calculated. It is shown that these dependences 
are fundamentally different for various models of the ionospheric 
layer. The necessity of measuring the irregular refraction at 
small zenith angles, in order to establish a model of the layer, is in- 


dicated. 


Application of Chaplygin’s Method for the Solution of Some 
Nonlinear Problems of Radio Engineering, V. I. Kaganov, pp. 
54-58. 

Abstract: The advantages of using Chaplygin’s method for the 
solution of the dynamic equations of systems, containing non- 
linear inertia-free elements are shown. This method enables us 
to find a solution quickly and to evaluate the error. As an ex- 
ample we examine the transient process in an automatic fre- 
quency control system. 


Method of Solving Boundary Problems of Electrodynamics, 
la. N. Fel’d, pp. 74-90. 

Abstract: A method is developed for finding the distribution of 
currents in a system of N metal bodies excited by an incident 
wave or external emf, applied to their surfaces. The solution is 
given in the form of series in special orthonormalized functions 
with finally defined coefficients. The method is illustrated by a 
problem of exciting a system of N finite parallel tubular dipoles 
and infinite strips. 


Spectral Analysis of Quasi-Static Electric Field Space Har- 
monics, G. M. Gershtein and A. V. Khokhlov, pp. 127-137. 

Abstract: A new method proposed earlier for modelling static 
fields, based on the Shockley-Ramo induced-current theorem 
is used for transforming the space harmonics of a quasi-static 
electric field in a periodic structure to time harmonics of induced 
current, with subsequent analysis of the spectrum of the latter. 
Using a low frequency spectrum analyzer the experimental anal- 
ysis of the space-harmonie spectrum of the azimuthal field com- 
ponent of a four-segment structure rotating about a charged probe 
is carried out. Comparisons of the spectra obtained in this way 
with theoretical spectra and with spectra obtained using an elec- 
trolytic tank, show satisfactory agreement. 
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Method for Estimating the Temperature of Hot Stars From Their | 
N III Emission Spectrum, A. A. Nikitin, pp. 461-462. 


System of Programing Concepts, R. I. Podlovchenko, pp. 483- 
485. 

Introduction: The analysis of logical schemes arising in the 
composition of programs has led to the introduction of the con- 
cepts of operator, logical condition and operator sequence. In- 
this paper, one of the possible methods of defining these concepts 
is investigated, and certain relations between operator sequences 
are established on the basis of their functions. 


New Principle for the Construction of a Memory Device, G. G. 
Stetsyura, pp. 486-488. 

Introduction: In this paper, a principle for the construction of a 
parallel-type memory device is set forth. This type of device. 
will allow rapid access to information stored in it when an | 
arbitrary portion of the information itself is given as a guide to- 
the search, rather than its storage location. 


Analyzing the Flow of an Aerated Fluid Assuming a Two- 
Parameter Permeability Characteristic, D. A. Efros and I. F. 
Kuranovy, pp. 497-500. 

Introduction: In hydrodynamic calculations of the flows of an 
aerated fluid, it is the usual procedure to make use of relative 
permeabilities specified in the form of functions Fy (p) and F7(p), 
which are determined in filtration experiments with a mixture of 
liquid and gas introduced from without. The data of recent 
investigations show, however, that, in the stationary filtration of 
a gas-impregnated fluid, the process is described with considerably 
more completeness by two-parameter relations of the form 

kg k 


F4(p,P*) = F\(P*,\) 


Here kg and k; are the permeabilities for the gas and for the 
liquid, respectively, k is the initial permeability, p is the satura- 
tion of the liquid, P* = P/P;, where P; is the pressure at the in- 
flow: \ = T'3*/S(P), here S(P) and I* are the solubilitycoefficient 
and gas factor, both reduced to the pressure P;. 


Motion of an Axially Symmetric Solid in Space About a Stationary 
Point Under the Action of Moments Slowly Varying With Time, 
G. E. Kuzmak, pp. 526-529. 


Speeding up the Limiting Process for Obtaining Roots to Secular 
Equations by the Method of Mayants, B. L. Livshits, pp. 575- 
579. 


Apparatus for High Pressure and High Temperature With a 
Conical Piston, L. F. Vereshchagin, V. A. Galaktionov, A. A. 
Semerchan and V. N. Slesarev, pp. 602-604. 

Introduction: In order to carry out physical investigations, the 
authors have made several different types of apparatus for high 
pressure and high temperature. An apparatus with conical 
pistons is described here. This setup, similar to the recently 
described ‘‘belt’’ setup, was made independently in the develop- 
ment of an idea published earlier and differs from the ‘‘belt’’ in 
that it has a large working volume. 


Fig. 1 shows the essentials of the apparatus. The working 


Fig. 1. Essentials of the apparatus, 
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Tela). 


~ substance G, which is subjected to the high pressure, is placed in 

the cylinder of a matrix C, and is compressed by two pistons 4, 

which have the form of truncated cones. Pyrophyllite sealing 

material D is packed between the conical surfaces of the matrix 
and the plungers. The purpose of this packing is to support, and 
therefore to strengthen, the narrow part of the conical pistons, at 
the cross section where the stress exceeds the elastic limit of the 
piston material. In addition, the pressure gradient along the 
generators of the cone, which arises in the packing when the 
plunger is loaded, provides confinement of the pressure in the 
volume G, not permitting flow of the working substance in the gap 
between the conical surfaces of tne pistons and the matrix. To 
reduce the deformation of the matrix when the pressure is applied, 
and to strengthen it, support is provided by three backing rings B, 
pressed sequentially into one another. Another variation was 
also used, in which support was accomplished by the rods of four 
hydraulic presses, whose axes were situated in a horizontal plane, 
passing through the middle of the working volume and at an angle 
+ of 90 deg to one another. 


_ SOVIET PHYSICS-SOLID STATE (Fizika Tver- 
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The Dispersion and Re-Establishment of Contacts Between 
Microblocks During Plastic Deformation, R. I. Garber, I. A. 
_ Gindin and L. M. Polyakov, pp. 984-990. 

Abstract: In order to explain the reason for the weak strength 
of real solids we must assume that there arises during the deforma- 
tion process an accumulation of dislocations, ruptures and micro- 
cracks leading to large strain concentrations reaching into the 
microregion of the theoretical values of strength. a 


Cross-section 
of valve A 


Fig. 1. Schematic of set-up used to control the porosity of metallic plates, 
1) Mercury manometer, 2) differential oil manometer, 3 and 5) reservoirs 

4) trap for the drying and purification of air entering the set-up, 6) valve 
(junction) for monitoring the samples, The porosity is computed using the 
formula y = /n(Ap,/Ap,)/kt where k is a constant determined by the set- 
up, 7 the time interval, and Ap, and Ap» are the initial and final pressures 
registered by the manometer 2, 


In the case of plastic deformation it is usually considered that 
during all stages of this complex process the complexity is pre- 
served even though crushing and disorientation of the micro- 
blocks are observed. In the present paper we examine the disper- 
sion phenomenon and the processes involved in the re-establish- 
ment of contacts between microblocks which insure the macro- 
scopic integrity of the plastically deformed material. 


The Effect of Heat Treatment on the Photo-EMF in Cuprous 
Oxide, I. D. Kirvalidze, pp. 1043-1045. 

Introduction: Cuprous oxide which is well known for its semi- 
conducting properties exhibits an exceptionally wide scatter in 
conductivity values in dependence on the excess oxygen concen- 
tration. The stoichiometric excess oxygen value depends on the 
previous heat treatment of the specimen. Excess oxygen also 
affects the photo-activity of cuprous oxide. 

The aim of the present investigation was to explain the effect of 
heat treatment in vacuo on the photo-emf in the cuprous oxide 
surface layer. 
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Temperature Dependent Quantum Green’s Functions, Sh. M. 
Kogan, pp. 1074-1084. 


Introduction: The effectiveness of the method of quantum 
Green’s functions in application to the problem of many bodies 
has been very well demonstrated in quite a number of papers. In 
papers by Bonch-Bruevich, it has been shown that quantum 
Green’s functions contain extremely full information about the 
energy spectrum of the quasi-particles in a many-body system 
The connection of Green’s functions with the density matrix has 
been established, and a spectral theorem has been proved on the 
connection of the poles of the Green’s function with the energy 
spectrum of the system. The method of quantum Green’s 
functions has been used to solve problems on plasma oscillations 
and on the screening of an external field in a degenerate electron 
or electron-hole gas for an arbitrary type of anisotropy of the iso- 
energy surfaces. This method has been successfully used for the 
approximate solution of the problem of chemical adsorption on a 
metal. The method of quantum Green’s functions is also used in 
the latest papers of Migdal, Galitskii and Migdal, Galitskii, and 
Belyaev to study the spectrum of the quasi-particles in non-ideal 
Fermi and Bose gases. 

In all of these papers the systems studied are either in the 
ground state or in the immediate neighborhood of the ground 
state. This situation is associated with the use of Green’s func- 
tions defined as averages of 7'-products of field operators over the 
ground state of the system in question. In particular, to study 
the thermodynamic properties of quantum systems it Is necessary 
to generalize the method of Green’s functions to systems that are 
at arbitrary temperatures 7’ > 0. The present paper is devoted 
to the realization of such a generalization 


Theory of Infrared Absorption in Crystals, L.. Ef. Gurevich and Z 
I. Uritskii, pp. 1123-1133. 

Abstract: Absorption of long wave radiation by crystals was 
investigated at frequencies w < w, (w is the photoelectric thresh- 
old) and also in the fundamental absorption region. 

In the first case, phonon absorption was investigated by induc- 
ing with light in crystals virtual excitations which decayed with 
the production of phonons. It was shown that in the continuous 
and discrete absorption regions the most important is absorption 
accompanied by the production of two phonons. Absorption by 
free carriers in a magnetic ficld was also investigated. The ab- 
sorption coefficient was found to oscillate when degeneracy was 
present and also in the case when there was no degeneracy, but 
in a strong magnetie field h Q > 7(Q is the Larmor frequency and 
T is temperature expressed in energy units). 

For the fundamental absorption in a magnetic field, oscillations 
of the absorption coefficient were a function of (w — wo)/Q2 with a 
period equal to 1. In the case of degeneracy the absorption edge 
was displaced, and oscillated owing to the oscillations of the 
Fermi level. 


Ejection of Electrons From Metals by Ions, N. N. Petrov, pp. 
1182-1188. 

Introduction: In order to explain the nature of the kinetic 
ejection of electrons, the present paper treats the excitation of 
bound electrons in a metal. On the basis of the discussions offered, 
recent experimental data on ion-electron emission are evaluated. 

The temperature dependence of potential ejection of electrons, 
the dependence of y on the energy expended in the formation of 
the bombarding ion, and the role of energy losses by the excited 
electrons are also discussed 

Despite the fact that the study of ion-electron emission has 
drawn a great deal of attention in the last few years, there remains 
yet quite a number of unanswered questions. The present paper 
represents an attempt to explain the kinetic ejection of electrons 
by ions, and on the basis of recent data to discuss certain special 
aspects of ion-electron ejection. 


Interaction of Acoustic Oscillations in Ion and Electron-Ion 
Plasmas, P. 8. Zyryanov and FE. G. Skrotskaya, pp. 1196-1200 

Abstract: The role of nonlinear effects in ion and electron-ion 
plasmas in the case of acoustic oscillations is considered. The 
mean free path of Debye phonons is computed, the phonon ther- 
mal conductivity and relaxation time in the phonon system are 
determined. 


SOVIET PHYSICS-TECHNICAL PHYSICS (Zhur- 
nal Tekhnicheskoi Fiziki). Published by Ameri- 
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Vol. 5, no. 2, Aug. 1960. 


Optical Properties of Axially Symmetrical Magnetic Fields With 


Central Location of the Source of Charged Particles, 8S. A. Kuchai, 
pp. 131-140. 
Abstract: In this paper it is shown that axially symmetrical 


fields with central location of the source can be used for isotope 


separation of wide angle and high beams. 


The Mechanism of Space Charge Fluctuations in Quasi-Com- 
pensated Ion Beams, M. V. Nezlin, pp. 154-164. 


>? mm 
2 i 


Constant | 
cure ntam 
| 4 plifier 
Fig. 1. The experimental apparatus, I) lonbeam; 


II)secondary plasma; 1,2) diaphragms; 3) collector. 
The dimensions are in mm. R = 500,000 ohms, 


Abstract: This paper describes a study of the state of compensa- 
tion of the space charge in an intense ion beam passing through a 
rarefied gas in a strong magnetic field. The mechanism is deter- 
mined that causes the fluctuations that produce strong decompo 
sition of the space charge in the beam, so that the electric field in 
the beam reaches some tens of volts per centimeter. 


Method of Investigating the State of Explosion Products by 
Measuring Shock Wave Parameters, ©. A. Tsukhanova, pp. 
219-224. 

Abstract: A method of determining the state of the reaction 
products and the range of the chemical process by the experimen- 
tal investigation of the parameters of the shock waves formed 
during the explosion of gaseous mixtures in closed pipes is investi- 
gated. The results obtained for the mixture 20. + H» are given 
and the heat emitted by the reaction products during nonstation- 
ary flow in pipes is determined 

Summary: On the basis of the 
experimental measurement of 
shock wave velocity, and of the 
pressure and velocity of gas flow 
in front of and behind a shock 
Wave propagating in the explosion 
products of a gaseous mixture in 
a cylindrical chamber, we can 
determine the density distribu- 
tion along the length of the cham- 
ber and compute the value of p/p 
= R7T/u as functions of time. 
These data also permit us to 
judge the degree of completion of 
the chemical process and to com- 
pute the heat transfer from the 
combustion products and the 
sides of the chamber. Processing 
of the experimental results for ex- 
plosions of the mixture 202 + Hy» 
has shown that the equilibrium 
condition is reached after a time 
interval of the order of 10 ~ sec. 


Vol. 5, no. 3, Sept. 1960. 


Study of a Magnetic Trap, K. D. 
Sinel’nikov, V. D. Fedorchenko, 
B. N. Rutkevich, B. M. Chernyi 
and B. G. Safronov, pp. 236-240. 


[See illustration at right.| ator. 
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Abstract: Experiments are described that indicate the aecumu- 
lation of particles in a magnetic trap with a spatially periodic 
magnetic field and the formation in it of a sizable potential well 
for positive ions. 


Experiments on Electrodynamic Acceleration of Plasmas, I. F. 
Kvartskhava, R. D. Meladze and K. V. Suladze, pp. 266-273. 
3 4 
a 2 < 
$ 


Fig. 1. a) Diagram of the coaxial accelerator; b) 
diagram of the induction accelerator; 1) vacuum 
chamber; 2) outer electrode; 3) inner electrode; 
4) porcelain insulator; 5) trigger electrode; 6) dis- 
charge unit; 7) condenser bank; 8) plasma bunch; 
9) vacuum chamber; 10) current-carrying turn. 


Abstract: The results of experiments on electrodynamic ac- 
celeration of plasmas in coaxial and induction accelerators are pre- 
sented. It is found that three different kinds of plasma bunches 
ean be produced in a coaxial accelerator. The maximum bunch 
velocities have been determined. An estimate of the total mass 
of the bunches accelerated during one discharge cycle of the con- 
densers in a coaxial device is made. A qualitative description of 
the described effects is given. 


H 


Fig. 1, Diagram of experimental apparatus: 1) electron gun; 2) collector; 3) reson- 
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Magnetohydrodynamic Effects Observed in Pulsed Contraction 
of a Plasma, I. F. Kvartskhava, K. H. Kervalidze and Yu. S 
Gvaladze, pp. 274-281. 

Abstract: Streak photographs of a pulsed contraction of a 
plasma in strong magnetic fields are discussed. These photo- 
graphs indicate an instability which involves the expulsion of 
plasma formations from the surface of the pinch. The expulsion 
of plasma apparently results from the interaction of radial shock 
waves which are reflected at the pinch axis and the surrounding 
magnetic field. This eruptive instability can be observed in both 
induction and linear pinches Probe measurements of the mag- 
netic field are discussed 


High Frequency Oscillations in a Bounded Plasma, R. A. Demirk- 
hanov, A. K. Gevorkov, A. F. Popov and G. I. Zverev, pp. 282- 
289. 

Abstract: An investigation has been made of oscillations excited cee 
in a bounded plasma in the frequency range from 100-10° to . osc, 
1000-10° eps. It is shown that these plasma oscillations are due 
to oscillations of electrons in the potential well created in the 
plasma. Several oscillations at different frequencies are excited 
simultaneously; these frequencies are lower than the plasma fre- 
quency. The oscillations are accompanied by electromagnetic 
radiation. A possible mechanism for the excitation of these oscil- : ; dyne osc. 
lations is proposed. [See illustration at right.] ' 
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Radiation Energy Losses in a Plasma, V. D. Kirillov, pp. 295- aa 
304. Amplifier 

Abstract: Experiments are described which have been under- weal 
taken to investigate the mechanism for the loss of energy from a c— 
stable plasma pinch which is isolated from the chamber walls. 
The experiments have been carried out in a cylindrical porcelain _| Hetero- Mixer III Lad Amplifier an Amplifier 
system (with electrodes) at discharge currents of tens of kilo- Ill, LF. 

amperes, a longitudinal magnetic field of 24,000 oe and a deute- 
rium pressure of 2-:10-! to 10-2 mm Hg. The duration of the Fig. 2. Block diagram of the UHF receiver with a sensitivity of 10-4 w. 
half-cycle of the discharge current is approximately 500 usec. 
It is found that the flow of charged particles to the walls of the 
chamber is small. Experiments with an ionization chamber, 
vacuum spectrograph and thermal phosphors show that a consid- 
erable portion of the energy loss from the plasma is due to ultra- 
violet radiation from the impurities. 

Conclusion: The results of this group of measurements indicate Y Detector 


that under the present experimental conditions most of the energy 
in the plasma is lost as impurity radiation. Even under unstable 
discharge conditions, in which case the plasma is in contact with : . , 


the walls of the chamber the charged particles carry only a small wll 

part of the Joule heat to the walls. 
The problem of the amount of impurities and the evolution of Fig. 3. Block diagram of the apparatus used for detecting oscillations between 

these impurities in the discharge process in the present experi- 200 and 1000 Meps. 

ments remains open. However, the uniformly low electrical con- 

ductivity in the majority of experiments described in the litera- It would be difficult to achieve any success in heating a deute- 

ture leads us to believe that the level of the impurities in those rium plasma by Joule heating without first removing these im- 

experiments (as well as in the present ones) is still very high. purities. 
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Faced with an almost unsolvable problem, 
Martin engineer McKay Goode helped par- 
lay a unique concept and three years of re- 
search into a spectacular advancement in 
radio communications. © The problem turned 
over to Goode and a team of specialists was to 
design a flexible, interference-free, mobile 
communications system that would reliably 
perform in high-density military operations. 
@® Their answer: a system called RACEP.* 
Here, speech signals are pulse-coded into 
millionth-of-a-second fragments, combined at 
random and transmitted all at onee over the 
same frequency. Only receivers preset for the 


proper code can receive and reconstruct the 
fragments into normal speech. © Result: 
Engineers working in the creative environ- 
ment at Martin in Florida have given the — 
almost saturated rf spectrum a new life. 


*Random Access and Correlation for Extended Performance 


Breakthroughs like this are achieved when 
men with talent are allowed to seek their own 
solutions. At The Martin Company in Or- 
lando, Florida, there is a continuing need for | 
engineers and scientists who can help push 
forward technological frontiers. © For in- — 
formation about your place in one of our six 
major missile and electronics systems, or in- 
one of our future programs, write: C. H. 
Lang, The Martin Company, Orlando 62 


Immediate technical staff openings include: 

COMMUNICATIONS / DIGITAL COMPUTER RESEARCH / HUMAN FACTORS / OPERATIONS RESEARCH / WEAPONS SYSTEMS ANALYSIS / 
GUIDANCE & CONTROL SYSTEMS / HIGH RESOLUTION RADAR / FLUID DYNAMICS / AEROSPACE DYNAMICS / INFORMATION THEORY 
All qualified applicants will receive consideration for employment without regard to race, creed, color or national origin. — 
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Aerojet-General's glass filament winding 
process, AeroROVE, is a simple technique 
applicable wherever outstandin 

reliability and high load carrying capacity 
per unit weight are required. 


Aerojet has built and tested the world's 


| 


largest filament wound rocket chamber. 


i 


This experimental Polaris motor, made of 


i 


a million miles of glass thread, is three 


times the size of any glass case ever 


CORPORATION 
Azusa, California 


THE \ A SUBSIDIARY OF THE 
GENERAL, * °” 


\ tire / GENERAL TIRE AND RUBBER COMPAN 


Engineers, scientists: investigate outstanding opportunities at Aerojet. 
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